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INTRODUZIONE 


Se si sfogliano i piu’ recenti libri di matematica per le superiori 
e li si confrontano con quelli di alcuni anni fa, si nota un’aria di 
novità’. Per i capitoli dedicati alla probabilità’ e alla statistica, si 
tratta addirittura di una vera rivoluzione: i contenuti, il linguag¬ 
gio, le dimensioni, la collocazione e la veste tipografica hanno su¬ 
bito sostanziali modifiche. In effetti, e’ auspicabile che uno stu¬ 
dente della media superiore sappia qualcosa di probabilità’ e di sta¬ 
tistica (qualcosa che non si riduca a frasette del tipo "rapporto tra 
casi favorevoli e casi possibili’’ o "raccolta ed elaborazione di da¬ 
ti) e questi temi, che già’ compaiono nei programmi delle scuole 
elementari e medie inferiori, vengono ora proposti, ad un livello 
adeguato, nel "Programma di matematica per il biennio degli isti¬ 
tuti ad indirizzo scientifico e tecnico". 

I primi due capitoli di questo libro sono dedicati alla probabili¬ 
tà’: "Elementi di calcolo delle probabilità’” e "Problemi di calco¬ 
lo delle probabilità”’. Inizialmente mi ero proposto di mescolarne 
il contenuto, nel tentativo di introdurre elementi di teoria solo quan¬ 
do questi servivano alla soluzione di un certo problema (Problema 
- Costruzione di un modello matematico - Acquisizione degli stru¬ 
menti matematici necessari alla soluzione - Risoluzione - Simu¬ 
lazione). 

Poiché’ ne risultava un discorso piuttosto frammentario, ho de¬ 
ciso di sintetizzare nel primo capitolo le nozioni indispensabili e 
di introdurre i concetti piu’ avanzati nel secondo col procedere de¬ 
gli esercizi. 

II lettore e’ dunque invitato a dare almeno un’occhiata alla pri¬ 
ma parte, non fosse altro che per adeguarsi alla simbologia adot¬ 
tata, e a scegliere quindi fra i 17 esercizi proposti nella seconda quelli 
che ritiene piu’ interessanti. 

I testi dei vari quesiti (spesso presentati come giochi) sono di im- 


mediata comprensione e, nonostante nessun esercizio sia di ovvia 
soluzione, il "bagaglio matematico’’ di uno studente del biennio 
e’ quanto basta per seguirne gli sviluppi (derivate ed integrali com¬ 
paiono molto raramente). 

Molti dei programmi Basic associati agli esercizi operano delle 
simulazioni il cui scopo e’ quello di verificare "sul campo la vali¬ 
dità’ delle conclusioni raggiunte teoricamente (ed e' con un certo 
piacere che si scopre come, alla lunga, le cose vadano secondo le 
nostre previsioni). 

Il terzo capitolo e’ dedicato alla statistica. Si tratta di un capito¬ 
lo piuttosto breve dato che, in questa stessa collana, e’ pubblicato 
"Probabilità’ Statistica e Termodinamica”, che contiene, su que¬ 
sto argomento, molto materiale (e programmi). 

I temi proposti nel terzo capitolo dovrebbero bastare per uno 
studente del biennio. La presentazione di argomenti piu' comples¬ 
si richiederebbe conoscenze matematiche troppo avanzate e, in man¬ 
canza di queste, apparirebbe poco convincente. 

Questo libro ha come destinatari sia studenti che insegnanti; le 
parti piu’ formali sono spesso inserite per completezza e si rivol¬ 
gono principalmente agli insegnanti. Di conseguenza, lo studente 
potrà’ trovare un po’ difficili alcuni paragrafi del primo capitolo: 
senza preoccuparsene eccessivamente, cerchi di assimilare i concetti 
e di familiarizzare con le idee fondamentali. Potrà' rivedere a suo 
tempo le parti piu’ formali, eventualmente con l'aiuto dell’inse¬ 
gnante. 

Sono molto grato a Sergio Casiraghi, che ha curato la versione 
MS DOS dei programmi, e ad Antonella Vitalini che ha paziente- 
mente letto una prima stesura del testo. 

C’e’ una persona che voglio ringraziare in modo particolare: Ni¬ 
colo’ Pintacuda, docente di calcolo delle probabilità' presso l’uni¬ 
versità’ di Pavia. Ho avuto la fortuna di assistere alle sue lezioni 
durante gli anni dell’università’ e gli sono debitore non solo per 
i consigli e suggerimenti relativi al contenuto di questo libro, ma, 
soprattutto, per aver suscitato il mio interesse per uno dei piu’ af¬ 
fascinanti capitoli della matematica. 


Eugenio Rapella 
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CAPITOLO 1 

ELEMENTI DI CALCOLO 
DELLE PROBABILITÀ 


In questo primo capitolo vengono riassunte le nozioni base di 
cui si fara’ ampiamente uso nel secondo, dedicato agli esercizi. 

A chi desidera approfondire questa parte consiglio i testi [14] e 
[15] citati nella bibliografia, si tratta di libri didatticamente prege¬ 
voli e di piacevole lettura. 

Il calcolo delle probabilità’ viene egregiamente trattato in diver¬ 
si libri di testo per le superiori; fra quelli a mia conoscenza merita¬ 
no una particolare segnalazione [4], [8], [19] in bibliografia. 


1.1 LA PROBABILITÀ’ 


Fornire una definizione elementare di probabilità’ non e’ sem¬ 
plice e, di norma, i libri di testo procedono ”per piccoli passi” di¬ 
stinguendo diverse formulazioni. 

In questo paragrafo inizieremo con l’impostazione assiomatica, 
dove la probabilità’ viene definita attraverso alcune proprietà’ di 
cui deve ragionevolmente godere. Successivamente si vedrà’ come, 
a seconda delle informazioni sul fenomeno casuale, si giustifichi¬ 
no la concezione "classica”, quella "frequentata” e quella ”sog- 
gettiva”. 

Giungeremo alla definizione attraverso una serie di osservazioni: 

1 - E’ facile trovare esempi di fatti che possono accadere oppure 
no; li chiameremo EVENTI ALEATORI. Li identificheremo con 
degli insiemi, sottoinsiemi di un insieme piu’ fi vasto che compren¬ 
de tutte le eventualità’; fi sara’ ”1’EVENTO CERTO”. 

Un evento aleatorio sara’ indicato con una lettera maiuscola e, 
quando necessario, ne sara’ data una descrizione, tra parentesi graf- 
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Unione ed Intersezione 
di eventi 


Eventi incompatibili 


fe, con una frase o con un elenco di eventi elementari che caratte¬ 
rizzano la sua realizzazione. 

Esempio 

Consideriamo la situazione in cui viene lanciato un normale da¬ 
do a sei facce. L’insieme fi sara’ l’uscita di uno dei valori da 1 a 6 
fi - 11,2,3,4,5,6} 

Esempi di eventi aleatori sono: 

A = {uscita di un valore minore di 5} = {1,2,3,4} 

B = {uscita di un valore tra 3 e 5) = {3,4,5} 

Come si nota, A c fì e B cr fi. 

2 - Se A e B sono due eventi, indicheremo con: 

A U B l’evento che ha luogo se si verifica A oppure si verifica 
B (o entrambi); 

A fi B l’evento che ha luogo se si verifica sia A che B; 

^(A) l’evento che si verifica quando non si verifica A (l’evento 
contrario). Corrisponde all’insieme complementare di A rispetto 
a fi. 

Riprendendo l’esempio precedente, avremo: 

A U B = {1,2,3,4,5} 

A fi B = {3,4} 

^(A) = {5,6} 

^(B) = {1,2,6} 

3 - L’insieme vuoto 0 rappresenterà’ un EVENTO IMPOSSIBI¬ 
LE. Nel lancio del dado, l’evento 

C = {uscita di un valore maggiore di 6} 
e’ un esempio di evento impossibile: C = 0. 

4 - Diremo INCOMPATIBILI due eventi A e B per i quali e’ AITB 
= 0 (e’ impossibile che si verifichino entrambi). Nel lancio del 
dado, posto X = {1,2} , Y = {4,5,6} si ha XfìY = 0. 

Qualunque sia l’evento A e’ facile verificare che: 

AU A = A; 

ADA = A; 

AU0 = A; 

AD 0 = 0; 

AU^(A) = fi; 

AfW(A) = 0 

5 - Dire che A e’ un sottoinsieme di B (A 9 B) equivarrà’ a dire 
che, quando si verifica A, si verifica anche B, ovvero che ”A im¬ 
plica B”. Nel caso del dado ciò’ avviene, per esempio, per 

A = {2} 
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Probabilità 


Impostazione 

assiomatica 


e 

B = {uscita di un numero pari}. 

6 - Consideriamo una scatola che contiene 100 palline di ugual pe¬ 
so e dimensioni: 99 sono nere, una e’ bianca. 

Viene estratta una pallina, sia 
A = {viene estratta una pallina nera) 

B = {viene estratta una pallina bianca} = ^(A). 

A e B non sono ne’ certi ne’ impossibili, ma l’evento A si verifi¬ 
ca ”piu’ facilmente” dell’evento B. 

Volendo quantificare diversi gradi di incertezza, possiamo pen¬ 
sare di associare ad un evento A un numero reale P(A) - PROBA¬ 
BILITÀ’ DI A - che esprima una misura di "quanto” l’evento A 
possa accadere. Diremo MISURABILI quegli eventi a cui e’ asso¬ 
ciato un valore di probabilità’. 

7 - L’insieme di tutti gli eventi misurabili costituisce una famiglia 
di ja/sottoinsiemi di fi. 

Non e’ indispensabile che jafsia formata da tutti i sottoinsiemi 
di fi (insieme delle parti ^(fi)); e’ comunque ragionevole postula¬ 
re che, se A e B e ja^, anche A U B e ja/; AflBf <tf(A) e 
In altri termini, se agli eventi A e B e’ associato un valore di pro¬ 
babilità’, avra’ senso misurare la probabilità’ che accadano entram¬ 
bi o uno almeno dei due ecc. 

8 - Se richiediamo che: 

• v A e ja/sia P(A) ^ 0 

(la probabilità’ di un evento e’ non negativa) 

• P(fi) = 1 

(la probabilità’ dell’evento certo e’ 1) 

• V A e j/, vBe ja^con A fi B = 0 sia P(A U B) = P(A) 
+ P(B) 

(postulato di additivita’) 

la probabilità’ ha tutti i requisiti per essere una misura sensata del¬ 
l’incertezza di un evento. Dagli assiomi discendono infatti diverse 
proprietà’ che era logico attendersi; eccone alcune: 

P(0) = 0 : l’evento impossibile ha probabilità’ nulla (infatti 
fi = 0 U fi disgiunti quindi P(fi) = P(0) + P(fi) ...) 

A Q B => P(A) ^ P(B) : se un evento ne implica un altro, la 
probabilità’ del primo non supera quella del secondo (infatti, indi- 
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cato con B\A l’evento che si verifica quando si verifica B e non 
si verifica A, si ha: 


B = A U (B \ A) disgiunti quindi 
P(B) = P(A) + P(B\A)^P(A) ) 

Qualunque sia l’evento A e’ 0 £ a c: Q per cui 0 ^ P(A) ^ 
1 : la probabilità’ di un evento e’ un numero compreso tra 0 e 1 

P(^(A)) = 1 - P(A) (infatti fi = A U ^(A) disgiunti quindi 
P(fi) = P(A) + P(^(A)) ...) 


Riassumendo, fissato l’insieme fl di tutte le eventualità’ e una 
famiglia ~o/di sottoinsiemi di Q chiusa per unione (*), intersezione 
(*), passaggio al complementare; la probabilità’ viene definita co¬ 
me una applicazione da jz /in [0,1] tale che: 

- P(fi) = 1 

- v A,B e Jà^con A fi B = 0 

sia P(AUB) = P(A) + P(B) (*) 

[ In generale le proprietà’ (*) vengono estese ad una infinita’ nu¬ 
merabile di eventi, ovvero: 

Se A 1 ,A 2 ,...,A n ,.. e’ una successione di elementi di ja/si ha 
UjA^ jzft HjAie ja/; se gli eventi Aj sono a due a due disgiunti 
(A,flAj = 0 per i^j), la serie £ n P(A n ) converge e la sua som¬ 
ma e’ PdJAJ.] 


La definizione non dice dunque come misurare la probabilità’ 
di un generico evento (ne’ ci si poteva aspettare qualcosa di diver¬ 
so); vai la pena comunque di sottolineare come la probabilità’ di 
un evento possa cambiare al variare delle informazioni in nostro 
possesso sul fenomeno in discussione. 

Se 

A = {La squadra vincitrice del campionato italiano di calcio del 
1975 e’ la Juventus), 
alla domanda 

”Qual e’ la probabilità’ di A” ? 
posso rispondere P(A) = 1, visto che ho appena controllato su 
un’enciclopedia (di cui mi fido). Se la stessa domanda mi fosse stata 
rivolta mezz’ora fa, avrei risposto con un valore minore di uno (dato 
che non mi interesso di calcio) e cosi’, a maggior ragione, se mi 
fosse stata rivolta nel 1974. 


Ritorniamo ad un esempio descritto in precedenza: un normale 
dado a sei facce viene lanciato; ci si chiede la probabilità’ dell’e¬ 
vento A = {uscita di un valore minore di 5}. 

Costruiamo fi e 

fi = {1,2,3,4,5,61 ; j^= &(Q). 

Consideriamo i sei eventi elementari Bj = {i}, i = 1,2,..,6 e no¬ 
tiamo che: 

fi = B,UB 2 UB,UB 4 UB 5 UB 6 

(e’ certo che uscirà’ un numero tra 1 e 6) 


Bj fi Bj = 0 se i j* j 

(l’uscita di un certo numero esclude l’uscita di uno degli altri cinque) 
Per il postulato di additivita’ e’: 


P(fi) = P(B,) + P(B 2 ) + P(B 3 ) + P(B 4 ) + P(B 5 ) + P(B 6 ). 

Si noti come, nella costruzione di fi, abbiamo già’ tenuto conto 
di una informazione: il dado ha sei facce contrassegnate dai valori 
da 1 a 6. Abbiamo un’ulteriore informazione: si tratta di un ”nor- 
male” dado; con questo intendiamo che l’uscita di un certo nume¬ 
ro non e’ ne’ ”piu’ facile” ne’ ”piu’ difficile” dell’uscita di un altro. 

Ciò’ significa che, per un dado di questo tipo che chiameremo 
"regolare”, le probabilità’ degli eventi Bj (i= 1,2,..,6) sono uguali 
tra loro. Se indichiamo con ”x” il loro comune valore, otteniamo: 


P(fi) = x + x + x + x + x + x 


ovvero 6x = 1 e quindi x = 1/6 . 

A questo punto, per calcolare P(A), basterà’ notare che A e’ l’u¬ 
nione dei quattro eventi incompatibili B,, B 2 , B 3 , B 4 . 
Applicando nuovamente il postulato di additivita’, avremo: 


P(A) = (1/6)+ (1/6)+ (1/6)+ (1/6) = 2/3 . 

Sfruttando dunque in modo opportuno le informazioni, si e’ 
giunti alla valutazione effettiva della probabilità’ di un evento. 

In questo esempio, le condizioni che hanno reso possibile que¬ 
sta operazione sono due: 

il fatto che l’insieme fi fosse costituito da un numero finito di eventi 
elementari 

l’aver assegnato ad ogni evento elementare la stessa probabilità’. 
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Impostazione classica 


In situazioni di questo tipo, la probabilità’ di un evento A risul¬ 
ta pari al rapporto tra il numero di elementi di A e il numero di 
elementi di fi. Nell’impostazione classica, la probabilità’ di un even¬ 
to viene definita proprio in questi termini: 

”La probabilità’ di un evento A e’ il rapporto tra il nu¬ 
mero m di casi favorevoli al verificarsi di A e il numero 
n di casi possibili, giudicati ugualmente possibili: 

P(A) = m/n ” 

Questa definizione e’ stata giustamente criticata, sia perche’ può’ 
essere impiegata solo in casi particolari (non sempre sussistono le 
condizioni 1 e 2), sia perche’ non viene indicato quando i ”casi 
possibili” siano da giudicarsi "ugualmente possibili” (ne’ si può’ 
dire che gli eventi elementari sono "ugualmente possibili” se han¬ 
no la stessa probabilità’, dato che la probabilità’ e’ ciò’ che si vuol 
definire). 

Non e’ neppure detto che una stima sensata della probabilità’ 
di un evento sia possibile solo attraverso la definizione classica: negli 
esercizi 8 e 9 del secondo capitolo si valutano delle probabilità’ in 
situazioni in cui i casi favorevoli e quelli possibili non sono in nu¬ 
mero finito. 

Per quanto riguarda la costruzione di fi e si consideri il lan¬ 
cio di una moneta regolare. Posto 

T = {uscita della faccia "testa”); 

C = {uscita della faccia ”croce”), 

costruiamo 

fi = {T,C} e 

= {0»fi»T,C). 

Ragionando come nell’esempio del dado, si avra’ P(T) 
= P(C) = 1/2 . 

Supponiamo ora di voler trattare il problema del lancio di una 
moneta regolare, ma di disporre di un dado regolare. Possiamo uti¬ 
lizzare il dado in luogo della moneta costruendo fi e ja/in modo 
opportuno; basterà’ considerare i due eventi 

A = {uscita di un numero pari ) 

B = {uscita di un numero dispari ) 

e porre 
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fi = {A,B) , ^(ft). 


Impostazione 

frequentista 


Si avra’ P(A) = P(B) = 1/2 ottenendo un modello identico 
al precedente. Con ciò’ si vede come, per uno stesso esperimento 
aleatorio, la coppia (fi. Sì/) possa essere costruita in modi differen¬ 
ti per soddisfare esigenze diverse. 


Consideriamo la seguente situazione: 

”Una scatola contiene palline bianche e nere, indistingui¬ 
bili al tatto, per un totale di 100 palline. Si vuole estrarre 
una pallina; sia 

A = {estrazione di una pallina bianca J. 

Ci proponiamo di stimare P(A). L’impostazione classica non e’ 
di grande aiuto dato che non e’ noto il numero di ”casi favorevo¬ 
li” all’evento A; P(A) può’ essere un qualsiasi valore del tipo k/100 
dove k rappresenta un intero tra 0 e 100. 

Integriamo il testo con una nuova informazione: 

”Dalla scatola sono state in precedenza eseguite 1000 estra¬ 
zioni rimettendo ogni volta nella scatola la pallina estrat¬ 
ta: 950 volte e’ stata estratta una pallina bianca; 50 volte 
una pallina nera”. 

Se chiamiamo '‘FREQUENZA DI A” il rapporto tra il numero 
di prove nelie quali A si e’ verificato e il numero delle prove ese¬ 
guite, potremo dire che la frequenza di A e’ pari a 950/1000 = 
0.95 . 

L’informazione supplementare ha cambiato la situazione e l’i¬ 
potesi che A si verifichi piu’ facilmente di ^(A) ne esce rafforzata. 

Come si vedrà’ nel paragrafo 1.6, sotto certe ipotesi, la frequen¬ 
za può’ essere considerata una buona stima della probabilità’; nel¬ 
l’impostazione frequentista, la probabilità’ di un evento viene ad¬ 
dirittura identificata con la frequenza relativa ad un numero di pro¬ 
ve, eseguite nelle stesse condizioni, ritenuto sufficientemente 
elevato. 

Senza entrare nel merito delle critiche che si possono muovere 
a questa definizione (che, fra l’altro, renderebbe misurabili solo 
quegli eventi per i quali e’ concepibile una ripetizione dell’esperi- 
mento aleatorio in condizioni immutate), si potrà’ accettare la fre¬ 
quenza come stima della probabilità’. In altri termini, nell’esem¬ 
pio della scatola con le 100 palline, non diremo che P(A) e’ 0.95, 
ma che il valore 0.95 e’ quanto di meglio abbiamo per stimare P(A). 
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“Probabilità 

statistiche” 


Impostazione 

soggettiva 


Il calcolo delle probabilità’ e la statistica forniscono gli strumenti 
per valutare i rischi di approssimazioni di questo tipo ed e’ per que¬ 
sto che "PROBABILITÀ’ STATISTICHE” vengono comunemen¬ 
te utilizzate quando si dispone di dati statistici relativi a fenomeni 
che, in linea di massima, si ripresentano in condizioni analoghe (pro¬ 
babilità’ di sopravvivenza, probabilità’ di incidenti, probabilità’ di 
particolari situazioni meteorologiche ecc.). 

Abbiamo visto che la valutazione della probabilità’ di un evento 
applicando la definizione classica o frequentista e’ possibile in am¬ 
biti molto ristretti. Nonostante ciò’, ciascuno di noi e’ chiamato 
ogni giorno a prendere delle decisioni in condizioni di incertezza 
e, piu’ o meno consciamente, a valutare delle probabilità’. 

Supponiamo che il pugile X, detentore del titolo mondiale di una 
certa categoria, stia per incontrare lo sfidante Y; sia 

A = {dopo l’incontro X sara’ riconfermato campione del 
mondo j. 

La probabilità’ di A non si presta ad essere valutata ne’ secondo 
lo schema classico (non c’e’ ragione di pensare che A e ^(A) siano 
equiprobabili), ne’ secondo quello frequentista (i due pugili potreb¬ 
bero non essersi mai incontrati e, anche in caso contrario, le "pro¬ 
ve” non sarebbero ne’ sufficientemente numerose ne’ "eseguite nel¬ 
le stesse condizioni"). 

C’e’ di piu’; l’insieme delle informazioni relative a questa situa¬ 
zione e’ talmente articolato che due persone diverse possono be¬ 
nissimo attribuire ad A diverse valutazioni di probabilità’: c’e’ chi 
scommette su X, c’e’ chi scommette su Y. 

Nell’impostazione soggettiva, la probabilità’ di un evento A viene 
proprio definita ipotizzando una scommessa: 

"La probabilità’ di un evento A e’ quella somma P(A) che: 

1 accetterei preventivamente di pagare per ricevere 1 
lira al verificarsi di A 

e, coerentemente, 

2 accetterei preventivamente di ricevere, impegnandomi 
a pagare 1 lira se A si verifica.” 

La disponibilità’ ad accettare indifferentemente una o l’altra delle 
scommesse, fa’ si’ che V A sia 0 ^ P(A)<1 (se fosse P(A)> 1, la 
prima scommessa risulterebbe comunque svantaggiosa; per 
P(A)<0, lo sarebbe comunque la seconda). 

Vediamo come valutare la probabilità’ di un evento ritenuto certo 
(fi). Con la scommessa (1) sono sicuro di ricevere una lira: per gio¬ 
care sono disposto a pagare una qualsiasi somma x ^ 1. Accettan- 
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Probabilità 

“geometriche” 


Fig. 1 Esempi di pro¬ 
babilità geometriche 


do la scommessa (2) dovro’ certamente pagare una lira; giocare mi 
conviene solo se mi viene preventivamente versata una somma x ^ 1. 
Si ha dunque x ^ 1 e x ^ 1 per cui P(fì) = 1. In modo analogo si può’ 
verificare che P(0) = O. 

L’impostazione soggettiva formalizza la concezione di probabi¬ 
lità’ che si ha comunemente, compresa quella possibilità’ di diver¬ 
se valutazioni da parte di diversi individui che, di fatto, si verifica 
per molti eventi. 

Quest’ultima definizione si può’ applicare praticamente a qual¬ 
siasi evento. 

Concludiamo questo paragrafo con due esempi: 


Es.l Un punto X e’ scelto a caso (nel senso che ogni punto ha 
la stessa probabilità’ di essere scelto) nella regione a (fig. la); qual 
e’ la probabilità’ che X appartenga alla regione /3? 
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I ”casi possibili” sono gli infiniti punti di a , quelli ”favorevo¬ 
li” gli infiniti punti di (3 e la definizione classica non può’ essere 
utilizzata. In un certo senso, l’area di una data regione ne misura 
il numero di punti ed e’ ragionevole porre P[Xe/3] = 

Area(/3)/Area(a). 

Es.2 Si sceglie a caso un numero reale x compreso tra 0 e 1. Qual 
e’ la probabilità’ che sia 1/3 1/2? E che sia x = 3/4 ? 

Ragionando come nell’esempio precedente (vedi fig. lb e le), 
avremo: 


P[l/3^x^ 1/2] 


CD 

AB 



1 

6 


mentre P[x = 3/4] = 0. 

Si noti come l’evento (x = 3/4), pur non essendo rigorosamente 
impossibile, abbia probabilità’ nulla (si usa dire che l’evento e’ 
"quasi impossibile” e il suo complementare "quasi certo”). La cosa 
non deve sorprendere piu’ di tanto; come fa notare F.Speranza 
([18], pag.60), quando si lavora su insiemi infiniti, si presentano 
delle "stranezze” e il calcolo delle probabilità’ non fa eccezione. 

Probabilità’ valutate attraverso modelli geometrici saranno uti¬ 
lizzate negli esercizi 8 e 9 del 2’ capitolo. Naturalmente l’estrazio¬ 
ne di un elemento da un insieme che ne contiene infiniti non e’ at¬ 
tuabile praticamente e le simulazioni su calcolatore forniranno una 
realizzazione approssimata deH’esperimento. Cosi’, se viene gene¬ 
rato un valore pseudocasuale (vedi paragrafo 1.8) dell’intervallo 
(0,1), la probabilità che x sia uguale a 3/4 risulta essere maggiore 
di zero poiché’, nel calcolatore, il numero di casi possibili e’ finito. 


1.2 PROBABILITÀ’ CONDIZIONALE 


Dati due eventi A e B con P(B)>0, chiamiamo PROBABILI¬ 
TÀ’ CONDIZIONALE DI A DATO B (o, piu’ semplicemente, 
PROBABILITÀ’ DI A DATO B) la probabilità’ che si verifichi 
l’evento A nell’ipotesi che B si sia verificato. 

Esempio. Nel lancio di un dado regolare sia: 

A = [1,2,3] 

B = [2,3,4,5,6} 

La probabilità’ di A dato B, che indicheremo con P(A/B), e’ 
la probabilità’ che si verifichi A sapendo che e’ uscito un numero 
maggiore di uno. 

Per effetto della nuova informazione, i casi possibili si riducono 
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a 5 (2,3,4,5,6) di cui solo 2 favorevoli al verificarsi di A (2,3). Si 
ha dunque P(A/B) = 2/5. 

L’ipotesi che B si sia verificato, restringe l’insieme delle even¬ 
tualità’ ai soli elementi di B; contribuiranno alla probabilità’ di A 
solo quegli elementi di A compatibili con B. E’ logico quindi porre 

P(A/B) = P(A fi B)/P(B) 

quantità’ che misura l’incidenza relativa di A in B. 

Per un evento B fissato, la probabilità’ condizionale soddisfa agli 
assiomi della probabilità’; sono di immediata verifica le seguenti 
proprietà’: 

v A e’ P(A/A) = 1 

(E’ certo che A si verifichi, se ci si pone nell’ottica che si sia veri¬ 
ficato) 

Se A fi B = 0 e’ P(A/B)= 0 

(Poiché’ A e’ incompatibile con B, e’ impossibile che A si verifichi 
se supponiamo verificato B) 

V A e’ P(A/fi) = P(A). 

Sia 

H|,H 2 ,..,H n una partizione di Q 

(Hj / 0 v i; H, fi H, = 0 se i / j; H, U H 2 U.. U H n = fi); 
sia A c Q e H, f ja/vi, A e ja/. 

Poiché’: 

A = (AnH 1 )U(AnH 2 )U...U(AnH D ) 

A OH, e’ incompatibile con AflHj per i^j; 
e’ possibile applicare il postulato di additivita’: 

P(A) = P(AD H,) + P(A fi LL) +.. + P(A fi H n ). 

Dalla definizione di probabilità’ condizionale si ha che 
P( A fi H,) = P(A/H j )P(H l ) per cui 

P(A) = P(A/H,)P(H,) + P(A/H 2 )P(H ; ) +.. + P(A/H„)P(H n ) . 

Questa uguaglianza e’ utilissima perche’, molto spesso, risulta 
piu’ semplice il calcolo delle diverse probabilità’ condizionali al se¬ 
condo membro piuttosto che il computo diretto della probabilità’ 
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Fig. 2 Matrice per la 
soluzione del problema 
delle 7 palline 
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1* estrazione 


2“ estrazione 


al primo membro; si confrontino i due metodi di soluzione del se¬ 
guente problema: 

”Una scatola contiene 7 palline: 4 rosse, 3 bianche. Ne 
vengono estratte 2; calcolare la probabilità’ dell’evento 
A = {le palline estratte sono di colore diverso) 

Soluzione 1 

Pensiamo di numerare le palline rosse: RI, R2, R3, R4 e quelle 
bianche: Bl, B2, B3. 

Nella matrice di fig.2 

- le colonne corrispondono ai risultati della prima estrazione 

- le righe corrispondono ai risultati della seconda estrazione. 

Ogni elemento della matrice, ad eccezione di quelli della diago¬ 
nale principale (la stessa pallina non può’ essere estratta due vol¬ 
te), corrisponde ad uno dei casi possibili. Le caselle contrassegna¬ 
te con ”X” identificano i casi favorevoli. 

Procedendo al conteggio, si ottiene: 

”casi possibili” : 42 
”casi favorevoli”: 24 
P(A) = 24/42 = 4/7 

Soluzione 2 

Sia H = {rossa alla 1’ estrazione) 
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K = {bianca alla 1* estrazione) 
e’ HfìK = 0 e HUK=fi; 

la coppia H,K costituisce una partizione di fi per cui 
P(A) = P(A/H)P(H) + P(A/K)P(K) 

P(H) = 4/7. Supponiamo che H si sia verificato, l’urna conter¬ 
rà’ 3 rosse 3 bianche. Perche’ si realizzi A si deve ora estrarre una 
bianca: P(A/H) = l/2. 

P(K) = 3/7. Supposto verificato K, l’urna contiene 4 rosse e 2 
bianche. ”A” si verifica se ora viene estratta una rossa, per cui: 
P(A/K) = 4/6 = 2/3. 

In definitiva 

P(A) = (4/7)(l/2) + (3/7)(2/3) = 4/7 

Il calcolo diviene molto semplice con l’ausilio di un GRAFO 
(fig.3): ad ogni nodo e’ associata una situazione e da esso si stac¬ 
cano i rami delle situazioni ottenibili da quella relativa al nodo. 


Fig. 3 Grafo per la so¬ 
luzione del problema 
delle 7 palline 



( R)R ( R ) B (B)R (B)B 


R = pallina rossa ; B = pallina bianca 


contenuto dell'urna ; (..) estrazione precedente 


Sono sottolineate le situazioni che concorrono al 
verificarsi dell’evento A — { palline di colore diverso } 


P (A) = P[(R)B] + P[(B)R] 
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A fianco di ogni ramo e’ espressa la probabilità’ del verificarsi del¬ 
l’evento ”di arrivo” a condizione che si sia verificato quello ”di 
partenza”. La probabilità’ di un evento X si ottiene moltiplicando 
le probabilità’ dei rami che conducono a situazioni "favorevoli a 
X”, sommando poi i vari contributi. 


1.3 INDIPENDENZA 


Dati due eventi A e B, può’ essere P(A/B)<P(A), P(A/B)>P(A) 
o P(A/B)=P(A). 

L’ultimo caso e’ particolarmente interessante: il fatto di supporre 
B verificato non ha cambiato la probabilità’ di A. Quando 
Eventi indipendenti P(A/B) = P(A), diremo che l’evento A e’ stocasticamente INDI¬ 
PENDENTE dalfevento B. 

Esempio 1 

Vengono lanciati due dadi regolari, uno rosso e uno bianco. Sia 

A = {uscita del 6 sul dado rosso) 

B = {uscita dei 6 sul dado bianco ) 

e’ P(A) = P(B)= 1/6, P(AHB) = 1/36 . 

P(A/B) = P(A fi B)/P(B) = (l/36)/(l/6) = 1 /6 
quindi P(A/B) = P(A) : la probabilità’ che esca 6 sul dado rosso 
non viene alterata dall’informazione che e’ uscito un 6 su quello 
bianco; A e’ indipendente da B. 

Esempio 2 

Lancio di un dado regolare, sia 

A = (1,2) 

B = {2,4,6) . 

E’ P(A)= 1/3 ; P(A/B) = P(AOB)/P(B) = (l/6)/(l/2) =1/3. 

A e’ indipendente da B (anche se i due eventi si riferiscono allo 
stesso esperimento aleatorio). 

Se A e’ indipendente da B, anche B risulta indipendente da A, 
infatti 

P(B/A) = P(B D A)/P(A) = P(A/B)P(B)/P(A) = 

= P(A)P(B)/P(A) = P(B). 

Data la simmetria, basterà’ dire che gli eventi A e B sono indi- 


14 




pendenti (uno dall’altro) e la situazione di indipendenza potrà’ es¬ 
sere caratterizzata dalla proprietà’: 


Legge di una 

variabile 

aleatoria 


A e B sono indipendenti » P(ADB) = P(A)P(B) . 

Una famiglia A, , A 2 A„ di eventi si dira’ COMPLETA¬ 
MENTE INDIPENDENTE se, comunque se ne scelgano k, la pro¬ 
babilità’ che si verifichino tutti e’ pari al prodotto delle probabili¬ 
tà’ di ciascuno di essi. In simboli: 

P(A M n A i? n.. n A ik )=p(A il )P(A i2 )..P(A ik ) 
per qualsiasi scelta degli indici ij . 


1.4 VARIABILI ALEATORIE 


Esempio 1 

Un giocatore lancia un dado regolare e riceve un numero di dol¬ 
lari pari al doppio del punteggio ottenuto. Detta X la somma gua¬ 
dagnata dal giocatore, X può’ assumere uno dei sei valori 2, 4, 6, 
8, 10, 12 e, piu’ precisamente: 

X vale 2 se e’ uscito il numero 1; 

X vale 4 se e’ uscito il numero 2; 
ecc. 

”X” e’ dunque una variabile che "dipende dal caso" nel senso che 
assume certi valori al verificarsi di certi eventi. 

Gli esempi di variabili di questo tipo sono innumerevoli e il loro 
studio e’ di fondamentale importanza nel calcolo delle probabili¬ 
tà’; ci limitiamo a considerare il caso in cui la variabile assume un 
numero finito o una infinita’ numerabile di valori reali: VARIA¬ 
BILI ALEATORIE REALI DISCRETE. 

Fissato l’insieme fi delle eventualità’, sia H,,H 2 ,..,H n ,.. una par¬ 
tizione finita o numerabile costituita da sottoinsiemi misurabili di 
fi. Chiameremo VARIABILE ALEATORIA una funzione X che 
assume il valore reale x ( su H ; . 

Poiché’ gli eventi Hj costituiscono una partizione di fi, posto Pi 
= P(Hj), si avra’ Ejjj = 1. 

Nell’esempio 1 e’ fi = {1,2,3,4,5,6} 

H| = (i) ; Pi = 1/6 ; Xj = 2i per i=l,2,..,6 

L’insieme delle coppie (Xj.Pj) definisce la LEGGE di una varia¬ 
bile aleatoria discreta; variabili aleatorie aventi la stessa legge si 
dicono EQUIDISTRIBUITE. 
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Media di 
una variabile 
aleatoria 


Varianza di 
una variabile 
aleatoria 


Esempio 2 

Lancio un dado regolare; se esce un numero pari vinco 100 lire; 
se esce un numero dispari perdo 100 lire. Sia Y = "somma vin¬ 
ta”; Y e’ una variabile aleatoria di legge: 

Yi +100 -100 

Pi 1/2 1/2 

Esempio 3 

Lancio una moneta regolare; vinco 100 lire se esce ”testa”, per¬ 
do 100 lire se esce ”croce”. Sia Z = "somma vinta”; la legge di Z e’ 

Z; + 100 -100 

Pi 1/2 1/2 

Le variabili Y e Z sono equidistribuite. Come si nota, la legge 
fornisce, molto spesso, un’adeguata descrizione della variabile alea¬ 
toria senza che sia necessaria una costruzione dettagliata di 0 e della 
relativa partizione (il gioco proposto negli esempi 2 e 3 e’ fonda¬ 
mentalmente lo stesso: con probabilità’ 0.5 si vincono 100 lire, con 
probabilità’ 0.5 si perdono 100 lire). 

Se X e’ una variabile aleatoria di legge (x t ,Pj), si definisce ME¬ 
DIA (o SPERANZA MATEMATICA o VALORE ATTESO) di 
X il valore 

E[X] = EjPiXj 

(nel caso di una infinita’ numerabile di x^ la media e’ definita solo 
se la serie converge). 

La media fornisce una stima globale di una variabile aleatoria: 
ogni valore x ( contribuisce al computo di E[X] con un peso pari 
alla probabilità’ di essere assunto. 

Per quanto riguarda gli esempi 1, 2 e 3, avremo: 

E[X] = (l/6)(2 + 4 + 6 + 8+10+12) = 7 

E[Y] = 0 ; E[Z] = 0. 

Variabili aleatorie equidistribuite hanno evidentemente la stessa 
media, ma l’uguaglianza delle medie non implica, in generale, l’u¬ 
guaglianza delle leggi (se, nell’esempio 3, si sostituiscono i valori 
+ 100 e -100 con + 1000 e -1000, si ottiene una variabile aleatoria 
di ugual media e legge diversa). 

1 valori Xj di una variabile aleatoria X possono essere piu’ o me¬ 
no concentrati attorno al valore medio n = E[X]; una misura del 
”grado di dispersione” e’ fornita dalla VARIANZA di X: 
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Giochi 


o*(X) = E[(X - ^) 2 ] = £ (x, - n) 2 Pi . 

ff(X), radice quadrata della varianza, e’ detta SCARTO QUA¬ 
DRATICO MEDIO (o DEVIAZIONE STANDARD) della varia¬ 
bile aleatoria X. 

Esempio 4 

Sia X di legge: 

Xj 4 10 20 

Pi 1/4 1/2 1/4. 

Calcolare E[X], <r 2 (X). 

Avremo: 

fi = E[X] = 4(1/4)+10(1/2)+ 20(1/4) = 11 

a 2 (X) = (4-1 l) 2 (l/4) + (10-1 l) 2 (l/2) + (20-1 l) 2 (l/4) = 

= 12.25 + 0.5 + 20.25 = 33 
a(X) = V33 = 5.74 

Se, come negli esempi 1, 2 e 3, i valori di una variabile aleatoria 
X rappresentano i guadagni (positivi o negativi) associati al realiz¬ 
zarsi di certi eventi, E[X] rappresenta il "guadagno medio teori¬ 
co” ovvero la somma che un giocatore può’ aspettarsi di guada¬ 
gnare mediamente per ogni prova eseguita (se le prove vengono ese¬ 
guite nelle stesse condizioni e il loro numero e’ sufficientemente 
elevato). 

Il gioco e’ detto: 

VANTAGGIOSO se E[X] > 0 
SVANTAGGIOSO se E[X] < 0 
EQUO se E[X] = 0 

Il gioco descritto nell’esempio 1 e’ vantaggioso (in effetti si può’ 
soltanto vincere), quello delTesempio 2 e’ equo. 

Se due variabili aleatorie hanno la stessa media, quella con va¬ 
rianza maggiore e’ ”piu’ aleatoria”: lo scarto dal valore atteso e’, 
mediamente, piu’ grande. 

Dovendo scegliere tra due investimenti 

Invest.l X: guadagno -100 +200 
con prob. 1/2 1/2 

lnvest.2 Y: guadagno -100 +14900 

con prob 99/100 1/100 


noterò’ che E[X] = E[Y] = 50, ma a (X) < <r(Y): il secondo investi- 


Speranza 

condizionale 


Schema di 
Bernoulli 


mento e’ piu’ "rischioso” del primo. 

Se, nella formula della speranza, sostituiamo i valori 
Pi = P[X = Xj] 

con le probabilità’ condizionali rispetto ad un evento A, otterre¬ 
mo una stima di X sotto l’ipotesi del verificarsi di A. Chiameremo 
SPERANZA CONDIZIONALE DI X DATO A la quantità’ 

E[X/A] = Ex ( P[(X = Xj)/A]. 

Se H,, H 2 ,.., H n e’ una partizione di fi in eventi misurabili e X 
e’ una variabile aleatoria, e’ 

E[X] = E[X/H,]P(H,) + E[X/H 2 ]P(H 2 ) +.. + E[X/H n ]P(H n ) 

formula analoga a quella riportata alla fine del paragrafo 1.2 e al¬ 
trettanto utile. 


1.5 SCHEMA DELLE PROVE RIPETUTE 


Consideriamo un "esperimento” aleatorio (es. lancio di un da¬ 
do regolare) nell’ambito del quale si possa realizzare un evento S 
che chiameremo "Successo” (es. "uscita del numero 6”). Sia ”p" 
la probabilità’ di successo e q = 1 - p la probabilità’ di #(S) = 
F ovvero la probabilità’ di "Fallimento”. 

Se pensiamo di ripetere l’esperimento piu’ volte, sempre nelle 
stesse condizioni, cosi’ che l’esito di una prova non venga influen¬ 
zato da quanto successo nelle precedenti, realizziamo quello che 
nel calcolo delle probabilità’ e’ detto SCHEMA DELLE PROVE 
RIPETUTE (s.p.r.) o SCHEMA DI BERNOULLI di parametro 

"p”. 

In uno s.p.r. si considera l’indefinita ripetizione dell’esperimen¬ 
to, nel senso che non viene posto, a priori, alcun limite al numero 
delle prove da eseguire. 

Esempi di s.p.r. sono le estrazioni dei numeri del lotto nelle di¬ 
verse settimane, o i successivi "giochi” alla roulette. 

Se indichiamo con Sj l’evento ("Successo" alla i-esima prova), 
la successione S,,S 2 ,..,S n ,.. e’ una famiglia (numerabile) comple¬ 
tamente indipendente. 

PRIMO PRINCIPIO DELLE PROVE RIPETUTE 

In uno s.p.r. di parametro p > 0, prima o poi si verificherà’ 
un successo. 


In effetti, l’evento 




Paradosso 
di Borei 


E = {sempre e solo fallimenti ) <= {fallimenti nelle prime N pro¬ 
ve) = E n , 
quindi 

P(E)^P(E N ) = P[FDFn..nF]= qq..q =q N . 

N volte N volte 

Se p>0 e’ q<l e q N può’ essere reso arbitrariamente piccolo 
(Ve>0 3N: q N <e) per cui 

P(E) = 0 e 

P(^(E)) = P [prima o poi si verificherà’ un successo] = 1. 

Il risultato non dipende dal particolare valore di ”p”, purché’ 
positivo, ed e’ possibile concepire delle situazioni che hanno del¬ 
l’incredibile: 

PARADOSSO DI BOREL: una scimmia che batta a caso i tasti 
di una telescrivente, scriverà’, prima o poi, il testo esatto delle Di¬ 
vina Commedia. (Si considera come prova la battitura di un nu¬ 
mero di caratteri pari a quelli della Divina Commedia -compresi 
spazi, interlinee ecc.- e come ”Successo” la redazione esatta del 
testo. La probabilità’ di successo in ogni singola prova e’ assai ri¬ 
dotta, ma non nulla). 

Certi di ottenere un successo, possiamo chiederci quanto dovre¬ 
mo attenderlo. Indichiamo con T il numero d’ordine della prima 
prova in cui si realizza l’evento S; T e’ una variabile aleatoria: 

T vale se Io schema inizia con 

1 S 

2 FS 

3 FFS 

4 FFFS 


Di questa variabile aleatoria avra’ senso calcolare il valore atte¬ 
so E[T]. 

Posto H = [Successo alla 1’ prova) 

K = [Fallimento alla 1’ prova) 

(H,K) e’ una partizione di ft, per cui (vedi paragrafo 1.4) 

E[T] = E[T/H] • P(H) + E[T/K] • P(K). 

Ma, 

E[T/H] = 1 : nell’ipotesi He’T=l 

E[T/K] = 1 +E[T] : nell’ipotesi K, la prima prova viene ”spre- 
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cata”, dalla seconda bisognerà’ attendere mediamente altre E[T] 
prove; 

P(H) = p ; P(K) = q 
quindi: 

E[T] = lp + (1 + E[T])q 
da cui, ricordando che p + q= 1, si ricava: 

E[T] = 1/p. 

SECONDO PRINCIPIO DELLE PROVE RIPETUTE 

In uno s.p.r. di parametro p>0, il numero medio delle prove 
da eseguire per avere il primo successo e’ 1/p. 

Ad esempio, per ottenere il primo 6 nel lancio del dado, occor¬ 
rerà’ lanciarlo, mediamente, 6 volte. 

Il paradosso di Borei diviene ora piu’ accettabile: la scimmia scri¬ 
verà’ la Divina Commedia, ma p e’ talmente piccola che il tempo 
medio di attesa e’ astronomico. 

Se, di uno s.p.r., si eseguono esattamente ”n” prove, ci si può’ 
chiedere quale sia la probabilità’ di avere un numero prefissato ”k” 
(O^k^n) di successi. Poiché’ il calcolo di tali probabilità’ richie¬ 
de alcune nozioni di calcolo combinatorio, a questa domanda si 
risponderà’ nel 4’ capitolo dove si affrontano appunto argomenti 
di calcolo combinatorio e la distribuzione binomiale (si veda an¬ 
che il problema "Compito in classe’’ nel 2’ capitolo). 


1.6 LEGGE DEI GRANDI NUMERI 


In questo paragrafo si affrontano tre argomenti: 

(a) Disuguaglianza di Cebiscev 

(b) Calcolo di media e varianza della frequenza di successi in uno 
s.p.r. 

(c) Legge dei grandi numeri (in relazione ad uno s.p.r.) 

Le dimostrazioni relative ai punti (a) e (b) sono piuttosto diffici¬ 
li e possono essere ignorate senza che ciò’ pregiudichi la compren¬ 
sione del discorso. 

Sia X una variabile aleatoria di legge (x^p;), di media n = E[X] 
e varianza a 2 (X). 


Disuguaglianza 
di Cebiscev 


Se X e’ un numero reale positivo arbitrario, vale la seguente di¬ 
suguaglianza (DISUGUAGLIANZA DI CEBISCEV): 

o*(X) 

P[|x - n |^x]< 

A 

In termini discorsivi, la probabilità’ che la variabile aleatoria as¬ 
suma valori che si discostano dalla media ”piu’ di X” non supera 
tF(X)/X 2 . 

Ciò’ significa che, se <F(X) e’ "piccola”, e’ poco probabile che 
X si discosti "molto” da ix. Se, ad esempio, poniamo n = 100 e 
o 2 (X)= 1, per X = 10 si ha: 

P[|X -100|^ 10] ^ 1/100 : 

la probabilità’ che X assuma valori inferiori a 90 o superiori a 110 
e’ minore di 1/100. 


DIMOSTRAZIONE 

Fissato X>0, indichiamo con x’j i valori x ( per i quali e’ 
jx’i - ix |^X ; si avra’ anche (x’j - fx) 2 ^\ 2 . 

L’insieme (|X - /x|^X) e’ l’unione disgiunta di {x’j} quindi 
P[|X - n\ ^X] = E p’j. 

Si ha dunque 

(F(X) = E (xj - a *) 2 Pi ^ Z: (Xj 1 - /*) 2 Pi - ^X 2 E Pi ’ = X 2 -P[|X-/tj^X]. 
Dividendo per X 2 si ottiene la disuguaglianza di Cebiscev. 


In uno s.p.r. di parametro ”p”, sia F n la variabile aleatoria "fre¬ 
quenza di successo nelle prime n prove" = Numero di successi/n. 
F n ha media p e varianza pq/n. 

DIMOSTRAZIONE 

Se indichiamo I k la variabile aleatoria che vale: 

1 se alla prova k si e’ avuto successo 
0 se alla prova k si e’ avuto fallimento 

1 

sara’ F n = — (I, +1 2 +.. + I n ). 
n 
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Legge dei 
grandi numeri 
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Avremo: 

E[I k ] = 1 ■ p +0-q = p e 

E[F„] = — (E[I,]+E[IJ + .. + E[I„]) = — = p 
n n 

ove si e’ tenuto conto che E e’ un’operatore lineare: 
E[aX + bY] = aE[X] + bE[Y] con X,Y var. aleatorie e a,b numeri 
reali. 

<F(F„) = E[(F n -p) 2 ] = E[F n 2 + p 2 -2F„ • p] = 

= E[F 2 „] + p 2 -2pE[FJ = E[F 2 „] + p 2 -2p 2 = E[F 2 n J-p 2 


Vediamo come calcolare E[F 2 „] : 


FÌ » = ^(I l + I 2 + - I n) 2 = 4d 2 , + I 2 2+- + I 2 „+ h § l I h I k ) 


Tì z 




ma I 2 k = I k (02 = 0 e 1 2 =1) quindi E[I 2 k ] = E[I k ] =p. 

Per h^k, I h I k vale 1 solo se si verifica un successo sia alla pro¬ 
va numero h sia alla prova numero k. Per l’indipendenza, ciò’ av¬ 
viene con probabilità’ p p = p 2 , quindi 

E[I h I k ] = 1-p 2 + 0(1 -p 2 ) = p 2 . 

Allora 

E[F 2 n ] = -J r [E[I 2 1 ]+E[I 2 2 ) + .. + E[I 2 n ] + ( E k E[I h I k ]] = 

1 

[p + p + .. + p + E p 2 ] 

n 2 * 

l’ultima somma contiene n(n-l) addendi uguali a p 2 per cui: 

E[F 2 n ]=— [np + n(n-1 )p 2 ] = — [np + n 2 p 2 -np 2 ] = — 1 + p 2 . 

n 2 n 2 n 


In definitiva: 


^(F n ) 


PO-P) 

n 


+ P 2 - P 2 


pq 

n 


Poiché’ <r(F n ) = pq/n , la varianza diminuisce al crescere del nu¬ 
mero delle prove; applicando la disuguaglianza di Cebiscev a F n , 
si ha (LEGGE DEI GRANDI NUMERI): 

VX>0 P[|F„-p|>AK x~ : 

A~n 



Osservazioni 


comunque si scelga un numero X>0, la probabilità’ che la frequenza 
di successo si discosti da ”p” per piu’ di X, può’ essere resa arbi¬ 
trariamente piccola, pur di considerare un numero di prove abba¬ 
stanza grande. 

Passando al complementare, e’ : 

pq 

vX>0 E [p-X<F n <p + X]> 1- - 7 -=-. 

X 2 -n 

Ciò’ significa che, comunque si scelga un intervallo (p-X,p + X), 
anche molto piccolo, la probabilità’ che la frequenza di successo 
cada in questo intervallo può’ essere resa arbitrariamente prossi¬ 
ma ad 1 scegliendo un numero di prove sufficientemente elevato. 

Attraverso una serie di osservazioni, rivediamo i risultati degli 
ultimi due paragrafi applicandoli al seguente esempio: 

Un dado regolare viene lanciato piu’ volte; chiamiamo "Succes- 
so” l’evento 

S = {uscita del numero 61 

I successivi lanci forniscono un esempio di s.p.r. di parametro 
p = P(S)= 1/6. 

1. Poiché’ p>0, se ripetiamo l’esperimento, prima o poi otterre¬ 
mo un 6 (primo principio delle prove ripetute) 

2. Quante volte occorre tirare il dado per avere il primo successo? 
Si potrebbe avere un 6 al primo lancio, oppure dover attendere il 
secondo o il terzo o il quarto ecc. Il numero di lanci "dipende dal 
caso’’ ma, per il secondo principio, occorreranno mediamente 1 /p 
= 6 tiri. 

3. Pensiamo di lanciare il dado un numero prefissato ”n” di vol¬ 
te; contiamo il numero ”k” di successi ottenuti; sara’ 0 ^k^n. 

Calcoliamo F n = k/n che chiamiamo "frequenza di successo” 
nelle n prove; avremo 0^F n ^ 1 (F n = 0 se non si e’ verificato al¬ 
cun successo; F n = 1 se in tutti i lanci e’ uscito 6 ). F n e’ una varia¬ 
bile aleatoria; calcoliamone il valore medio e la varianza: 

E[F n ] = p = 1/6 ; <r(F n ) = pq/n = 5/(36n) 

4. Cosa si può’ dire di F n quando n diviene molto grande? La leg¬ 
ge dei grandi numeri applicata a questo esempio diviene: 

VX > 0 E j^y-X < F n < y + xj > 1 - 5/(36X 2 n). 
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Guadagno 
medio teorico 


Prendiamo X = 1/100; e’ 1/6 =0.167; 1/6 - X 0.157; 1/6 + X = 
0.177 e 5/(36X 2 )= 1389 
quindi, approssimativamente, 

P[0.157 < F n < 0.177] > 1 - 1389/n. 

Supponiamo di eseguire 100000 prove (n= 100000), avremo: 

P[0.157 < F.oo .000 < 0.177] ^ 0.986. 

Eseguendo 100000 prove, la probabilità’ che la frequenza di suc¬ 
cesso sia compresa tra 0.157 e 0.177 e’ piu’ del 98% ! 

All’aumentare del numero delle prove, la frequenza ”si avvici¬ 
na” alla probabilità’ nel senso che "grandi differenze” tra frequen¬ 
za e probabilità’ divengono sempre meno probabili (attenzione, que¬ 
sto non significa che P(F n = p) si avvicini a 1; al contrario, la pro¬ 
babilità’ che la frequenza sia ESATTAMENTE uguale a ”p” ten¬ 
de a zero al crescere del numero delle prove). 

La stretta parentela tra "frequenza” e "probabilità” giustifica, 
almeno in parte, la concezione frequentista di cui si discuteva nel 
paragrafo 1 anche se, definire la probabilità’ come particolare fre¬ 
quenza, crea comunque delle difficolta’. Con le dovute cautele, si 
potrà’ utilizzare la frequenza come stima di una probabilità’ inco¬ 
gnita anche perche’, quando si ricorre a probabilità’ di questo ti¬ 
po, piu’ che il valore esatto interessa ”non sbagliarsi di molto”. 

Se, in uno s.p.r., all’evento "Successo” e’ associato un guada¬ 
gno ”x”, eseguendo ”n” prove ed ottenendo ”k” successi, si gua¬ 
dagnerà’ complessivamente G = kx. Il guadagno medio nella sin¬ 
gola prova e’ dunque GM = kx/n = F n x . 

Sostituendo F n con p, si ottiene px = E[X] dove X e’ la variabi¬ 
le aleatoria di legge 

X; : x 0 

Pi : P 1-P • 

Si spiega ora l’interpretazione di E[X] come "guadagno medio 
teorico” nei termini descritti nel paragrafo 1.4. 


1.7 SIMULAZIONI 


La maggior parte dei programmi associati agli esercizi contenuti 
nel prossimo capitolo realizzano delle simulazioni. 

Ciascuno dei problemi proposti si risolve calcolando certe quan¬ 
tità’; generalmente si tratta delle probabilità’ di particolari eventi 
e/o dei valori medi di certe variabili aleatorie. 
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Dopo aver risolto l’esercizio ”per via teorica”, perche’ non met¬ 
tere alla prova le nostre conclusioni? Probabilità’ e valori attesi pos¬ 
sono essere approssimati dalle frequenze in un opportuno s.p.r.; 
purché’ si eseguano molte prove in condizioni, per quanto possibi¬ 
le, immutate. 

La ripetizione di un esperimento aleatorio (es. lancio di un da¬ 
do) per centinaia o migliaia di volte e’ un lavoro lungo e noioso, 
se eseguito manualmente; un invito a nozze per l’instancabile cal¬ 
colatore. Quasi tutti i calcolatori dispongono di un generatore di 
numeri pseudocasuali (vedi paragrafo 1.8) che può’ essere utiliz¬ 
zato, con le opportune manipolazioni, per simulare l’esito di una 
prova di uno s.p.r. 

Un programma di simulazione sara’ dunque un programma nel 
quale viene ripetuto l’esperimento aleatorio (simulato) e nel quale 
vengono calcolate certe frequenze da confrontarsi con i valori attesi. 

Quando il problema riguarda un gioco in cui chi partecipa può’ 
scegliere tra differenti strategie, e’ possibile simulare contempora¬ 
neamente il comportamento dei diversi giocatori. Se abbiamo cal¬ 
colato che la strategia migliore e’ quella adottata dal sig. Rossi, 
ci fara’ piacere vedere che, di fatto, il nostro eroe si piazza in testa 
alla classifica, superando, magari di misura, avversari che adotta¬ 
vano strategie molto promettenti. 

Naturalmente non ci aspetteremo una perfetta sovrapposizione 
tra risultati teorici e valori simulati, ma quella ”giusta approssi¬ 
mazione” tra frequenza e probabilità’ nei termini previsti dalla legge 
dei grandi numeri (cosi’, se ci aspettiamo che A vinca piu’ spesso 
di B e succede il contrario, non e’ il caso di buttare il dischetto dal¬ 
la finestra). 

Vediamo, con un semplicissimo esempio, come si sviluppa il di¬ 
scorso problema - soluzione teorica - simulazione. 

Problema 

Flavio e Daniela, due bambini di sette anni, decidono di ripas¬ 
sare le tabelline giocando con due dadi, uno rosso e uno bianco. 
Flavio dice a Daniela: ''lanciamo i dadi e moltiplichiamo i due nu¬ 
meri; se il risultato e’ pari tu dai una caramella a me, se e’ dispari 
la do io a te”. 

Fa bene Daniela ad accettare la proposta di Flavio ? 

Soluzione 

Daniela vince se il prodotto dei valori e’ dispari e ciò’ avviene 
solo se su entrambi i dadi esce un numero dispari; se 

DR = {uscita dispari sul dado rosso j e 

DB = {uscita dispari sul dado bianco} e’ 

P(DR) = P(DB) = 1/2. 
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I due lanci sono indipendenti per cui 
P[Daniela vince] = E (DRfl DB) = P(DR)P(DB) = (l/2)(l/2) = 1/4. 
Quando Daniela perde, Flavio vince, quindi 
P[Flavio vince] = l-l/4 =3/4 

(oppure, dei 4 casi ugualmente possibili: PP, PD, DP, DD, i primi 
3 sono favorevoli a Flavio e solo l’ultimo a Daniela). 

Daniela farebbe bene a rifiutare la proposta di Flavio. 

Simulazione 

Cosa succede se i due bambini si mettono a giocare? 

Un programma di simulazione andra’ cosi’ strutturato: 

1- Richiesta del numero N di partite che si vogliono simulare 

2- Ciclo principale: per ogni partita, cioè’ da 1 a N, 

2.1 generazione casuale (o, meglio, pseudocasuale) del valore DR. 
risultato del lancio del dado rosso (DR sara’ un numero da 
1 a 6 costruito in modo che ogni valore abbia la stessa pro¬ 
babilità’ di essere scelto) 

2.2 operazione analoga per il dado bianco, valore DB 

2.3 calcolo del prodotto PR = DR*DB 

2.4 se PR e’ pari, viene incrementato di una unita’ il contatore 
VF, contatore delle partite vinte da Flavio 

2.5 se PR e’ dispari, viene incrementato di una unita’ il contato¬ 
re VD (partite vinte da Daniela) 

2.6 se le partite simulate sono meno di N, si ritorna al punto 2.1 

3- A questo punto VF e VD contengono il totale delle partite vinte 
da Flavio e Daniela rispettivamente; le frequenze saranno VF/N 
e VD/N, valori da confrontare con quelli attesi: 0.75 e 0.25 

Per chi ha un minimo di familiarità’ con il Basic, non dovrebbe 
essere difficile realizzare in pratica quanto descritto; la redazione 
del programma viene lasciata come esercizio. 


1.8 NUMERI PSEUDQCASUALI 

L’esigenza di disporre di lunghe sequenze di numeri casuali ha 
spinto i ricercatori a studiare una serie di procedure per la loro ge¬ 
nerazione. La maggior parte dei generatori lavora, grosso modo, 
secondo il seguente algoritmo: 


1 - viene richiesto un numero x (seme) appartenente ad un ceno 


Test di 
casualità’ 


Congruenza 

lineare 


intervallo (a,b) 

2 - viene calcolato un numero y = f(x) dove ”f” e’ costruita in mo¬ 
do che y e (a,b) 

3- il valore y viene sostituito a x 

4 - si rientra al punto 2 per il calcolo di un nuovo numero y. 

I numeri cosi’ generati vengono detti ”pseudocasuali”, in quanto, 
non solo la sequenza si ripete identica dopo un certo numero di 
passi (PERIODO), ma ogni valore risulta perfettamente determi¬ 
nato dal precedente. 

Praticamente pero’, quello che interessa e’ che la successione ”si 
comporti” come un’autentica successione casuale; il generatore vie¬ 
ne ritenuto soddisfacente se la sequenza che produce supera un certo 
numero di test di casualità’ (test che ne misurano appunto il "gra¬ 
do di casualità”’). 

Supponiamo di voler simulare i lanci di un dado regolare e di 
dover scegliere una delle seguenti successioni: 

Sequenza 1 

1234561234 

5612345612 

3456123456 

Sequenza 2 

1415265353 
2346264333 
25241 16335 

Mentre la prima sequenza appare a prima vista insoddisfacente 
(perche’, ad esempio, x i + 1 ^XjVi), la seconda sembra del tutto ca¬ 
suale, nonostante sia semplicemente costituita dalle prime 30 cifre 
dello sviluppo decimale di t (dove sono state ignorate quelle non 
comprese tra 1 e 6). 

Una trattazione dei test statistici di casualità’ richiede conoscen¬ 
ze matematiche superiori ed esula dagli scopi di questo libro; il let¬ 
tore interessato potrà’ consultare [7] e [9] dove sono descritti i test 
piu’ noti (test di uniformità’, test seriale, test del gap, test del po¬ 
ker, test della distanza). 

Per quanto riguarda la generazione di numeri pseudocasuali uni¬ 
formemente distribuiti (ovvero scelti indifferentemente ”a caso” 
in un intervallo), il metodo piu’ usato e’ quello della CONGRUEN¬ 
ZA LINEARE dove la relazione ricorrente e’ data da: 

X n + I = (A-X n + C) mod M 
(X n+1 e’ uguale al resto della divisione tra A-X n + C e M) 
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La funzione 
RND 
in Basic 


in cui A,C,M sono numeri interi e X 0 e’ compreso tra 0 e M. Per 
ottenere una buona successione con periodo sufficientemente este¬ 
so, A, C ed M vanno scelti in modo opportuno (si veda [7], pag.60), 
in particolare M deve essere piuttosto elevato, dato che il periodo 
non può’ superare M. Ecco un esempio: 

1 C = 99991:A = 24298:M= 199017:INPUT X 

2 N = A*X + C:X = N-INT(N/M)*M:PRINT X:GO TO 2 
Dividendo per M i valori cosi’ generati, si otterranno numeri 

pseudocasuali deU’intervallo [0,1). 

In Basic, la funzione RND fornisce appunto un valore pseudo¬ 
casuale compreso tra 0 e 1; generalmente X 0 non viene richiesto 
dall’esterno, ma viene automaticamente assegnato utilizzando va¬ 
lori forniti dall’orologio interno del calcolatore. La sintassi dell’i¬ 
struzione prevede anche la possibilità’ di ottenere la stessa sequen¬ 
za in momenti diversi (basta iniziare con lo stesso seme); un’op¬ 
zione che consente di ripetere in modo identico una simulazione 
(e, in effetti, la RIPRODUCIBILITÀ’ e’ un requisito importante 
di un buon generatore di numeri pseudocasuali). 



CAPITOLO 2 


PROBLEMI DI CALCOLO 
DELLE PROBABILITÀ 


In questo capitolo sono raccolti 17 problemi; a ciascun proble¬ 
ma e’ collegato un programma (a volte piu’ di uno) che calcola certi 
valori o realizza delle simulazioni, come specificato alla fine del 
capitolo precedente. 

Ho cercato, per quanto possibile e a costo di qualche ripetizio¬ 
ne, di mantenere gli esercizi "indipendenti” l’uno dall’altro: ogni 
problema e’ un paragrafo a se’ e contiene i richiami agli argomenti 
teorici necessari alla comprensione della soluzione (e, dove serve, 
qualcosa di piu’). 

In questo modo il lettore può’ scegliere liberamente gli esercizi 
che ritiene piu’ interessanti; gioverà’ sapere che sono disposti in 
ordine crescente di difficolta’ (i numeri 13, 14 e 15 sono piuttosto 
difficili). 

A volte il problema consiste nel determinare la miglior strategia 
per un certo gioco. In questi casi, e’ bene soffermarsi sul testo e 
procedere alla lettura della soluzione solo dopo aver ben compre¬ 
so il meccanismo del gioco; provate sempre a formulare una rispo¬ 
sta che, intuitivamente, vi sembra plausibile e, perche’ no, a gio¬ 
care qualche partita. 

Gli ultimi due esercizi (”I1 linguaggio” e "L’opinione pubbli¬ 
ca”) non sono dei problemi veri e propri, ma, per cosi’ dire, delle 
curiosità’: un microesempio di come modelli probabilistici possa¬ 
no essere utilizzati (e di fatto lo sono) in campi davvero insospettati. 

Per quanto riguarda i programmi di simulazione, va segnalato 
che il tempo di elaborazione può’ crescere considerevolmente al- 
l’aumentare del numero delle prove: nei primi esperimenti, con¬ 
verrà’ non esagerare nelle richieste. 

Con pochissime eccezioni, i programmi registrati ricalcano quelli 
inseriti nel testo (la velocita’ e’ stata leggermente sacrificata a fa¬ 
vore della "leggibilità”’ dei listati); una routine sonora si incarica 
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di segnalare la conclusione dell’elaborazione. 


2.1 CARLO E ANDREA 


Un programma associato: 

”CARLO E ANDREA’’ 

Un dado regolare a sei facce e’ stato modificato: su due facce e' 
rappresentata la cifra 1, su altre due la cifra 2 e sulle rimanenti 
due la cifra 3. 

Carlo e Andrea si alternano nel lancio del dado: il primo che 
non ottiene un valore maggiore di quello tirato dall’avversario nel 
turno precedente, perde il gioco (e l’altro vince). Inizia Carlo. 

Qual e’ la probabilità’ di vittoria di Carlo? Quale quella di An¬ 
drea? Qual e’ il numero medio di lanci per ogni singola partita? 
(durata media del gioco)? 

ESEMPI 

(C = tiro di Carlo; A 
NL = numero lanci) 

1 . C = 1; A = 1 

2. C = 3 

3. C= 1; A = 3 

4. C= 1; A = 2; C = 3 

SOLUZIONE 

Il problema, proprio perche’ semplicissimo, suggerisce immedia¬ 
tamente queste considerazioni: / 

- Il gioco non puo’terminare in parita’(si gioca ”al rialzo”) 

- Carlo e’ decisamente avvantaggiato 

- Il numero complessivo di lanci necessari a ”chiudere” una parti¬ 
ta va da 1 (esempio 2) a 3 (esempio 4). 

Il limitato numero di casi consente un elenco di tutte le partite 
possibili: tre lanci del ”nostro” dado originano 27 partite; a fian¬ 
co di ciascuna e’ scritta l’iniziale del vincitore e il numero di lanci 
effettuati per concludere la gara. Ad esempio, ”132 A 2” significa 
che Carlo ottiene 1 al 1’ lancio, Andrea 3 al 2’ lancio; la partita 
si conclude quindi dopo 2 lanci con la vittoria di Andrea (A). 


= tiro di Andrea; V = vincitore C/A; 

V = C; NL = 2 

V = C; NL= 1 

V = A; NL = 2 

V = C; NL = 3 


111 C 2 211 C 2 

112 C 2 212 C 2 


311 C 1 

312 C 1 


113 C 2 

121 A 3 

122 A 3 

123 C 3 

131 A 2 

132 A 2 

133 A 2 


213 

C 

2 

313 

C 

1 

221 

C 

2 

321 

c 

1 

222 

c 

2 

322 

c 

1 

223 

c 

2 

323 

c 

1 

231 

A 

2 

331 

c 

1 

232 

A 

2 

332 

c 

1 

233 

A 

2 

333 

c 

1 


Le "situazioni favorevoli" a Carlo sono dunque 19, quelle fa¬ 
vorevoli ad Andrea 8. Le probabilità’ richieste sono: 
p = P[Carlo vince il gioco] = 19/27 (= 0.704) 
q = P[Andrea vince il gioco] = 8/27 (- 0.296) 

La "durata media del gioco" (intesa come numero di lanci) e’ 
data dalla media dei valori della terza colonna ovvero 48/27 = 16/9 
= 1.778. 


Fig. 4 Grafo per la so¬ 
luzione del problema 
2.1 
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Un approccio di questo tipo e’ attuabile solo se il numero di "casi 
possibili” si mantiene moderato; alla soluzione si giunge molto piu’ 
rapidamente utilizzando il grafo di fig.4 (vedi paragrafo 1.2), ot¬ 
tenendo: 

P [vince Carlo] = 1/3 + 3/9 + 1/27 = 19/27 

P [vince Andrea] = 2/9 + 2/27 = 8/27. 

Un discorso analogo può’ essere fatto per il numero medio di 
lanci (paragrafo 1.4): 

1 lancio con probabilità’ 1/3 

2 lanci con probabilità’ 5/9 

3 lanci con probabilità’ 1/9 

da cui: 

1 (1/3) + 2 (5/9) + 3 (1/9) = 16/9 . 

SIMULAZIONE 

I risultati precedenti suggeriscono che, su 100 partite, Carlo to¬ 
talizzerà’ circa 70 vittorie contro le 30 di Andrea; il dado verrà’ 
lanciato mediamente 1.78 volte per concludere ogni singola gara. 

Nel programma che segue, RND(l) fornisce un numero casuale 
tra 0 e 1 e la funzione introdotta alla linea 10 determina un intero 
pseudo-casuale tra 1 e 3. L’utilizzo di variabili (U, D, T) in luogo 
di costanti (1,2, 3) riduce il tempo di elaborazione nel ciclo princi¬ 
pale (30-90). 

II programma richiede il numero di partite da giocare e fornisce 
la frequenza di successo di Carlo (che ”tendera’” a 0.704), quella 
di Andrea e il numero medio di lanci per partita. 


Le variabili utilizzate sono : 

NP = numero di partite che si intendono simulare 
LI, L2, L3 = risultato del 1’, 2’ e 3’ lancio 
C = contatore vittorie di Carlo 
A = contatore vittorie di Andrea 
L = contatore numero totale dei lanci eseguiti. 

Esempi : 

QUANTE PARTITE? 500 
FREQ. DI SUCCESSO PER CARLO: 0.704 
FREQ. DI SUCCESSO PER ANDREA: 0.296 
NUMERO MEDIO DI LANCI PER PARTITA: 1.782 



QUANTE PARTITE? 2000 
FREQ. DI SUCCESSO PER CARLO: 0.7045 
FREQ. DI SUCCESSO PER ANDREA: 0.2955 
NUMERO MEDIO DI LANCI PER PARTITA: 1.772 

Listato 

5 U=1:D=2:T=3 

10 DEF FND(X)=INT(RND(U)*T)+U 
20 INPUT "QUANTE PARTITE ";NP 
30 FOR I=U TO NP 

40 L1=FND(U):IF L1=T THEN C=C+U:L=L+U:G0‘T0 90 
50 L2=FND(U) : IF L2<=L1 THEN C=C+U:L=L+D:G0 T0 90 
60 IF L2=T THEN A=A+U:L=L+D:G0 T0 90 
70 L=L+T:L3=FND(U):IF L3=T THEN C=C+U:G0 T0 90 
80 A=A+U 
90 NEXT I 

100 PRINT''FREQ. CARLO'';C/NP 

110 PRINT''FREQ. ANDREA'';A/NP 

120 PRINT"NUM. MEDIO LANCI PER PARTITA" ;L/NP 


2.2 IL COLLEZIONISTA 


Un programma associato: 

”COLLEZIONISTA” 

Ogni scatola di detersivo contiene una figurina da incollare su 
un album in cui trovano posto N differenti figurine. 

Quante scatole di detersivo si dovranno comperare mediamente 
per completare la raccolta? Qual e’ la probabilità’ di cavarsela con 
esattamente N acquisti? 

Si tratta di un problema classico (si veda [15], pag. 72), che offre 
lo spunto per due programmi: un breve programma per il calcolo 
della soluzione e un programma di simulazione. 

Si ricorda che in uno schema di Bernoulli (s.p.r.) di parametro 
p>0, il ”tempo medio di attesa” per il primo successo e’ 1/p . 

SOLUZIONE 

La prima scatola va sempre bene, dato che la figurina troverà’ 
posto sull’album. Il secondo acquisto e’ la prima prova di uno sche¬ 
ma delle prove ripetute in cui “successo” e’ trovare una figurina 
diversa dalla prima. La probabilità’ di successo e’ p = (N-l)/N e 
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il numero medio di scatole da acquistare e’ allora N/(N-1). 
In modo analogo si procede per gli ulteriori acquisti: 


FIGURINA 

UTILE 

CASI 

FAVOREVOLI 

PROB. DI 

SUCCESSO 

NUM. MEDIO DI SCA¬ 
TOLE DA ACQUISTARE 

r 

tutti 

1 

1 

2 ’ 

N - 1 

(N - 1)/N 

N/(N - 1) 

3’ ‘ 

N - 2 

(N - 2)/N 

N/(N - 2) 

N-esima 

1 

1/N 

N 


Il numero medio di scatole da acquistare e’: 


M = 1 


N N 

-+-+ 

N-l N-2 



N(1 



1 

+ —). 
N 


Completare la raccolta con esattamente N acquisti significa col¬ 
lezionare N successi nelle prime N prove. La probabilità’ di tale 
evento e’ data da: 

N-l N-2 1 N! 

P = 1 • - --...-- —— 

N N N N N 


quantità’ questa che, anche per modesti valori di N, e’ assai picco¬ 
la (per esempio, per N=10 e’ P = 0.000359). 

Il programma 1 calcola il numero medio di acquisti necessari per 
il completamento di un album di N figurine. 

Listato programma 1 

10 INPUT "QUANTE FIGURINE'';N 

20 F0R 1=1 T0 N:S=S+1/I:NEXT I 

30 PRINT''NUMERO MEDIO DI ACQUISTI'';S#N 

Ecco alcuni valori ottenuti con il programma 1: 

N M circa 
10 29 

(*) 52 236 

100 519 

200 1176 

1000 7485 

(*) Molti anni fa erano in vendita scatole di formaggini con una carta da gioco di 
un mazzo di 52 carte in ogni confezione 
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Simulazione 


Nel programma 2 si realizza la simulazione. All’istruzione 5 vie¬ 
ne richiesto da quante figurine e’ composto l’album (NF). Il vetto¬ 
re A(..), dimensionato alla linea 10, rappresenta l’album: se A(I) = 0 
la casella I e’ vuota, se A(I) = 1 la casella I e’ riempita dalla relati¬ 
va figurina. 

Sempre nella 10, viene richiesto il numero di simulazioni, ovve¬ 
ro quanti collezionisti partecipano alla raccolta. 

Per ogni collezionista (ciclo 30-90), la variabile A contiene, in 
ogni momento, il numero di acquisti effettuati; la variabile N il 
numero di figurine presenti sull’album. Inizialmente (istruzioni 35 
e 36) questi valori vengono azzerati e l’album "svuotato”. 

All’istruzione 40, F contiene il numero della Figurina estratta (al 
solito RND(l) fornisce un valore pseudo-casuale tra 0 e 1) e, alla 
50, viene incrementato il numero degli acquisti. Se la figurina e’ 
utilizzabile, viene incrementato N e la figurina viene inserita nel¬ 
l’album (istruzione 60). 

L’istruzione 70 si occupa di controllare se la raccolta e’ comple¬ 
ta: in caso negativo si ritorna alla linea 40 per un nuovo acquisto, 
in caso affermativo viene stampato il numero del collezionista e 
il numero di acquisti sostenuti. 

La variabile T dell’istruzione 90 contiene il totale di tutti gli ac¬ 
quisti, che viene aggiornato. 

All’uscita del ciclo 30-90, viene stampato il numero medio di ac¬ 
quisti effettuati da ogni collezionista. Tale valore sara’, per un nu¬ 
mero di simulazioni sufficientemente elevato, abbastanza prossi¬ 
mo a quello previsto dal programma 1. 


Listato programma 2 

5 INPUT''NUMERO FIGURINE" ;NF 
10 DIM A(NF):Z=0:U=1:INPUT S 
30 FOR K=U T0 S 

35 FOR L=U T0 NF:A(L)=Z:NEXT L 

36 A=Z:N=Z 

40 F=INT(RND(U)*NF)+U 
50 A=A+U 

60 IF A(F)=Z THEN N=N+U:A(F)=U 
70 IF N=NF THEN 90 
80 G0 T0 40 

90 PRINT K,A:T=T+A:NEXT K 
120 PRINT "VALORE MEDIO"; T/S 
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Esempio di output del programma 2 


NUMERO FIGURINE : 200 
NUMERO SIMULAZIONI : 10 


COLLEZIONISTA N. 

ACQUISTI 

1 

1059 

2 

1602 

3 

1016 

4 

947 

5 

1360 

6 

1178 

7 

1246 

8 

904 

9 

1085 

10 

1205 


VALORE MEDIO : 1166.2 

(Il valore teorico per un album di 200 figurine e’ 1176). 

Nel programma "COLLEZIONISTA”, registrato sul supporto 
magnetico, sono condensati i due programmi: inizialmente viene 
chiesto il numero delle figurine; il programma stampa il numero 
medio di acquisti e richiede il numero di simulazioni procedendo 
quindi come nel programma 2. 


2.3 IL GIOCO DEL MINI-LOTTO 


Un programma associato: 

"MINI-LOTTO” 

Immagino che quanto segue risulterà’ ovvio per la maggior par¬ 
te di coloro che leggono queste pagine, ma la pubblicazione sui gior¬ 
nali di quei numeri del lotto che accumulano molte settimane di 
"ritardo” mi fa pensare che qualche considerazione possa risulta¬ 
re utile. 

Le conclusioni a cui giungeremo tramite il modello semplificato 
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rimangono valide per il normale gioco del lotto (liberandoci dalle 
complicazioni di carattere combinatorio). 

Ecco la favoletta: 

Molti cittadini di Bingotown giocano ogni settimana al Mini-lotto. 
Dal lunedi’ al venerdi’ e’ possibile puntare su uno dei dieci numeri 
0 , 1 , 2 ,..., 9 ; al sabato viene eseguita un’estrazione e coloro che 
hanno giocato il numero sorteggiato vincono. La somma per una 
giocata e’ fissa e piuttosto elevata; normalmente una persona non 
gioca piu’ di un numero a settimana. 

Il sig. Rossi e il sig. Bianchi sono fra i piu’ incalliti giocatori: 
ogni settimana si recano al botteghino. 

Rossi tiene aggiornata una ”tabella dei ritardi” Ri: ad ogni nu¬ 
mero ”i” (i —0,1,..,9) associa il numero di settimane trascorse dal¬ 
l’ultima uscita (Es. Se le uscite sono state 3/0/4/5/1/0/9, sara’ 
R0=1, Rl=2, R2 = 7, R3 = 6, ..., R9 = 0). Ogni venerdì’ Rossi gio¬ 
ca il numero con maggior ritardo (in caso di parità’, sceglie il mi¬ 
nore; nell’esempio precedente giocherà’ il 2 piuttosto del 6 o del 
7 o 8). 

Il sig. Bianchi agisce diversamente: ogni volta, con l’aiuto di un 
mazzo di carte ben mescolato, sceglie un numero a caso e lo gioca. 

Interminabili discussioni nascono tra i due: Rossi sostiene che 
la propria strategia e’ decisamente migliore e accusa Bianchi di ”non 
usare il cervello”; Bianchi ribatte che il suo metodo e’ altrettanto 
valido e che il "ritardato” non e’ il numero, ma una persona di 
sua conoscenza... 

Il programma MINI-LOTTO simula il comportamento dei due: 
Istr. 5: Inizializzazione; richiesta numero simulazioni (NS) 
Istr. 30-90: Loop principale 

- 30-40: in base al valore di massimo ritardo R(K), Rossi decide 
il numero X da giocare (inizialmente e’ R(K) = 0 per ogni K). 

- 40: Y e’ il numero che Bianchi intende giocare; E e’ il numero 
estratto. 

- 50-60: TR e TB (Totale vincite Rossi/Bianchi) vengono aggior¬ 
nati in caso di viincita. 

- 70-90: Aggiornamento della "tabella dei ritardi” (ad uso di Rossi). 

Istr. 100-130: Stampa risultati (numero di giocate vinte da Ros¬ 
si/Bianchi e relative percentuali) 


Listato del programma MINI-LOTTO 

5 Z=0:U=1:N=9:D=10:DIM R(N):INPUT NS 
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30 FOR I=U TO NS:M=R(Z):X=Z:FOR K= 

U TO N:IF R(K) >M THEN M=R(K):X=K 
40 NEXT K:Y=INT(RND(U)*D):E=INT(RND(U)*D) 

50 IF E=X THEN TR=TR+U 
60 IF E=Y THEN TB=TB+U 

70 FOR L=Z TO N:IF L=E THEN R(L)=Z:G0 TO 90 
80 R(L)=R(L)+U 
90 NEXT LrNEXT I 
100 PRINT 1 ' V.ROSSI " ;TR 
110 PRINT'» % ROSSI '';100*TR/NS 
120 PRINT" V.BIANCHI 1 ' ;TB 
130 PRINT" % BIANCHI" ;100*TB/NS 


Esempi di output 


QUANTE SIMUL.7500 

V. ROSSI 46 
<7o ROSSI 9.2 

V. BIANCHI 55 
% BIANCHI 11 


QUANTE SIMUL.? 1000 

V. ROSSI 119 
% ROSSI 11.9 

V. BIANCHI 100 
% BIANCHI 10 


Morale della favola: i risultati forniti dalla simulazione confer¬ 
mano che non esiste differenza fra le due strategie; sia Rossi che 
Bianchi vincono, mediamente, una volta ogni dieci settimane. 

A chi vuol saperne di piu’, consiglio la lettura dell’articolo ”Una 
tassa sulle illusioni” di Vinicio Villani, apparso sul numero di ot¬ 
tobre 1985 della rivista ”Archimede”. 


2.4 TIRO AL PIATTELLO 


Un programma associato: 



” PIATTELLO” 


Simulazione 


Nelle gare di tiro al piattello, Rambo e’ un vero fuoriclasse: in ogni 
prova colpisce il bersaglio con probabilità’ p = 0.97. 

Prima di iniziare, Rambo deve dichiarare quante prove intende 
eseguire: se non commetterà’ neppure un errore, vincerà’ un nu¬ 
mero di dollari pari al numero delle prove dichiarate (nulla in caso 
contrario). 

Qual e’ la dichiarazione piu’ conveniente? 

[si suppongono indipendenti gli esiti delle varie prove] 

SOLUZIONE 

Se ”x” e’ il numero delle prove dichiarate e ”p” la probabilità’ 
di successo nella singola prova, Rambo guadagna: 

”x” dollari con probabilità’ p x 

nulla con probabilità’ l-p x . 

Fissato pe(0,l) - i casi p = 0 e p = 1 sono di ovvia soluzione -, il 
guadagno medio e’ G(x) = x • p x e la dichiarazione ideale e’ quel¬ 
l’intero x’ che rende massima G(x). Consideriamo la differenza 

G(x+ l)-G(x) = (x+l)p x+1 - xp x = 

= (x+l)p x+1 -xp x+1 + xp x + 1 - xp x = 

= p x + l - xp x (1-p) = p x (p - xq) dove q = 1-p. 

Il valore da dichiarare sara’ il piu’ piccolo intero per cui 
G(x + l)^G(x) (dichiarare un solo tiro in piu’ peggiora le cose), 
dunque il primo intero per il quale p-xq ^ 0: x’ = -INT(-p/q) (la fun¬ 
zione INT(X) fornisce il piu’ grande intero minore o uguale a X; 
INT(-12.2) = -13). 

In alternativa, considerando x reale, si ha G’(x) = p x (l 4-x-log p) 
per cui G(x) ammette massimo per x= 1/jlog p|. Si tratterà’ allora 
di scegliere l’intero migliore tra INT(x) e INT(x)+l. 

Il programma PIATTELLO determina, in funzione di P, la di¬ 
chiarazione ideale X ed il relativo guadagno G(X); si passa quindi 
alla simulazione. 

In ogni gara (il loro numero e’ NS e N e’ il relativo contatore) 
vengono eseguite X prove (M contatore): La variabile S, numero 
di successi, viene incrementata solo quando le X prove hanno tut¬ 
te fornito esito positivo. 

Il rapporto X*S/NS esprime il guadagno medio realizzato da 
Rambo nella singola gara. Per P = 0.97, la dichiarazione ideale e’ 
33 associata ad un guadagno medio teorico di 12.0776 dollari (con¬ 
fronta fig.5). 
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Fig. 5 Andamento del guadagno in funzione dei tiri dichiarati (problema 2.4) 


Listato del programma PIATTELLO 

15 INPUT''PROB. SUCCESSO " ;P:X=-INT(-P/(l-P)) 

60 PRINT " CONVIENE DICHIARARE" ;X:PRI NT "GUADAGNO 
MEDIO'';X*PtX 

70 INPUT NS:FOR N=1 TO NS:M=X 
100 IF M THEN IF RND(l)<P THEN M=M-1:G0 TO 100 
110 IF M=0 THEN S=S+1 
120 NEXT N 

140 PRINT "GUAD. MEDIO SIMULATO " ;X*S/NS 


Esempi di output del programma PIATTELLO 


PROB. SUCCESSO? 0.97 
CONVIENE DICHIARARE 33 
GUADAGNO MEDIO PREVISTO 12.07761 
NUM. SIMULAZIONI? 1000 
GUAD. MEDIO OTTENUTO 12.606 

PROB. SUCCESSO? 0.5 
CONVIENE DICHIARARE 1 
GUADAGNO MEDIOPREVISTO 0.5 
NUM. SIMULAZIONI? 1000 
GUAD. MEDIO OTTENUTO 0.5 
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2.5 COMPITO IN CLASSE 


Distribuzione 

binomiale 


Due programmi associati: 

"BINOMIALE” 

"TEST” 

Una prova scritta di matematica e’ costituita da 10 quesiti ad ognu¬ 
no dei quali e’ associato un gruppo di quattro possibili risposte di 
cui una esatta e tre errate. 

Al candidato e’ richiesto di segnalare qual e’ l’alternativa che 
ritiene corretta; l’insegnante valuterà’ il compito assegnando un 
punto per ogni risposta esatta. 

Augusto, del tutto impreparato sugli argomenti oggetto della pro¬ 
va, sceglie a caso le risposte (in ogni quesito, ciascuna alternativa 
ha la stessa probabilità’ di essere scelta indipendentemente dalla 
domanda proposta e dalle risposte fornite alle domande precedenti). 

Detto X il numero di risposte esatte del compito di Augusto, cal¬ 
colare: 

1. P[X = k] per k = 0,l,2,...,10 

2. P[X ^ 6] probabilità’ che, nonostante tutto, la preparazio¬ 
ne di Augusto venga considerata ”non insufficiente”. 

SOLUZIONE 

Sia Ai = {Augusto risponde correttamente alla i-esima doman¬ 
da) i= 1,2,...,10; e’ p = P(Ai)=l/4 per ogni i. 

Le risposte fornite da Augusto realizzano uno schema delle pro¬ 
ve ripetute di parametro p= 1/4; il problema richiede di calcolare 
la probabilità’ di ottenere esattamente k successi (k = 0,l,...,10) nelle 
n = 10 prove eseguite. Tali probabilità’ sono (vedi 4.5 DISTRIBU¬ 
ZIONE BINOMIALE): 

pk q»-k dove q = 1-p 

Il programma BINOMIALE richiede i valori ”p” e ”n” e for¬ 
nisce b(k,n:p) per k = 0,l,...,n. 

Per comodità’ di lettura, i valori della tabella vengono presenta¬ 
ti con 4 cifre decimali (si tratta della forma tipica utilizzata nelle 
tavole). 

Successivamente si possono richiedere, per uno o piu’ k, i valori: 
P[X = k] con un miglior grado di approssimazione 

P[X^k] pari a y E [X = i] 

i = k 


b(k,n;p) 
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P[X<k] pari a l-P[X^k] 

Sintesi del listato del programma B1NOMIALE: 

60 INPUT P,N:DIM CB(N),T(N):Q=1-P:M=N:X=1: 
J=0:R=N:S=1:K=0 

70 IF M THEN X=X*Q:M=M-1:G0 T0 70 
80 CB(K)=X:T(K)=U:U=U+X 

: IF K<N THEN X=X*R/S*P/Q:R=R-1 :S=S+1 :K= 

K+1:G0 TO 80 

90 F0R K=0 TO N:PRINT K,CB(K),T(K):NEXT K 

La relazione ricorrente utilizzata per il calcolo di b(k,n;p) e’ ri¬ 
portata nel paragrafo 4.6; al termine dell’elaborazione CB(K) con¬ 
tiene P[X = k] e T(K) il valore di P[X<k] per k = 0,l,..,N. 

Con l’aiuto del programma BINOMIALE, risolviamo il proble¬ 
ma del compito in classe. Scegliendo P = 0.25 e N = 10 otterremo: 

K B(K,N;P) 

0 0.0563 

1 0.1877 

2 0.2816 

3 0.2503 

4 0.1460 

5 0.0584 

6 0.0162 

7 0.0031 

8 0.0004 

9 0.0000 

10 0.0000 

proseguendo con 
K = ? 6 
si ha: 

P[X = 6] =0.0162220001 
P[X> = 6] =0.0197277069 
P[X < 6] =0.980272293 

Cattive notizie per Augusto: la probabilità’ di superare la prova 
(P[X> =6]) e’ meno del 2°7o ... non sempre la fortuna aiuta gli 
audaci. 

Naturalmente, scegliendo opportunamente i valori di N e P, si 
possono inventare e risolvere diversi problemi legati allo schema 
delle prove ripetute (esempio: Cosa succede se, ad ogni quesito, 
vengono proposte 5 alternative - 4 errate + 1 esatta - ?... si utiliz¬ 
zerà’ il programma BINOMIALE con N= 10 e P = 0.2). 



SIMULAZIONE 

Nel programma TEST vengono simulate le prove di un gruppo 
di studenti che si comportano come Augusto. 

Inizialmente viene richiesto il numero di studenti (ovvero il nu¬ 
mero NS di simulazioni), quindi viene costruito il vettore R%(I) 
(RISPOSTE) 1 = 0,1,...,9 che contiene la soluzione del test. Suc¬ 
cessivamente, per ogni studente, viene calcolato il numero E di ri¬ 
sposte esatte (il tentativo dello studente INT(RND(1)*4) + 1 viene 
confrontato con R%(I) per I = 0,1,...,9). 

Il vettore A(E) contiene, in ogni momento, il numero degli stu¬ 
denti che hanno fornito esattamente E risposte corrette 
(E = 0,1 ,...,10); il rapporto A(E)/NS esprime la relativa frequen¬ 
za. Il programma prosegue calcolando la "frequenza delle prove 
sufficienti" ovvero (A(6) + A(7) + A(8) + A(9) + A(10))/NS. 

Listato del programma TEST 

20 DIM R*(9),A(10):INPUT NS 

40 FOR 1=0 T0 9:R$(I)=INT(RND(1)*4)+1:NEXT I 

50 FOR K=1 T0 NS:E=0 

60 FOR 1=0 T0 9:IF INT(RND(1)«4)+1=R*(I) THEN E=E+1 
70 NEXT I 

80 A(E)=A(E)+1:NEXT K 

90 FOR K=6 T0 10::W=W+A(K):NEXT K 

95 FOR 1=0 T0 10:PRINT I,A(I)/NS:NEXT I:PRINT W/NS 

Si confrontino i valori ottenuti dalla simulazione con quelli for¬ 
niti dal programma BINOMIALE per P = 0.25 e N=10. 

Esempio di output del programma TEST 

SIMULAZIONI ? 70 
K FREQ. 

0 0.0571428571 

1 0.2 

2 0.242857143 

3 0.257142857 

4 0.171428571 

5 0.0428571429 

6 0.0142857143 

7 0.0142857143 

8 0 

9 0 

10 0 

FREQ.SUFF. 0.0285714286 
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2.6 L’ESPLORATORE 


Due programmi associati: 

” ESPLORATORE 1” 

”ESPLORATORE 2” 

Un esploratore risulta disperso. Questi si e’ recato a ovest con pro¬ 
babilità’ Pw, oppure a est con probabilità’ Pe=l-Pw. 

Ogni squadra di soccorso, operando nella ”zona ovest”, lo tro¬ 
va con probabilità’ Rw (ammesso che si sia recato da quella par¬ 
te); operando nella ”zona est” lo troverebbe con probabilità’ Re. 

Le squadre disponibili, che operano indipendentemente l’una dal¬ 
l’altra, sono in numero di ”n”. Quante inviarne a ovest e quante 
a est per massimizzare la probabilità’ di ritrovamento? 

SOLUZIONE 

Prima di procedere con i calcoli, e’ buona cosa affrontare un 
esempio in modo da chiarire i termini del quesito e da evidenziar¬ 
ne la ”non banalità’” (il problema risulta banale se, ad esempio, 
Pw = 0; in tal caso tutte le squadre disponibili verranno indirizzate 
a est; altri casi ovvi si hanno se Re o Rw sono nulli: viene allora 
esclusa la possibilità’ di ritrovamento, in una certa zona, indipen¬ 
dentemente dal numero di squadre operanti). 

Supponiamo che sia Pw = 0.6 e Pe = 0.4 (questi potrebbero esse¬ 
re valori da noi stimati in base alle conoscenze delle intenzioni, de¬ 
gli interessi, delle predilezioni dell’esploratore), per cui e’ ”piu’ fa¬ 
cile” (ma non di molto) che il soggetto si sia recato a ovest. 

Purtroppo a ovest la vegetazione e’ assai intricata: una squadra 
di soccorso ritroverebbe l’esploratore con probabilità’ (Rw) di so¬ 
lo 0.1 (se c’e’!). 

A est le cose vanno meglio: la visibilità’ e’ migliore e una squa¬ 
dra lo troverebbe con probabilità’ (Re) pari a 0.2. 

L’operare indipendente delle squadre significa che la probabili¬ 
tà’ di ritrovamento da parte di due differenti squadre e’ pari al pro¬ 
dotto delle probabilità’ di ritrovamento di ciascuna di esse. 

A questo punto dobbiamo decidere come distribuire le squadre, 
diciamo 13, fra le due zone (vedremo che, nella situazione qui de¬ 
scritta, converrà’ mandarne 8 a ovest e 5 a est con probabilità’ di 
successo del 61%). 

Torniamo al caso generale: sia x il numero delle squadre inviate 
a ovest (O^x^n) e n-x il numero delle squadre inviate a est. 

Indichiamo con P(Ax) la probabilità’ che l’impresa abbia suc¬ 
cesso inviando x squadre a ovest. Si calcolerà’ P(Ax) variando x 
da 0 a n e si sceglierà’ il valore di x che rende massima tale quantità’. 




Soluzione 1 


Nel nostro caso risulta piu’ agevole il calcolo della probabilità’ 
dell’evento Bx, complementare di Ax: 

P(Bx) = P[l’esploratore NON viene ritrovato inviando x squadre 
a ovest (e, quindi, n-x a est)]. 

Sia H = ( l’esploratore si e’ recato a ovest}, 

K = {l’esploratore si e’ recato a est}; 

(H,K) e’ una partizione di fì per cui (vedi 2.1): 

P(Bx) = P(Bx/H)P(H) + P(Bx/K)P(K). 

Per il calcolo degli addendi, poniamo Fw= 1-Rw; Fw rappresenta 
la probabilità’ di "fallimento nella ricerca” di una singola squa¬ 
dra operante a ovest. Analogamente sia Fe= 1-Re. 

Si ha P(Bx/H) = (Fw) x : supponendo che l’esploratore si sia re¬ 
cato a ovest, l’evento Bx si verifica se fallisce sia la squadra 1, sia 
la squadra 2, sia la 3,..., sia la numero ”x” (la funzione vale an¬ 
che per x = 0: P(B 0 /H) = (Fw)°= 1: se nessuna squadra e’ inviata 
a ovest e l’esploratore si e’ recato proprio a ovest, e’ certo che non 
sara’ ritrovato). Per l’indipendenza, la probabilità’ e’ (Fw)\ 
Analogamente P(Bx/K) e’ (Fe) nx in quanto tutte le n-x squadre 
devono fallire. 

Poiché’ P(H) = Pw e P(K) = Pe si ottiene: 

P(Bx) = Pw(Fw) x + Pe(Fe) nx . 

Il numero delle squadre da inviare a ovest può’ essere calcolato 
senza "scomodare l’analisi”; e’ infatti quell’intero x, con O^x^n, 
che rende minima P(Bx). 

Il programma ESPLORATORE 1 richiede: 

- il numero di squadre a disposizione (N:N> 1) 

- la probabilità’ che l’esploratore si sia recato a ovest (PW: 
O^PWs^l) 

- la probabilità’ che una singola squadra operante nella zona 
ovest riesca a localizzarlo (RW: O^RW^l) 

- come sopra, per la zona est (RE: O^RE^ 1) 

determina: 

la ripartizione delle n squadre che massimizza la probabilità’ di ri¬ 
trovamento e il valore di tale probabilità’ (su richiesta, vengono 
stampate le probabilità’ di ritrovamento per le altre ripartizioni nel 
formato "numero di squadre a ovest / probabilità’ di successo"). 

Listato del programma ESPLORATORE 1 
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10 INPUT'’N.SQUADRE’ 1 ;N 



30 INPUT''PROB.ESPL. A OVEST''jPW 
55 INPUT' 'PROB.RITROVAMENTO SQUADRA IN 
OVEST ";RW 

75 INPUT' 'PROB.RITROVAMENTO SQUADRA IN EST " ;RE 
100 PE=1-PW:FW=1-RW:FE=1-RE:MN=2 
110 FOR X=0 TO N 

120 PB=PW*FWtX+PE*FEt(N-X) :PRINT X,1-PB 
140 IF PB<MN THEN Y=X:MN=PB 
150 NEXT X 

160 PRINT' 'SQUADRE A OVEST' ' ;Y 
170 PRINT' 'SQUADRE A EST" ;N-Y 
180 PRINT''PROB. DI SUCCESSO'» jl-MN 

Nel ciclo principale (istr. 110-150), P(Bx) viene calcolato per tutte 
le possibili ripartizioni (X rappresenta il numero di squadre invia¬ 
te a ovest); quando, per un certo X, la probabilità’ di fallimento 
e’ minore di quella ottenuta nei casi precedenti, il valore di X vie¬ 
ne caricato in Y. A ciclo concluso, Y contiene la soluzione. 

Soluzione 2 II problema può’ anche essere risolto con ”un po’ di analisi”. 

Sia k = Pw; h = Pe; a = Fw; b = Fe; si tratta allora di determinare 
il minimo di 


f(x) = k-a x + h b nx per xe[0;n] 
con k + h = 1 . 

Se consideriamo a, b, h, k strettamente compresi tra 0 e 1, si ha: 
f’(x) = k -a x -log a - h-b nx -log b 
per cui f’(x)^0 se k-a x -log a ^h-b nx -log b. 

Poiché’ a < 1 e k > 0, e’ k-log a < 0 mentre b n x > 0; divi¬ 
dendo per b nx -k-log a (<0), si ha: 

a x h-log b 
b nx "" k-log a 

Posto M = (h-log b)/(k-log a), passando al logaritmo di en¬ 
trambi i membri, e’: 

x-log a - (n-x)-log b ^ log M 
ovvero 

x(log a + log b) ^ log M + n-log b. 

Poiché’ log a + log b < 0 (essendo a,be(0,l)), e’ 


log M + n-log b 

X $5 
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log a -l- log b 


= x’ 



In definitiva, per x = x’ si ha un minimo di f(x). 

Dato che f(x) e’ definita e continua per ogni x reale e x’ e’ l’uni¬ 
co minimo relativo, volendo limitare il calcolo del minimo in [0,n], 
si avra’ 

se 0 ^ x’ ^ n : minimo in x’ 
se x’ < 0 : minimo in 0 

se x’ > n : minimo in n. 

Visto che interessa un valore intero, prenderemo l’intero piu’ vi¬ 
cino a x’ (istr. 150). 

Nel programma ESPLORATORE 2, la soluzione viene trovata 
in questo modo senza che sia necessario passare in rassegna le va¬ 
rie ripartizioni come avveniva nel programma precedente. 

Listato del programma ESPLORATORE 2 

(le variabili N, PW ecc. corrispondono alle analoghe di ESPLO¬ 
RATORE 1) 

100 INPUT N,PW,RW,RE :PE=1-PW:FW=1-RW:FE=1-RE 
110 M=PE*LOG(FE)/(PW#LOG(FW)) 

120 Y1=(L0G(M)+N*L0G(FE))/(L0G(FW)+L0G(FE)) 

130 IF Y1<0 THEN Y1=0 
140 IF Y1>N THEN Y1=N 

150 Y1=INT(Y1+.5):PF=PW*FWtYl+PE#FEt(N-Yl) 

160 PRINT Y1,N-Y1,1-PF 

Esempi 


esempio 1 


Numero squadre 

10 

Prob. esploratore a W 

0.2 

Prob. ritrovamento a W 

0.9 

Prob. ritrovamento a E 

0.2 

Conviene inviare: 1 squadra a ovest 

9 squadre a est 

la probabilità’ di successo e’ 

circa 0.87 

esempio 2 


Numero squadre 

71 

Prob. esploratore a W 

0.45 

Prob. ritrovamento a W 

0.04 

Prob. ritrovamento a E 

0.047 

Conviene inviare: 34 squadre a ovest 


37 squadre a est 

la probabilità’ di successo e’ circa 0.795 
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2.7 IL DADO TRUCCATO 


Un programma associato: 

"BAYES” 

Un sacchetto contiene due dadi: 

- uno ”regolare”: uscite da 1 a 6 equiprobabili 

- uno "truccato”: lanciandolo si ottiene sempre 6. 

Un giocatore ha estratto uno dei dadi e, senza esaminarlo, l’ha 
tirato (una volta) ottenendo 6. 

Qual e’ la probabilità’ che abbia utilizzato il dado truccato? 

SOLUZIONE 

Posto H = {e’ stato lanciato il dado truccato) 

K = {e’ stato lanciato il dado regolare) 

A = {e’ uscito il 6), 
e’ HUK = fi e HflK = 0. 

Vogliamo calcolare P(H/A), cioè’ la probabilità’ che sia stato 
utilizzato il dado truccato, SAPENDO che e’ uscito il 6. 

Si tratta di una semplice applicazione del teorema di Bayes. Nel 
nostro esempio e’: 

P(A) = P(A/H)P(H) + P(A/K)P(K) 
e P(H/A) P(A) = P(A/H) P(H) = P(AflH) 
quindi 


Formula 
di Bayes 


P(H/A) = 


P (A/H) P (H) 

P (A/H) P(H) + P (A/K) • P(K) 


Poiché’ P(A/H)= 1; P(H) = P(K)= 1/2; P(A/K)= 1/6, si ottie¬ 
ne P(H/A) = 6/7 = 0.857142.. . 

Un confronto tra P(H) e P(H/A) mostra quanto l’ipotesi H venga 
avvalorata dall’informazione ”e’ uscito il 6” (si passa dal 50% a 
circa l’86%). 


SIMULAZIONE 

Nel breve programma che segue, le variabili T e B fungono da 
contatori: 

T conta quante volte viene estratto il dado truccato 

B conta il "numero di 6” ottenuti nelle N prove 
("A” rappresenta il valore uscito nella singola prova). 

Il rapporto T/B fornisce la frequenza di utilizzo del dado truc¬ 
cato nei lanci che hanno dato un 6. Nei termini della legge dei grandi 
numeri, T/B approssimerà’ P(H/A). 
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Listato del programma BAYES 

10 H=.5:U=1:S=6:INPUT"N. PROVE" ;N 
20 FOR I=U T0 N:IF RND(U)<H THEN T=T+U:A= 
S:G0 TO 40 

30 A=INT(RND(U)#S)+U 
40 IF A=S THEN B=B+U 
50 NEXT I 

60 PRINT' ’FREQ =";T/B 

Esempio di output del programma Bayes 

* BAYES * 

N. PROVE 1000 
FREQ. =0.847058824 

[Il valore teorico previsto e’ 6/7 = 0.857142] 


Metodo 
Monte Carlo 


2.8 CALCOLIAMO IL LOGARITMO 


Un programma associato: 

"LOGARITMO” 

Le tecniche matematiche che consentono la determinazione di 
un certo valore ”x” attraverso un un metodo aleatorio portano il 
nome di metodi Monte Carlo. In [3] viene proposto l’esempio clas¬ 
sico del calcolo di 7 r; con procedimento analogo viene qui mostra¬ 
to come valutare il logaritmo naturale di un numero maggiore di 
uno. 

In un sistema di assi cartesiani, consideriamo la funzione y = 1/x 
neH’intervallo [l,b] con b>l. 

Sia ABCD il rettangolo di vertici A = (1,0); B = (b,0); C = (b,l); 
D = (1,1) e ABED il trapezoide limitato dall’asse x, dal grafico di 
1/x e dalle rette x = 1 e x = b (quindi E = (b, 1/b) - vedi fig. 6 -). 

La probabilità’ che un punto P, scelto a caso nel rettangolo, ap¬ 
partenga anche al trapezoide e’ data dal rapporto delle aree delle 
due figure: 

p = S(trapezoide ABED)/S(rettangolo ABCD). 

Si ha: 


S(trapezoide ABED) = 


b dx 

! T 


= log b 
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Fig. 6 Calcolo di log b 
con il metodo Monte¬ 
carlo 



S(rettangolo ABCD) = (b - 1)1 = b - 1 
quindi p = (log b)/(b-l) ovvero log b = (b-l)p 

SIMULAZIONE 

Nel programma che segue, all’istruzione 10 viene richiesto il nu¬ 
mero B di cui si vuole calcolare il logaritmo e il numero di simula¬ 
zioni (NS). 

Il ciclo principale e’ costituito dalle istruzioni 30-50: ad ogni ite¬ 
razione viene scelto casualmente un punto di coordinate X,Y ap¬ 
partenente al rettangolo (RND(U) fornisce un valore pseudocasuale 
tra 0 e 1); alla linea 40 si controlla se tale punto appartiene anche 
al trapezoide, nel qual caso si incrementa il contatore S. 

Per NS sufficientemente elevato, S/NS e’ una buona approssi¬ 
mazione di p (nei termini della legge dei grandi numeri) e 
V = (B-1)*S/NS e’ il valore approssimato del logaritmo naturale 
di B (istr. 60). 

Il programma si conclude con la stampa del valore simulato (V), 
del valore fornito dalla macchina (LB), dell’errore assoluto e di 
quello relativo rispetto a LB. 

Listato del programma LOGARITMO 

10 INPUT B,NS:U=1:C=B-U:LB=L0G(B) 

30 F0R I=U T0 NS:Y=RND(U):X=C*RND(U)+U 
40 IF Y< =U/X THEN S=S+U * 

50 NEXT I 

60 V=C#S/NS:EA=ABS(V-LB):PRINT V,LB,EA,EA/LB 
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Esempio di output del programma LOGARITMO 


NUM. (> 1) 
NUM. SIM. 
SIMULATO 
LOG(B) 
ERR. ASS. 
ERR. REL. 


24 

1000 

3.243 

3.17805385 

0.0649 

0.0204 


(Il rapido appiattimento di 1/x per elevati valori della variabile e’ 
responsabile delle scarsa attendibilità’ dei risultati per "grandi” va¬ 
lori di b; qui si innesterebbe un discorso sulla precisione-macchina 
che esula da questa sede). 


2.9 L’EQUAZIONE ALEATORIA _ 

Un programma associato: 

"EQUAZIONE” 

I coefficienti a, b, c dell’equazione ax 2 + bx + c = 0 sono numeri 
reali scelti ”a caso” in un intervallo. 

Qual e’ la probabilità’ che l’equazione abbia soluzioni reali? 

SOLUZIONE 

Supponiamo a ^ 0 [e’ P(a = 0) = 0]. Senza perdere di generalità’, 
possiamo trasformare l’equazione in x 2 + dx + c = 0 dove ”d” e 
”c” sono due reali scelti a caso. 

L’equazione fornisce soluzioni reali se d 2 -4c^0 ovvero se: 

c^d 2 /4 (*). 

La probabilità’ che la coppia (d,c) di IR 2 soddisfi la (*) e’ pari 
alla probabilità’ che un punto P = (x,y), scelto a caso nel piano car¬ 
tesiano, sia situato ”al di sotto” del grafico di y = x 2 /4 (oppure sul 
grafico). 

Per ragioni di simmetria, basterà’ considerare il caso x^0. 
Sia r un numero reale positivo e t = r 2 /4. Scegliamo a caso un 
punto P appartenente al rettangolo R limitato dall’asse y, dalla retta 
x = r e dalle rette y = t e y = -t (vedi fig.7). 

La probabilità’ P(r) che P sia situato ”al di sotto” dell’arco OC 
e’ pari a 

Area(trapezoide OABC)/Area(rettangolo ABCD) = A(T)/A(R). 
E’ 
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Fig. 7 Probabilità 
" geometrica” 
(problema 2.9) 






1 

e A(R) = — r 
2 


3 


quindi P(r) = 2/3. 

Essendo ovviamente lim P(r) = 2/3, 2/3 e’ la probabilità’ ri¬ 
chiesta. r->oo 


SIMULAZIONE 

Nel programma EQUAZIONE vengono richiesti due valori: K 
e N. Per un fissato valore di K (> 1), il programma determina ”ca- 
sualmente” i tre coefficienti dell’equazione nell’intervallo (-K, + K) 
e, per ciascuna delle N simulazioni, controlla se l’equazione alea¬ 
toria fornisce o meno soluzioni reali. 

Per N sufficientemente elevato, ci si attende una frequenza pros¬ 
sima a 2/3; la frequenza ottenuta viene fornita al termine dell’ela¬ 
borazione in forma di rapporto (numero di ”successi”/numero pro¬ 
ve simulate) e in forma decimale. 

Dal punto di vista teorico, il valore K non e’ essenziale (si e’ vi¬ 
sto che P = 2/3 per ogni K > 0), ma può’ rivelarsi importante in re¬ 
lazione alla precisione-macchina. 
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Listato del programma EQUAZIONE 

30 INPUT " INTERVALLO (>1)";K 
35 DEF FNF(X)=2#K#RND(l)-K:INPUT''NUM.SIMUL. 1 ';N 
60 FOR 1=1 T0 N:A=FNF(l):B=FNF(l):C=FNF(l):D=B*B-4#A#C:IF 
D>= 0 THEN R=R+1 

70 NEXT I :PRINT''FREQ. = " ;R;''/•';N:PRINT''FREQ. = '';R/N 


Esempio di output del programma EQUAZIONE 

INTERVALLO (> 1)? 1000 

NUM. SIMULAZIONI? 300 
FREQ. = 185/300 
FREQ. = 0.616666667 
[P = 2/3 = 0.667 CIRCA] 


2.10 PING PONG 1 


Un programma associato: 

’ PING PONG 1” 

In un certo gioco, Aldo "guadagna un punto" contro Bruno con 
probabilità’ fissata ”p”. 

Vince la gara chi per primo raggiunge un vantaggio di due punti 
sull’altro. 

Calcolare: 

1. la probabilità’ che Aldo si aggiudichi la gara supponendo che 
le succcessive realizzazioni dei vari punti siano eventi mutualmen- 
te indipendenti; 

2. il numero medio di "giochi” necessari alla conclusione della par¬ 
tita (somma finale dei punteggi ottenuti da Aldo e Bruno). 

Questa "regola di torneo" viene utilizzata nella parte finale di 
giochi quali il tennis e il tennis da tavolo. 

Esempio: 


Azione 

Situazione 

Punti A. 

Punti B. 

— 

parità’ 

0 

0 

punto per A 

vantaggio A 

1 

0 

punto per B 

parità’ 

1 

1 

punto per B 

vantaggio B 

1 

2 

punto per B 

FINE, vince B 

1 

3 
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Fig. 8 Grafo per la soluzione del problema 2.10 

La gara si e’ conclusa con 4 "giochi”. 

SOLUZIONE 

Sia x la probabilità’ richiesta, ovvero 
x = P[Aldo vince la gara] = 

= P[Aldo vince con situazione iniziale di patita’] 
e sia q=l-p = P[Aldo perde un punto]. 

L’andamento del gioco può’ essere rappresentato dallo schema 
di fig.8. 

E’ 


x = p 2 + 2pqx 

(p 2 e’ la probabilità’ di raggiungere la situazione (**); 2pq e’ la 
probabilità’ di trovarsi nella situazione (*): da li’ la probabilità’ 
di aggiudicarsi la gara e’ nuovamente x). 

Ricavando x si ottiene: 

P 2 P : 

x = f(p) = -— - = ——“——~ 

1 - 2pq 2p- - 2p + 1 

Lo studio di f(p) in [0,1] (vedi fig.9), mostra come il meccani¬ 
smo di gioco amplifichi la probabilità’ di vittoria della gara quan- 
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Kig. 9 Probabilità di 
vincere la gara in fun¬ 
zione della probabilità 
di guadagnare il punto 
(problema 2. IO) 


Durata media 
della gara 



do la probabilità’ di guadagnare il punto sia maggiore di 0.5; que¬ 
sto significa che la "regola di torneo” contribuisce ad evidenziare 


il migliore: 


P 

f(p) 

0.5 

.50 

0.6 

0.69 

0.7 

0.84 

O 

OO 

0.94 

0.9 

0.99 


Sia D la "durata media della partita” ovvero il numero medio 
di prove neccessarie a chiudere la gara. Evidentemente e’ D^2, 
ricorrendo allo schema di fig.8, si ha: 

D = 2(p 2 + q 2 ) + (D + 2)2pq 
(E’ D = 2 se Aldo o Bruno realizzano subito 2 punti consecutivi, 
ciò’ avviene con probabilità’ p : + q 2 ; con probabilità’ 2pq si giunge 
ad una situazione di parila’: servono allora, mediamente, altre D 
prove per concludere). 

Ricavando D si ottiene: 


1 - 2pq 2p 2 - 2p + 1 
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Esempio 
di output 


con 1. 

Come prevedibile, g(0) = g(l) = 2 (se uno dei giocatori e’ imbat¬ 
tibile, la partita dura esattamente 2 giochi) e g(l/2 + t) = g(l/2 - 
t) per tf[0,1/2] (uno "scambio di abilita"’ tra Aldo e Bruno non 
altera la durata della gara). 

Derivando, si scopre che g’(p)^0 per p^ 1/2; la funzione am¬ 
mette un massimo nel punto M = (1/2,4): una partita tra giocatori 
di pari forza dura di piu’, mediamente 4 giochi. 

Il programma P1NG PONG 1 richiede P e fornisce f(P) e g(P); 
le istruzioni dalla 50 alla 90 realizzano la simulazione. 

Per ogni simulazione (1= 1,2,...,N), viene generato un valore 
pseudocasuale R dell’intervallo (0,1). 

Se R<P, Aldo guadagna un punto; in caso contrario lo guada¬ 
gna Bruno. 

La variabile S rappresenta il vantaggio relativo di Aldo su Bru¬ 
no; quando S*S = 4 la partita ha termine ed il totale delle vittorie 
di Aldo (T) viene incrementato se necessario; in questo caso si ri- 
pristina la situazione di parità’ e si procede ad una nuova simula¬ 
zione (istr. 90). 

Se S*S<4, la gara non c’ terminata: si rientra alla 60 per l'asse¬ 
gnazione di un nuovo punto; la variabile W conta il numero totale 
dei giochi effettuati. 

Alla 92 viene stampata la frequenza di vittoria di Aido (numero 
di gare vinte T sulle N simulate) da confrontare con il valore teori¬ 
co f(P); alla 94, operazione analoga per il numero medio di giochi 
per partita (W/N) da confrontare con g(P). 

Listato del programma PING PONG 1 

30 U=l:F=4:D=2:Z=0:INPUT f r PR0B. DI GUADAGNARE IL PUNTO 

, ,. P 

40 Q=1-P:FRINT f1 PR0B. DI VINCERE LA GARA 11 ;P*P/(l-2*P*Q) 
47 PRINT 1 'NUM. MEDIO TIRI PER GARA (PRE¬ 
VISTO) • 1 j2/(l-2*P*Q) 

50 INPUT 1 'NUM. SIMUL.* 1 ;N 
55 FOR I=U T0 N 

60 W=W+U:R=RND(U) :S=S-(R<P)+(R>=P) 

70 IF S*S=F THEN 90 
80 GO TO 60 

90 T=T-(S=D):S=Z:NEXT I 
92 PRINT 1 'FREQUENZA " jT/N 

94 PRINT 1 'NUM. MEDIO TIRI PER GARA (OTTENUTO)' 1 ;W/N 

Esempio di output del programma PING PONG 1 


PROB. DT GUADAGNARE IL PUNTO ? 


0.75 


PROB. DI VINCERE LA GARA 0.9 

NUM. MEDIO TIRI PER GARA (PREVISTO) 3.2 
NUM. SIMULAZIONI ? 517 

FREQUENZA 0.899419729 

NUM. MEDIO TIRI PER GARA (OTTENUTO)3.30754352 


2.11 PING PONG 2 


Un programma associato: 

”PING PONG 2” 

E' noto che nel gioco del tennis (anche nel tennis da tavolo, nella 
pallavolo ecc.), e’ preferibile trovarsi ”alìa battuta” (o ”in pos¬ 
sesso di palla”) piuttosto del contrario ed e’ realistico pensare che 
la probabilità’ dell’evento 

(il giocatore X guadagna un punto) 
sia maggiore quando il servizio tocca a X. 

Riprendiamo il problema proposto in PING PONG 1 : il regola¬ 
mento e’ lo stesso, ma Aldo e Bruno si alternano alla battuta. 

Aldo inizia: p e’ la sua probabilità’ di guadagnare il punto quando 
e’ in possesso di palla; p’ l’analoga probabilità’ nell’altro caso. 
Anche qui si vuol calcolare; 

1. la probabilità’ che .Xldo vinca la gara 

2. la durata media della gara. 

OSSERVAZIONE Se p = p’ si ritorna alla situazione descritta 
in PING PONG 1; converrà’ supporre p>p’. Per p<p’ si prospetta 
una situazione anomala, in cui Aldo gioca meglio quando non e’ 
in possesso di palla. 

SOLUZIONE 

Detto ”successo” la realizzazione di un punto da parte di Aldo, 
la situazione descritta determina uno schema delle prove ripetute 
di parametro p nelle prove dispari e parametro p’ nelle prove pari. 

Posto q = I-p e q’ = 1-p’, mantenendo le notazioni adottate in 
PING PONG 1, si ha lo schema di fig. 10. 

Se 

x = P[Aldo vince la gaia] 
e’ 

x = pp’ + (pq 1 + qp’)x 
da cui: 


PP 


PP’ 
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x 


1 - pq’ - qp’ 


2pp’ - (p + p’) + 1 



Fig. 10 Schema relati¬ 
vo aI problema 2.11 



Se D e’ la durata media della gara, e’ 

D = 2(pp’ + qq’) + (D + 2)(pq’ + qp’) 
da cui: 


1 - pq’ - qp’ 2pp’ - (p + p’) + 1 

(Se p = 1 ep’=0, o viceversa, la partita non ha fine: il denomina¬ 
tore si annulla e x e D perdono di significato). 

II programma PING PONG 2 richiede i valori p e p’ (variabili 
P e PP), fornisce x e D e realizza la simulazione in modo analogo 
a quanto avveniva in PING PONG 1. Nei vari giochi, la variabile 
L assume alternativamente i valori + 1 e -1: quando L = 1 la pro¬ 
babilità’ di successo e’ P, quando L = -l la probabilità’ di successo 
e’ PP. 

Listato del programma PING PONG 2 

10 U=1:F=4:D=2:Z=0:INPUT P,PP 
40 Q=1-P:QQ=1-PP:PRINT' ’x= f1 ;P*PP/(l-P*QQ-Q*PP) 

45 PRINT'’D=’’;2/(l-P*QQ-Q*PP) 

50 INPUT’’NUM. SIMUL. ”;N:F0R I=U T0 N:L=-U 
60 L=-L:W=W+U:R=RND(U) 

62 IF L=U THEN S=S-(R<P) + (R> =P) 

64 IF L=-U THEN S=S-(R<PP)+(R> =PP) 

70 IF S*S=F THEN 90 
80 G0 T0 60 

90 T=T-(S=D):S=Z:NEXT I 
95 PRINT 1 ’FREQUENZA " ;T/N 

97 PRINT "dJUM. MEDIO TIRI PER GARA (OTTENUTO) " ;W/N 
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Fig. 11 Probabilità di 
vincere la gara al va¬ 
riare delle probabilità 
eli guadagnare il pun¬ 
to della battuta p e 
alla risposta p’ (pro¬ 
blema 2.11) 


Fig. 12 Durata media 
della gara (problema 
2.H) 








Esempio di output de! programma PING PONG 2 


PROB. DI GUADAGNARE IL PUNTO 

- ALLA BATTUTA ? 0.7 

- ALLA RISPOSTA ? 0.4 

PROB. DI VINCERE LA GARA : 0.609 

NUM. MEDIO DI TIRI PER GARA (PREVISTO) :4.348 
* SIMULAZIONE* 

NUM. SIMULAZIONI ? 500 

FREQUENZA : 0.596 

NUM. MEDIO DI TIRI PER GARA (OTTENUTO) : 4.308 

Se p e’ prossimo a 1 e p’ prossimo a 0, Aldo, di norma, passa 
in vantaggio quando e’ al servizio e perde il vantaggio subito do¬ 
po. La durata media della gara diviene elevata e la simulazione ri¬ 
chiede tempo (per p = 0.99 c p' = 0.01 il numero medio di tiri per 
gara e’ 101.01). 


2.12 BOWLING _ 

Un programma associato 
"BOWLING” 

Ester sta concludendo una partita di bowling: due sono i birilli che 
le restano da abbattere. 

Da precedenti statistiche e’ noto che Ester, con un tiro, riesce 
ad abbattere: 

- 2 birilli con probabilità’ 0.1 

- 1 birillo con probabilità’ 0.4 

- nessun birillo con probabilità’ 0.5 

quando i birilli rimasti sono 2; 

- 1 birillo con probabilità’ 0.4 

- nessun birillo con probabilità’ 0.6 

quando resta da colpire l’ultimo birillo. 

Quanti tiri, mediamente, occorreranno ad Ester per terminare 
la partita? 

SOLUZIONE 

La situazione può’ essere rappresentata attraverso lo schema di 


Fig. 13 Grafo per la 
soluzione de! 
problema 2.12 



fig. 13 (il numero nella casella rappresenta il numero di birilli che 
restano). 

Sia X il numero medio di tiri che occorrono per completare la 
partita quando 2 sono i birilli da abbattere; Y l’analogo valore quan¬ 
do resta un solo birillo. 

Si ha: 

fY = 0.4*1 + 0.6*(1 + Y) 

(X = 0.1*1 + 0.4*(1 + Y) + 0.5*(1 + X) 

infatti: 

- se resta da abbattere un solo birillo, basterà’ un tiro con pro¬ 
babilità’ 0.4; con probabilità’ 0.6 il tiro viene ”sprecato” e, 
da li’, serviranno altri Y tiri; in totale (1 + Y) (equazione n. 

1 ) 

- se i birilli da abbattere sono 2, questi verranno eliminati con 
un tiro con probabilità’ 0.1; con probabilità’ 0.4 ne verrà’ 
abbattuto 1 e ne serviranno, mediamente, altri Y (quindi 

1 + Y); con probabilità’ 0.5 il tiro ”va a vuoto”: ne serviran¬ 
no altri X (equazione n. 2). 
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Risolvendo il sistema, si ottiene X = 4, Y = 10/4: ad Ester occor¬ 
reranno mediamente 4 tiri per terminare la partita. 

SIMULAZIONE 

Il programma BOWLING richiede il numero N di simulazioni 
e fornisce il numero medio di tiri che si sono resi necessari per com¬ 
pletare ciascuna delle N partite. 

Per ogni partita (1= 1,2,...,N), il numero di birilli B e’ inizial¬ 
mente posto uguale a 2, mentre viene azzerato il numero di tiri T 
(istruzione 60). 

Nell’istruzione 70 viene originato un numero casuale X dell’in¬ 
tervallo (0,1): 

- se B = 2 e X>0.5, B rimane 2 (nessun birillo abbattuto) 

- se B=2 e 0.1 <Xj$ 0.5, B diviene 1 (1 birillo abbattuto) 

(segue che, se B = 2 e X<0.1, B diviene 0: tutti abbattuti) 

- se B = 1 e X<0.6, B rimane 1 (non abbattuto) 

(segue che, con B=1 e X^O.6, B diviene 0: abbattuto). 

Finche’ Bf 0, la partita non e’ terminata, i tiri si susseguono 
(istruzione 70) e il loro numero viene opportunamente incremen¬ 
tato: T = T + 1. Quando B = 0, viene aggiornata la variabile S (to¬ 
tale dei tiri in tutte le simulazioni) e il programma procede ad una 
nuova simulazione ripristinando la situazione iniziale: B = 2 c T - 0. 
Al termine, S/N fornisce il valore cercato. 

Listato dei programma BOWLING 

50 INPUT N 

60 F0R 1=1 TO N:B=2:T=0 
70 X=RND(1) :B=(B=2)*(2*(X> .5) + (X> . 1 AND 
X< = .5))+(B=1)*(X< .6):T=T+1:IF BO0THEN 70 
90 S=S+'T :NEXT I:PRINT S/N 

Esempio di output del programma BOWLING 

QUANTE PARTITE VUOI SIMULARE ? 379 
NUMERO MEDIO OTTENUTO: 

3.92084433 TIRI 
PER PARTITA. 



2.13 DIVERSE USCITE 


Due programmi assodali: 

”1 TRE AMICI” 

’T SEI AMICI” 

Un giocatore lancia un dado regolare una o piu’ volte. Se, al ter¬ 
mine del gioco, avra’ ottenuto uscite tutte diverse, riceverà’ un nu¬ 
mero di dollari pari alla somma delle uscite stesse; nulla in caso 
contrario. 

Come deve comportarsi: 

(a) nel caso possa decidere di fermarsi quando lo desidera 

(b) nel caso debba dichiarare, prima di iniziare a giocare, il nume¬ 
ro di lanci che intende eseguire. 

SOLUZIONE (a) 

Si supponga di aver eseguito con successo ”n” lanci e siano 
u„u 2 ,...,u n i valori ottenuti (u^u. se i^j). 

Vogliamo considerare l’opportunità’ di procedere con un ulte¬ 
riore lancio (sara’ 1 <n^5; e’ ovvio che, se n = 6, ci si ferma gua¬ 
dagnando 21$). 

Sia A = 11,2,3,4,5,6) (valori possibili) 

U = {u,,u 2 ,..,uj (valori usciti) 

D = A\U (valori disponibili) 

Se ci si ferma, si guadagna la somma s = x 

Xrti 

se si prosegue, si guadagna: s + x se xeD 

0 se XeU. 

L’insieme D e’ formato da 6-n elementi e 

2 x - 2 x - 2 x = 21 - s - 

«t» x«A XfU 




Converrà’ continuare a giocare se G>s ovvero se 


5s - ns + 21 > 6s 
s + ns < 21 


s 


< 


21 

n + 1 




Calcolando l’espressione per i diversi valori di n, si ottiene il se¬ 
guente ”piano di gioco”. Nelle colonne ”min s” e ”Max s” sono 
riportati il minimo e il massimo valore di s per il corrispondente n. 

Un confronto tra questi valori e quelli ottenuti dalla (**) con¬ 
sente di raggiungere le conclusioni espresse nell’ultima colonna: 


land 

effettuati 

min 

s 

Max 

s 

conviene proseguire se 
“s” e’ minore di ... 

1 

1 

6 

10.5 ... sempre 

2 

3 

11 

7.0 /// 

3 

6 

15 

5.25 ... mai 

4 

10 

18 

4.2 ... mai 

5 

15 

20 

3.5 ... mai 


In definitiva, conviene lanciare comunque due volte e prosegui¬ 
re con un terzo, ultimo, lancio solo se, con i primi due, non si so¬ 
no raggiunti i 7$. 

Per il calcolo del guadagno medio realizzabile con questa strate¬ 
gia, si consideri la situazione dopo i primi due lanci; si avra’: 


s = 

con probabilità’ p s = 

# 3 

2/36 

# 4 

2/36 

# 5 

4/36 

# 6 

4/36 

7 

6/36 

8 

4/36 

9 

4/36 

10 

2/36 

11 

2/36 





Nei casi # si procede ad un terzo lancio e li’ il guadagno medio 
e’ dato dalla (*) con n = 2 (lanci eseguiti) ovvero (l/6)[3s + 21], Il 
valore richiesto e* dunque: 


6 


CM -J] 

s-3 


3s + 21 
6 



s-p s 


426 

216 


■ + 


152 

~36 


6.19444 


OSSERVAZIONE Operando 3 lanci, si hanno 216 casi possibili 
equiprobabili; il seguente programma calcola GM "senza andare 
troppo per il sottile” (la terna (X,Y,Z) viene fatta variare generan¬ 
do tutti i casi possibili) 


10 FOR X=1 T0 6 :FOR Y=1 T0 6:F0R Z=1 TO 6 
20 IF X=Y THEN 60 

30 A=X+Y:IF A> =7 THEN S=S+A:G0 TO 60 
40 IF Z=X OR Z=Y THEN 60 
50 S=S+A+Z 

60 NEXT Z : NEXT Y:NEXT X:PRINT M GM= S/216 


SOLUZIONE (b) 

Sia ”n” il numero di lanci dichiarato prima delPinizio del gioco 
(1 <n<6>; A= [1,2,3,4,5,61 
c 

B= {x = (x I ,x 2 ,..,x 11 ) con x s eA e x^Xj per i^j). 


L’insieme B e’ formato da D 6n elementi (vedi 4.1) e il guada¬ 
gno medio sara’: 




x k 


) 


Il numero complessivo degli addendi nella "doppia sommato¬ 
ria”, ovvero il numero totale degli x k> e’ n-D 6n e ogni elemento 
di A sara’ presente un ugual numero di volte N = n-D 6n /6. Ciò’ 
significa che: 


2 (È >.) 


1 ■ N + 2-N + 


+ 6-N = 21 -N 


da cui: 



21 n-D 


6,n 


2520 n 
6 n -(6 - n)! 
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[essendo D f>n = 6!/(6-n)! ] 
Calcolando al variare di n si ha: 


n 

Gn 

1 

3.50 

2 

5.83 

3 

5.83 

4 

3.88 

5 

1.62 

6 

0.32 


Conviene dunque dichiarare 2 oppure 3 lanci, guadagnando me¬ 
diamente 5.83$ (per scegliere tra 2 e 3, occorrerebbe calcolare la 
varianza). 

Come era prevedibile, il guadagno medio nella situazione (a) e’ 
maggiore di quello realizzabile nella situazione (b) (sempre adot¬ 
tando le strategie ottimali). 


LE SIMULAZIONI 
Caso (a): ”1 TRE AMICI” 

Il seguente programma di simulazione si propone due scopi: 

- valutare il guadagno medio effettivo realizzato da un giocatore 
(Bruno) che adotti la strategia ottimale (valore teorico: 6.19444.$) 

- confrontare i guadagni di Bruno con quelli ottenuti da Aldo e 
Carlo che adottano strategie non ottimali. 

Ecco come si comportano i tre partecipanti: 

ALDO. E’ un tipo molto prudente. Se al primo lancio ottiene 
6, si ferma (non volendo rischiare di perdere i 6 dollari); in caso 
contrario procede con un secondo, ultimo lancio. 

BRUNO. Fa uso della strategia ottimale descritta. 

CARLO. E’ amante del rischio. Lancia comunque fino a tre vol¬ 
te; se dopo i tre tiri ha accumulato solo 6 dollari (il minimo possi¬ 
bile), opta per un quarto, ultimo lancio. 

11 programma richiede il numero di simulazioni che si intendo¬ 
no effettuare c fornisce, per ogni giocatore: 




(1) il guadagno medio per partita - GM (X PAR) - 

(2) il relativo scostamento quadratico medio - DEV. STAN. - 

(3) il numero di partite risultate ”utili”, ovvero quante si sono con¬ 
cluse con un qualche guadagno - PAR. UTILI - 

(4) il "guadagno medio per partita utile”, ovvero il rapporto tra 
il totale guadagnato e il valore (3) -GM(XP.U.)- 

Principali variabili: 

Per 1 = 1,2,3 (rispettivamente Aldo, Bruno, Carlo); e’: 

G(I) : guadagno ottenuto nella partita in atto 

U(I) : contatore partite utili 

T(I) : somma dei dollari accumulati [ZG(I)] 

M(I) : somma dei quadrati dei guadagni [ZG(I)*G(I)] 

N$(I) : nome del giocatore 
inoltre: 

X,Y,Z,T : r,2’,3’,4’ uscita del dado 
N : numero delle simulazioni 

F : Flag: se F = 1 le prime 3 uscite del dado sono tutte 
differenti, in caso contrario e’ F = 0. 


Listato del programma 1 TRE AMICI 

5 DIM G(3),U(3),M(3),T(3),N$(3) 

7 N$(l) = "ALD0" :N$(2)=' 'BRUNO' ' :N$(3)=' 'CARLO" 

10 U=1:S=6:Q=4:A=0:D=2:H=3: INPUT "QUANTE PARTITE ";N 
30 FOR I=U T0 N 

40 X=INT(RND(1)*S)+U:Y=|INT(RND(U)*S)+U 
50 Z=INT(RND(l)*S)+U:T=INT(RND(U)*S)+U 
60 FOR K=U T0 H:G(K)=A:NEXT K:F=A 
70 IF X=S THEN G(U)=S 
80 IF X=Y THEN 140 
90 W=X+Y:G(U)=-W*(X<S)-S*(X=S) 

100 IF ZOX AND ZOY THEN F=U 

110 L=W+Z:G(D)=-W*(W>S)-L*(W<=S AND F=U) 

120 IF F=A THEN 140 

130 G(H)=-L*(L>S)-(L+T)*(T< >tX AND TOY AND TOlZ)* 

(l< =s) 

140 FOR K=U TO H:T(K)=T(K)+G(K) :U(K)=U(K)-(G(K) >0) : 

M(K)=M(K)+G(K)*G(K):NEXT K 
160 NEXT I 
220 FOR K=1 T0 3 

230 PRINT N$(K):ME=T(K)/N:DS=SQR(M(K)/N-ME*ME) 

235 GU=0:IF U(K)>0 THEN GU=T(K)/U(K) 

240 PRINT ME,DS,U(K),GU:NEXT K 
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Commenti al listato 


Istr. 5-10 : Introduzione (dimensionamenti ecc.) 

Istr. 30-160 : Loop principale : 

40-60 lancio dei dadi; azzeramento guadagni e flag 
70-90 calcolo del guadagno di Aldo 
100-110 calcolo del guadagno di Bruno 
120-130 calcolo del guadagno di Carlo 
140 aggiornamento T(..), U(..), M(..) 

Istr. 220-240: stampa risultati 

Esempio di output del programma I TRE AMICI 

NUM. PARTITE 1000 

** ALDO 

GM (X PAR) 5.498 
DEV. STAN. 2.89758451 
PAR.UTILI 852 

GM (X P.U.) 6.45305165 

** BRUNO 
GM (X PAR) 6.067 
DEV. STAN. 4.08907215 
PAR.UTILI 708 
GM (X P.U.) 8.56920904 

** CARLO 
GM (X PAR) 5.478 
DEV. STAN. 5.92735321 
PAR.UTILI 520 
GM (X P.U.) 10.5346154 

Un confronto tra la statistica di Aldo e quella di Carlo mostra 
come la strategia prudente del primo lo porti a guadagnare spesso 
(ben 852 partite utili su 1000) ma poco (6.45$ per ogni partita uti¬ 
le), al contrario di quanto avviene per Carlo che vince raramente 
(solo 520 partite utili), ma quando vince, vince bene (10.53$ per 
partita utile). 

I valori GM(XP.U.) si comprendono meglio se si nota che i mi¬ 
nimi e i massimi guadagni in una partita utile sono: 
per Aldo min = 3 $ (es. uscite 1-2) 

max =11$ (es. uscite 5-6) 
per Bruno min = 6 $ (es. uscite 1-2-3) 

max = 12 $ (es. uscite 4-2-6) 

per Carlo min = 7 $ (es. uscite 1-2-4) 

max = 15 $ (cs. uscite 4-5-6) 
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Naturalmente, il valore piu’ interessante e’ il primo — GM (X 
PAR) —,che, presumibilmente, risulterà’ maggiore per Bruno. 


Caso (b): ”1 SEI AMICI” 

In modo analogo a quanto visto per il caso (a), si vogliono con¬ 
frontare i risultati ottenuti da persone che si comportano diversa- 
mente. In questo caso simuleremo tutte le strategie possibili: pri¬ 
ma di iniziare il gioco, Aldo dichiara un lancio, Bruno ne dichiara 
2, Carlo 3, Dino 4, Ettore 5, Flavio 6. 

Dopo le assegnazioni e i dimensionamenti di rito, il programma 
richiede il numero di partite (isti*. 20). Le istruzioni dalla 30 alla 
200 costituiscono il ciclo principale; per ogni gara il dado viene lan¬ 
ciato un certo numero di volte: 
istr. 30 X(l) = uscita del V lancio 
istr. 60 T = uscita di un successivo lancio (2’,3’,..) 
istr. 80 si controlla se T e’ diverso dalle precedenti uscite... 
istr. 90 ...in caso negativo, viene posto X(K) = T e si ritorna alla 
60 per un nuovo lancio. 

L’istruzione 150 viene eseguita non appena T fornisce un valore 
identico a uno dei lanci precedenti (naturalmente non si possono 
avere piu’ di 6 uscite diverse); la variabile W contiene allora il mas¬ 
simo indice dei giocatori che realizzano un guadagno positivo. Per 
i giocatori 1,2,.. ,W viene valutato il guadagno G(L) e ne viene ag¬ 
giornata la situazione (T(L) = totale accumulato, U(L) = numero 
partite utili e Q(L) = somma dei quadrati dei guadagni). 

All’istruzione 180 i guadagni vengono azzerati e si procede ad 
una nuova gara. 

Esempio: 

1’ lancio 4: X = (4,.„.) 

2’ lancio T = 2: e’ TVX(l) allora X = (4,2,....) 

3’ lancio T = 6: e’ T^X(K) K= 1,2 quindi X = (4,2,6.) 

4’ lancio T = 2: e’ T = X(2): salto alla 150 con W = 3. 

Vincono Aldo, Bruno, Carlo c guadagnano: 

Aldo 4 $ : G(l) 

Bruno 6 $ : G(2) = G(1)+ X(2) 

Carlo 12 S : G(3) = G(2) + X(3) 

(Nel ciclo 150-170 viene utilizzato Tclcmcnto G(0) = 0). 

Per gli altri giocatori la situazione rimane immutata. 

Le istruzioni 400-420 realizzano la stampa dei risultati; con la 
successiva pressione di un tasto qualsiasi, vengono visualizzati i va¬ 
lori relativi a Bruno, poi Carlo ecc. 
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Listato del programma I SEI AMICI 

10 DIM X(6),G(6),T(6),U(6),Q(6),N$(6) 

11 FOR F=1 TO 6 :READN$(F):NEXT F 

12 DATA ALDO,BRUNO,CARLO,DINO,ETTORE,FLAVIO 
20 U=1:E=6:D=2:Z=0: INPUT "QUANTE PARTITE'';N 
30 FOR S=U TO N:X(U)=INT(KND(U)#E)+U 

60 FOR K=D TO S:T=INT(RND(U)*E)+U:W=K-U 
80 FOR H=U TO W:IF T=X(H) THEN 150 
90 NEXT H:X(K)=T:NEXT K:W=K-U 
150 FOR L=U TO W:G(L)=G(L-U)+X(L):T(L)=T(L)+G(L) 
170 U(L)=U(L)+U:Q(L)=Q(L)+G(L)*G(L):NEXT L 
180 FOR M=U TO E:G(M)=Z:NEXT M 
200 NEXT S 

400 FOR F=1 TO 6:PRINT F,N$(F):ME=T(F)/N:DS= 

SQR(Q(F)/N-ME#ME) 

415 GU=0 : IF U(F)>0 THEN GU=T(F)/U(F) 

420 PRINT ME,DS/J(F),GU:NEXT F 

Esempio di output del programma I SEI AMICI 

QUANTE PARTITE 1000 


1 ** ALDO 

GM (X PAR) 

3.526 

DEV. STAN. 

1.68858639 

PAR. UTILI 

1000 

GM (X P.U.) 

3.526 

2 ** BRUNO 

GM (X PAR) 

5.802 

DEV. STAN. 

3.21291083 

PAR. UTILI 

831 

GM (X P.U.) 

6.98194946 

3 ** CARLO 

GM (X PAR) 

5.895 

DEV. STAN. 

5.43378092 

PAR. UTILI 

563 

GM (X P.U.) 

10.4706927 

4 ** DINO 

GM (X PAR) 

3.94 

DEV. STAN. 

6.44394289 

PAR. UTILI 

278 

GM (X P.U.) 

14.1726619 



5 ** ETTORE 

GM (X PAR) 

1.53 

DEV. STAN. 

4.98268001 

PAR. UTILI 

87 

GM (X P.U.) 

17.58620069 

6 ** FLAVIO 

GM (X PAR) 

0.336 

DEV. STAN. 

2.63497704 

PAR. UTILI 

16 

GM (X P.U.) 

21 


Il "guadagno medio per partita" mostra, come ci si aspettava, 
che Bruno e Carlo adottano le strategie migliori. 

E’ possibile confrontare il numero di "partite utili" ottenute con 
il relativo valore atteso: per il giocatore di indice "k” (k = 1 Aldo, 
k = 2 Bruno,...), la probabilità’ che una partita si concluda con un 
guadagno positivo e’ p k = D ft k /6 k , per cui le partite utili saranno 
circa N-p k (N = numero simulazioni). 

Per N = 1000 si ha: 


ìndice giocatore 

valore atteso per P.U. (circa) 

1 

1000 (ovvio) 

2 

833 

3 

556 

4 

278 

5 

93 

6 

1 

15 


valori questi che ben si accordano con quelli simulati. 


2.14 AL RIALZO _ 

Un programma associato: 

"SFIDA" 

Un giocatore lancia un dado regolare una o piu’ volte. Fino a 
che il punteggio ottenuto e’ maggiore o uguale a quello dei lanci 
precedenti, egli riceve tanti dollari quanti sono i punti riportati, 
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« 


Strategie 


che si sommano al totale precedentemente guadagnato. Se ottiene 
un punteggio inferiore, il gioco ha termine e il giocatore perde tut¬ 
to quanto aveva accumulato. 

Egli può’ decidere di fermarsi quando lo desidera. 

Come deve comportarsi per massimizzare la sua speranza di 
guadagno? 

Esempio 1: 

esce ”1”: decide per un nuovo lancio 
..esce ”4”: decide per un nuovo lancio 
..esce ”4”: decide di fermarsi. 

Il giocatore ha guadagnato 9 $ 

Esempio 2: 

esce ”2”: decide per un nuovo lancio 
..esce ”1”: Il gioco ha termine con guadagno nullo. 

Il programma associato a questo problema consente all’utilizzato- 
re di definire la propria strategia e di misurarsi contro il calcolatore. 

La strategia ottimale, adottata dal calcolatore, verrà’ descritta 
nel paragrafo "soluzione”, paragrafo che andrà’ letto dopo aver 
giocato, magari piu’ volte, con il programma SFIDA (si eviti di 
esaminarne il listato). 

Nel problema proposto, una strategia e’ una "funzione di deci¬ 
sione” che, in base alle informazioni di cui disponiamo al momen¬ 
to, ci dice cosa dobbiamo fare: fermarci o continuare? 

Un esempio di strategia potrebbe essere: 

1. Lancio il dado la prima volta 

2. se oggi c’ lunedi’, lancio il dado una seconda volta e mi fermo; 
se non c’ lunedi’, mi fermo subito accettando l’uscita del primo 
lancio. 

E’ facile rendersi conto che (se ci si propone di guadagnare il 
piu’ possibile) quella descritta e’ una strategia insensata. 

Da tutte le strategie possibili isoliamo quelle piu’ promettenti. 
Due sono le informazioni che devono contribuire a determinare il 
nostro comportamento (dopo il primo lancio): 

- il valore ottenuto nell’ultimo lancio 

- la somma che si guadagna se si opta per ’Tarresto". 

Confrontiamo due esempi: 

Es.l A: 2/2/2/2/2/ 

B: 1/1/1/1/6/ 

Nei casi A e B sono stati eseguiti 5 lanci; in entrambi, se ci si 
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ferma, si guadagnano 10 S. Nel caso B, proseguire e’ assai piu’ ri¬ 
schioso che nel caso A: il risultalo dell’ultimo lancio e’ importante. 

F.s.2 A: 1/1/1/1/3/ 

B: 3/3/3/3/3/ 

Sia in A che in B, i lanci sono 5; l’ultima uscita c’ 3. Fermandosi 
si guadagnano 7 $ o 15 $ rispettivamente. Proseguendo nel caso 
B, si mette in gioco una somma maggiore che nel caso A: il ”tota1e 
accumulato” c’ importante. 

(Non interessa invece il numero di lanci eseguiti; se A: 3/3/ B; 
1/2/3/, i giocatori A e B si trovano in situazione identica). 

Una strategia che si rispetti saia’ cosi’ strutturala: 

1. Eseguo il primo lancio 

2. Se x e’ l’uscita dell’ultimo lancio eseguito c S e’ la somma che 
riceverei scegliendo di fermarmi, deciderò’ di 

PROSEGUIRE se S < f(x) (e rientro al punto 2) 
FERMARMI se S ^ f(x). 

Si tratta allora di scegliere i valori ”di soglia” associati ad ogni 
uscita. 

Avremo f(l) = ”infinito”: se l’ultimo lancio ha dato 1, conviene 
senz’altro proseguire. 

Inoltre sara’ f(2)^f(3)^ ... ^f(6): al crescere di x aumenta il ri¬ 
schio di perdere tutto e conviene accontentarsi di una somma 
minore. 

Un esempio di strategia ”sensata” potrebbe essere: 


X 

f(x) 

1 

infinito 

2 

30 

3 

20 

4 

10 

5 

5 

6 

2 


Nel programma SFIDA, l’utilizzatore e’ invitato ad impostare 
la propria strategia: per x = 2,3,...,6 si inserirà’ il valore f(x) che 
si ritiene piu’ conveniente. Il significato e’ ormai chiaro: porre 
f(2) = 30 equivale a dire ”se all’ultimo lancio ho ottenuto 2 e non 
ho ancora raggiunto i 30 $, intendo proseguire nel gioco”. 

Se, per qualche x, inserite un valore f(x)<x, intendete ”fermar- 
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vi” non appena il dado fornisce un’uscita pari a x. In questo caso, 
il valore effettivo di f(x) non ha importanza (l’effetto di f(6) = 2 
e f(6) = 3 e’ identico: non appena esce il 6 e’ S ^ 6 ed in entrambi 
i casi si opta per l’arresto). 

1 valori inseriti vengono caricati nelle componenti da 2 a 6 del 
vettore G(I) (istr. 40) precedentemente dimensionato (istr. 10) (G(0) 
e G(l) non vengono utilizzate). In modo analogo vengono caricati 
in C(I) i valori di soglia della strategia adottata dal calcolatore (istr. 
80). 

Dopo aver impostato il numero di partite che si intendono si¬ 
mulare (NS nella 100), ha inizio la sfida vera e propria al termine 
della quale il programma fornisce (istr. 600): 

- i II numero di dollari totalizzato dai giocatore 

- 2 II relativo guadagno medio per gara 

- 3 II numero di dollari totalizzato dal calcolatore 

- 4 II relativo guadagno medio per gara 

Poiché’ il calcolatore adotta la strategia migliore, c’e’ da aspet¬ 
tarsi che vinca la sfida (se vi batte ”di misura”, avete compreso 
perfettamente il meccanismo e sapete ben intuire le probabilità’ degli 
eventi in gioco; se il calcolatore ”stravinee”, evitate i giochi d’az¬ 
zardo; se siete voi a vincere, siete bravi ... e fortunati). 

Le variabili utilizzate nella simulazione sono: 

- U: valore fornito dal dado nell’ultimo lancio 

- G,C: flags. SeG =0, il giocatore intende ”proseguire”; se G= 1, 
il giocatore intende ”fermarsi”. ”C” e’ l’analogo valore.per il 
calcolatore 

- SG,SC: SG c’ la somma accumulata dal giocatore nella singola 
gara ovvero ciò’ che riceve se decide di fermarsi (cosi’ SC per il 
calcolatore) 

- TG,TC: totale dollari nelle varie partite per giocatore e calcola¬ 
tore rispettivamente 

- X: risultato del ”nuovo” lancio del dado 

- W: contatore delle simulazioni. 

La simulazione viene realizzata nel ciclo 120-540: 

200: viene eseguita una sola volta per ogni gara. U rappresenta qui 
l’uscita del primo lancio; le somme accumulate vengono poste ugua¬ 
li ad U; C e G azzerati. 

240-280: si controlla se il giocatore intende proseguire, in caso ne¬ 
gativo viene aggiornato il valore TG. 

300-340: come per 240-280 in relazione al calcolatore. 

360: se Giocatore e Computer sono entrambi ”fermi”, si procede 
ad una nuova gara. 



380: determina la nuova uscita del dado (X) 

400-440: se il giocatore ha deciso di proseguire e X^U, SG viene 
opportunamente incrementato; se X< U, la somma accumulata vie¬ 
ne azzerata e il giocatore ”si ferma” (SG = 0:G = 1). 

460-500: come per 400-440 in relazione al calcolatore 
520: U viene posto uguale a X; la gara prosegue (istr. 240). 

540: chiude il ciclo aperto alla 120. 

Il listato del programma SFIDA e’ riportato dopo la soluzione. 

Esempio di uilizzo del programma SFIDA 

SCEGLI I VALORI DI SOGLIA 
PER LA TUA STRATEGIA 
* 2 * 

? 30 

* 3 * 

? 20 

* 4 * 

? 10 

* 5 * 

? 5 

* 6 * 

? 2 


NUM. SIMULAZ.? 

500 

$ GIOC. 

2746 

G.M. GIOC. 

5.492 

$ COMP. 

2796 

G.M. COMP. 

5.592 

SOLUZIONE 



Supponiamo di aver effettuato con successo il lancio di indice 
”n” e di dover decidere se proccederc con l’n+l-esimo. Sia: 
”s” la somma accumulata con i primi n tiri 
”u” il valore ottenuto nell’ultimo lancio (n-csimo). 
RINUNCIANDO si guadagna ”s” 

PROSEGUENDO si guadagna: fs-f x se x ^ u 

: 0 se x < u 

dove x = valore uscito alTn-t- 1-esimo lancio. 

Nell’ultimo caso le uscite utili sono: u,u + che sono in nu¬ 
mero di 7-u. Poiché’ ognuno dei valori ha probabilità’ 1/6, deci¬ 
dendo di continuare si guadagnerà’ complessivamente: 

s + u con probabilità’ 1/6 

s + u + 1 con probabilità’ 1/6 
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s + 6 


con probabilità’ 1/6 


Il guadagno medio, per un giocatore che decida di proseguire, 
e’ dunque: 


1 


1 


1 


1 


g = -t(s + u) + — (s + u + 1) + ...+— (s + 6) = — 

6 6 6 6 l 


» 

(7-u)s + ^ 


E’: 




6-7 u(u - 1) 42 - u 2 + u 

2 2 = 2 


per cui 

r -, _ 1 14s - 2 us + 42 - u 3 + u 

6 ’ 2 “ 


Converrà continuare se G>s, ovvero se 

14s - 2us + 42 - u 2 + u > 12s 

2s(l - u) > u 2 - u - 42 
2s(u - 1) < (u + 6)0 - u) 


Per u — 1 la discguaglianza e’ vera per ogni s ed e’ ovvio che con¬ 
viene continuare. Per u>l conviene proseguire se: 

(u + 6) (7 - u) 
s < --2- i 

2 (u - 1) 

Calcolando l’espressione per u = 2,3,...,6 si ottiene il seguente 
piano di gioco” ovvero i valori ideali per f(u): 


se rultiino lancio 
ha dato... 

—--- 

... conviene proseguire 
se ’Y’ e’ minore di: 

1 

proseguire comunque 

2 

20 

3 

9 

4 

5 

(*)5 

2.75 

(*)6 

1.2 

- ——- - -_i 




Nei casi (*), si ha s ^ 5: non appena si ottiene 5 o 6 conviene 
sempre fermarsi. 

Listato del programma SFIDA 


10 DIM G(6),C(6) 

40 FOR 1=2 T0 6:INPUT G(I):NEXT I 

80 FOR 1=2 T0 6:READ C(I):NEXT I.-DATA 20,9,5,2.75,1.2 
100 INPUT''QUANTE SIMUL.";NS 
120 FOR W=1 T0 NS 

200 U=INT(RND(l)*6)+1:G=0:0=0:SG=U:SC=U 
240 IF G=1 THEN 300 

260 IF U=1 OR SG<G(U) THEN G=0:GO TO 300 
280 G=1:TG=TG+SG 
300 IF C=1 THEN 360 

320 IF U=1 OR SC<C(U) THEN C=0:GO TO 360 

340 C=1:TC=TC+SC 

360 IF G#C> 0 THEN 540 

380 X=INT(RND(1)*6)+1 

400 IF G=1 THEN 460 

420 IF X>=U THEN SG=SG+X:GO TO 460 

440 SG=0:G=1 

460 IF C=1 THEN 520 

480 IF X>=U THEN SC=SC+X:GO TO 520 

500 SC=0:C=1 

520 U=X:GO TO 240 

540 NEXT W 

600 PRINT TG,TG/NS,TC,TC/NS 


2.15 DADOMAX 


Un programma associato: 

’DADOMAX” 

Un giocatore dispone di un dado regolare a sei facce; dopo aver 
dichiarato quanti lanci intende effettuare (n), procederà’ con gli 
n lanci e riceverà’ un numero di dollari pari al massimo dei pun¬ 
teggi ottenuti. 

Per giocare, deve prima pagare una tassa proporzionale al nu¬ 
mero di lanci dichiarato: T = a-n. 

Esempio: sia a = 0.1 $, decido di lanciare 4 volte e quindi pago 
0.4$. 

Lancio: ottengo 2/472/3; ricevo 4$: ho guadagnato 3.6$. 
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Determinare, fissato ”a” (”costo di un lancio”), il valore n’ piu’ 
conveniente. 

SOLUZIONE 

Per ogni n fissato, sia x, il punteggio ottenuto all’i-esimo lan¬ 
cio ex’ = max x f . 

Vogliamo determinare la legge della variabile aleatoria x’ n , ov¬ 
vero P[x’ n = k] per k = 1,2,..,6. 

Fissati ”n” e ”k”, la quantità’ P[x’ n = k] può’ essere calcolata 
come rapporto ”casi favorevoli/casi possibili”. 

I ”casi possibili” sono tanti quanti le disposizioni con ripetizio¬ 
ne di 6 oggetti presi a gruppi di ”n”: 

D’m = 6 n 
(vedi 4.2). 

I ’ casi lavorevoli” sono costituiti da quelle disposizioni con ri¬ 
petizione che contengono ”k” come massima uscita: queste sono 
le disposizioni degli elementi < k eccettuate quelle che contengono 
solo elementi <k-l. 1 casi favorevoli sono dunque in numero di 

D Vn - D’ fcli „ = k n - (k-l)“. 

Allora 


Pr = P[x’ n = k] 


k n - (k - l) n 
6 n 


La variabile aleatoria x’ n assume quindi il valore k con proba¬ 
bilità’ p k con 1 $Ck^6. Il suo valore medio e’: 


E[x’j 


o 

-S 

k= ! 


k u - (k - l) 11 
6 n 


= 6 - 


§ ar 


[Il risultato si ottiene sviluppando la sommatoria e procedendo ad 
opportuni raccoglimenti, oppure, piu’ rapidamente, notando che 
E[x’ n ] = P(x’ n > 0) + P(x’„ > 1) + ... + P(x’ n > 5) 
con P(x’„^k) = (k/6) n per cui P(x’ n >k)= l-(k/6) n .] 

E[x’ n ] rappresenta dunque il ricavo atteso per un giocatore che 
decida di lanciare ”n” volte (come ci si aspettava, il valore si avvi¬ 
cina a 6 al crescere di n). 

Il guadagno medio e’: 

0(n, = 6-,g (|)")-a.n 



e il problema consiste nel determinare n’^l che rende massimo 
G(n). 

G(l)= 3.5-a>0 per a<3.5: se la tassa e’ maggiore o uguale a 
3.5 $, non conviene giocare. 

Per ogni a fissato minore di 3.5, e’ G(n)^6-a- n e il gioco perde 
di convenienza per n>6/a (risultato prevedibile dato che il ricavo 
non supera comunque i 6 $). 

Un modo per determinare n’ può’ essere quello di calcolare G(k) 
per k intero da 1 a INT(6/a) e scegliere il valore k che rende massi¬ 
ma G(..). 

Si può’ fare di meglio: consideriamo la differenza 


G(n+ 1) - G(n) 


È 

k= 1 


(*)■( 


k 

1 - — 
6 


- a = D, 


Poiché’ (k/6)< 1 e’ (k/6) n + 1 < (k/'6) n e anche D„,,- a < D n - 
a (vedi fig. 14) quindi: 

per n < n’ e’ G(n+ 1) > G(n) 

e per ogni n > n’ e’ G(n+ 1) < G(n). 

Questo significa che e’ possibile determinare la soluzione calco¬ 
lando G(k) k = 1,2,... fermandosi non appena G(..) inizia a di¬ 
minuire. 

Il programma DADOMAX determina, in questo modo, il nu¬ 
mero ideale di lanci in relazione alla tassa A fissata. 


Fig. 14 Grafico cìi 
G(n + l)-G(n) (pro¬ 
blema 2.15) 
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Listato del programma DADOMAX 


Processi di 
Markov 


10 INPUT’'TASSA PER LANCIO' 1 ;A:B=IKT(6/A):U=1:C=5:Z=0:Q=6 
30 POR I=’J TO B:S=Z 
50 FOR K=U TO C:S=S+(K/Q)tI:NEXT K 
60 G=Q-S-A#I 

70 IF G > M THEN M=G:N=I:G0 TO 80 
75 GO TO 90 
80 NEXT I 
90 PRINT N,M 

Esempi 


r 

a 

n’ 

G(n’) 

0.00001 

54 

5.999 

0.0001 

41 

5.995 

0.001 

29 

5.966 

0.01 

16 

5.784 

0.05 

9 

5.328 

0.1 

6 

4.960 

0.5 

2 

3.472 

1 

1 

2.500 


2.16 IL LINGUAGGIO 


Un programma associato: 

”LINGU AGGIO” 

PREMESSA 

11 programma che segue costituisce un esempio di applicazione 
dei PROCESSI DI MARKOV al calcolatore. 

Per la comprensione del programma, non e’ necessaria una co¬ 
noscenza completa dei processi di Markov; basterà’ pensare ad un 
”sistema” che può’ trovarsi in ”n” stati distinti e che può’ passa¬ 
re da uno stato all’altro (o rimanere nello stato attuale) con certe 
probabilità’. 

Quando la probabilità’ di passaggio dallo stato ”i” (presente) 
allo stato ”j” (futuro) dipende esclusivamente dalla situazione pre¬ 
sente e non da tutta ”la storia precedente” (ovvero da ”come vi 
si e’ giunti”), si ha un processo di Markov. 

Le probabilità’ di passaggio da uno stato all’altro vengono det¬ 
te probabilità’ di transizione. Se queste probabilità’ non variano 




al passare del tempo, il processo si dice stazionario. 

I valori delle probabilità’ di transizione vengono agevolmente rap¬ 
presentati in una matrice di transizione in cui Pelemento p^ rap¬ 
presenta la probabilità’ che il sistema, attualmente allo stato ”i”, 
passi, all’istante successivo, allo stato ”j” (il tempo viene qui con¬ 
siderato come grandezza discreta). 

Ad esempio, consideriamo un sistema con 3 possibili stati: A, 
B e C. 

La matrice di transizione potrebbe essere: 


a 

A 

B 

C al tempo t + 1 

da... 

A 

1/2 

1/4 

1/4 

al tempo t 

B 

0 

0 

1 


C 

1/2 

0 

1/2 


questo significa che, se il sistema si trova nello stato A, passera’, 
all’istante successivo, allo stato B con probabilità’ 0.25; allo stato 
C con probabilità’ 0.25, mentre con probabilità’ 0.5 rimarra’ nel¬ 
lo stato A. 

La seconda riga ci dice che, ogniqualvolta il sistema si trova nel¬ 
lo stato B, passa obbligatoriamente allo stato C; la terza ci infor¬ 
ma che dallo stato C si passa, con la stessa probabilità’, allo stato 
A o si rimane nello stato C. 

I valori della matrice, in quanto probabilità’, sono compresi tra 
0 e 1 e la somma dei valori di ogni riga e’ uguale a 1. 

Se il processo e’ stazionario, una possibile successione di stati 
del sistema proposto nell’esempio potrebbe essere AABCACC- 
CAAAAC..., mentre la successione ABACB... risulta inammissi¬ 
bile in quanto dallo stato B non si può’ passare allo stato A ne’ 
da C a B. 

IL LINGUAGGIO 

La successione delle lettere alfabetiche che formano le parole e’ 
strutturata in modo diverso da lingua a lingua. 

Nelle parole italiane, le coppie di lettere come ”ma”, ”ca”, ”le” 
sono senz’altro piu’ frequenti della coppia ”ps” (presente, ad esem¬ 
pio, in termini come ”lapsus” o "psicanalisi”); nella lingua ingle¬ 
se sono invece abbastanza frequenti le coppie di lettere come ”wh” 
(in ”when”, ”who”, ”why” ecc.), inesistenti nella lingua italiana. 

Consideriamo ora un sistema che può’ assumere 27 stati distin¬ 
ti: le 26 lettere dell’alfabeto + lo "spazio” inteso come carattere 
di separazione tra una parola e l’altra. 

Una frase sara’ intesa come una serie di passaggi da uno stato 
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all’altro, ad esempio nella frase ”10 CORRO” il sistema passa dallo 
stato iniziale ”1” allo stato ”0” quindi allo stato "SPAZIO”, poi 
allo stato ”C” eec. 

Nella prima parte del programma in esame, il calcolatore accet¬ 
ta in entrata un testo qualsiasi fino a quando non viene immesso 
il simbolo che segnala la fine del testo. In questa fase, simboli 
diversi dalle lettere alfabetiche vengono stampati sullo schermo, 
ma considerati come ”spazio”. 

Man mano che l’immissione del testo procede, il calcolatore ag¬ 
giorna i valori di una matrice quadrala T di 27x27 elementi: agli 
indici di linea e colonna corrispondono le lettere alfabetiche piu’ 
lo ” spazio” e ogni elemento della matrice rappresenta il numero 
di volte che una determinata lettera (indice I) e’ stata immediata¬ 
mente seguita da un’altra determinata getterà (indice J). 

Ad esempio, supponendo l’alfabeto costituito dalle sole tre let¬ 
tere A, B e C, l’immissione della stringa ABBACCCCACC* origi¬ 
nerebbe la matrice: 


’T’ 

A 

B 

c 

A 

0 

1 

2 

B 

1 

1 

0 

C 

1 

0 

4 


Una volta costruita, la matrice T viene trasformata in una ma¬ 
trice di transizione dividendo ogni elemento per la somma dei va¬ 
lori della propria linea (se i valori di una certa linea sono tutti nul¬ 
li, questa operazione non viene eseguita per detta linea). 

Proseguendo nell’esempio precedente c indicando con S(i) la som¬ 
ma dei valori dcll’i-esima linea, si ha S(l) = 3; S(2) = 2; S(3) = .5, per 
cui la matrice viene trasformata nella seguente: 


!- 

1 

A 

B 

c 

A 

0 

1/3 

2/3 

B 

1/2 

1/2 

0 

c 

1/5 

0 

4/5 


A questo punto il calcolatore stampa un testo casuale in cui la 
probabilità’ di successione di una data coppia di lettere viene for¬ 
nita dalla matrice T. L’output del testo prosegue in modo indefi¬ 
nito; premendo un tasto qualunque si ottiene la sospensione del- 






l’elaborazione che viene ripresa premendo ulteriormente un tasto. 

Se il testo in entrata e’ sufficientemente esteso (diciamo di alme¬ 
no 60 parole), quello in uscita risulta essere una buona simulazio¬ 
ne del lignaggio proposto: se inseriamo una pagina in lingua italia¬ 
na, otteniamo un testo casuale che ”sa di italiano”, cosi’ per il fran¬ 
cese, il tedesco ecc. 

Poiché’ lo ”spazio” e’ trattato alla stregua di qualsiasi altra let¬ 
tera, e' possibile ottenere in uscita un testo con qualche parola di 
lunghezza eccessiva; inoltre, se in entrata sono presenti parole con 
lettere doppie, es. ”io corro”, in uscita e’ possibile avere parole 
come ”corrro” o ”corrrro”. 

Nonostante queste limitazioni, il programma fornisce risultati 
interessanti che, naturalmente, sono tanto migliori quanto mag¬ 
giormente risulta esteso il testo in entrata. 

Listalo del programma LINGUAGGIO 

7 DIM A(26,26) 

9 T$= ' 1 ‘SPRINT CHR$( 17)h-CHRS( 29) +CHR$(29) + ' 1 -' 1 i 
10 GET A$ : IF A$= " '’ THEN 10 
20 IF A$='' THEN 205 
30 PRINT A$; :J=ASC(A$)-64:1=ASC(TS)-64 
60 IF J<1 OR J>26 THEN J=0 
70 IF I<1 OR I>26 THEN 1=0 
60 A(I,J)=A(I,J)+l:T$=A$:GO T0 10 
205 F0R 1=0 TO 26 
210 S=0 

220 F0R J=0 TO 26 :S=S+A(I,J):NEXT J 
240 IF S=0 THEN 280 

250 F0R J=0 TO 26 : A(I,J)=A(I,J)/S:NEXT 3 
280 NEXT I 
300 L=0 

310 X=RND(1):S=0:R$='' " :SS= r ’ " 

320 F0R 1=0 TO 26:S=S+A(L,I) 

340 IF S>X THEN 360 
350 NEXT I 

360 IF 1=0 OR 1=27 THEN PRINT" ";:G0 T0 380 
370 PRINT CHRS(1+64); 

380 L=I 

385 IF I>26 THEN L=26 
410 G0 T0 310 

Commenti al listato 

Istr. 10-30: accetta il testo, se terminato va alla 205 
Istr. 60-70: caratteri diversi dalle lettere alfabetiche vengono con¬ 
siderati ” spazi” 

Istr. 80: aggiornamento valori della matrice A 
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Istr. 205-280: trasforma A in matrice di transizione 
Istr. 310-370: sceglie il carattere da stampare in funzione dei va¬ 
lori di A - Nel ciclo 320-350, la variabile S, inizialmente nulla, vie¬ 
ne incrementata successivamente dai valori contenuti nella linea di 
indice L (ultima lettera stampata). Quando S supera X, l’indice I 
raggiunto definisce la nuova lettera da stampare cosicché’ la pro¬ 
babilità’ di stampa di una certa lettera e’ quella definita dalla linea 
L della matrice A-. 

Poiché’ la matrice A viene sistematicamente aggiornata, la cor¬ 
rezione di caratteri erroneamente introdotti risulterebbe impossi¬ 
bile. Nel programma registrato su disco, tuttavia, e’ consentito can¬ 
cellare e riscrivere i caratteri nell’ambito della stessa riga; questi 
vengono infatti caricati in una variabile di transito e la matrice A 
modificata solo a linea completa. 


2.17 L’OPINIONE PUBBLICA 


Due programmi associati: 

"OPINIONE/V” 

”OPINIONE/S” 

Nella rubrica "Ricreazioni al calcolatore" apparsa sul numero 
di giugno 1985 di "Le Scienze”, A.K.Dewdncy proponeva la re¬ 
dazione di cinque semplici programmi legati ad altrettanti proble¬ 
mi. Il programma che segue costituisce un ampliamento di uno dei 
temi proposti (VOTERS e ANTIVOTERS) e, con un’opportuna 
interpretazione dei due "partiti", un medico potrebbe vedervi un 
modello semplificato di fenomeni a lui familiari. 

Consideriamo una matrice di L linee e C colonne che identifi¬ 
chiamo con un rettangolo composto da TC = LxC "celle" (vedi 
fig .15). 

Pensiamo ora di incollare il lato AB con il lato CD, arrotolando 
il rettangolo in modo da avere un cilindro e di unire le circonferen¬ 
ze 7 1 e 7 2 come mostrato in fig. 16. 

La superficie del solido che si ottiene, ciò’ che i matematici chia¬ 
mano TORO (una specie di ciambella), e’ "priva di bordi". Pro¬ 
prio questa superficie quadrettata sara’ il teatro della nostra storia. 

Inizialmente, il programma OPINIONE assegna casualmente a 
ciascuna delle celle un elettore che può’ essere un CONSERVA¬ 
TORE (simbolo "*") o un RIFORMISTA (simbolo ”."); si co¬ 
mincia dunque da una situazione in cui i sostenitori dei due partiti 
sono ben mescolati sia come consistenza numerica sia come dispo¬ 
sizione geografica. 



Fig. 16 II rettangolo viene trasformato in toro 


A questo punto si apre il ciclo principale. Una coppia di celle 
contigue viene scelta a caso (viene generata una coppia di coordi¬ 
nate, diciamo (2,5) e viene selezionata, sempre casualmente, una 
delle quattro celle adiacenti; nell’esempio, la (1,5), la (3,5), la (2,4) 
o la (2,6)). 

Se i due elettori sono della stessa idea, non accade nulla e il pro¬ 
gramma rientra al punto precedente per la scelta di una nuova cop¬ 
pia. Quando vengono trovate opinioni differenti, la variabile "tem¬ 
po” viene incrementata e ciascuno degli occupanti cerca di con¬ 
vincere l’altro a "cambiare bandiera". 

I due elettori sono circondati da 10 vicini (parenti, amici, colle- 
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Fig. 17 Coppia di 
elettori circondata da 
10 “vicini” 
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ghi ...) la cui opinione può’ influenzare l’opera di convincimento 
di uno o dell'altro (figure 17 e 18). 

Nel programma OPINIONE, l’evento 

111 conservatore (*) convince il riformista (.) a passare dalla sua 
parte j 

ha probabilità’ 

P[*] = (l-2p)x/'10 + p 

dove ”x” e’ il numero di conservatori che circonda la coppia (x = 7 
in fig.17; x= 10 in fig. 18) e ”p” c’ una probabilità’ richiesta all’i¬ 
nizio del programma. 

Il valore assegnato a ”p” regola l’influenza di ”x” nel computo 
di P[*] (fig. 19): 

- Se p = Q e’ P[*j = x/10. In questo caso l’esito della disputa e’ for¬ 
temente influenzato dall’opinione del vicinato. Nei caso di fig. 18, 
e’ addirittura P[*] = 1 e l’elettore riformista viene certamente con¬ 
vinto a cambiare opinione (la pressione dell’opinione pubblica 
c’ schiacciante). 

- Quando 0<p< 1/2, c’e’, anche nel caso di fig. 18, la possibilità’ 
che sia l’elettore ”in minoranza” a convincere l’avversario. La 
probabilità’ che ciò’ si verifichi aumenta al crescere di ”p”. 

- Se p = 1/2 e’ P[*] = 1/2 qualunque sia x: l’opinione dei vicini non 
ha importanza. Con probabilità’ del 50% il conservatore convin¬ 
cerà’ il riformista, con probabilità’ del 50% succederà’ i! con¬ 
trario . 

In questo caso il programma OPINIONE si comporta come il pro¬ 
gramma VOTERS proposto sulla rivista ”Lc Scienze”. 

- Per p> 1/2 assistiamo ad una inversione di tendenza. 11 soggetto 
”piu’ forte” si lascia convincere facilmente a cambiar partito e 
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Fig. 19 Grafico di 
P[*/ = (l-2p)% 7 0+p 



la probabilità’ di questo evento cresce all’aumentare dei vicini che 
la pensavano come lui. 

Se assegniamo p= 1, si ha P[*] = l-x/10 e nel caso di fig. 18 sara’ 
l’elettore conservatore a trasformarsi in riformista. Per p = 1 il pro¬ 
gramma OPINIONE si comporta come il programma ANTIV'O- 
TERS di cui si parla nella rivista citata. 

Riassumendo 
INPUT (220-280) 

1 dimensioni del rettangolo (trasformato in toro) 

2 valore di ”p” (P) 

3 periodo di stampa W (ogni quante unita’ di tempo viene mostra- 
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numero 


ta la situazione) 

DETERMINAZIONE SCHEMA INIZIALE (310-360) 

CICLO PRINCIPALE 

4 scelta di una coppia di elettori contigui di diversa opinione; ag¬ 
giornamento variabile ”tempo” (380-440) 

5 calcolo di ”x” (450-580) e P[*] 

6 determinazione di un numero casuale R = RND(U) e [0,1]: se 
R > P[*] il riformista si trasforma in conservatore; se R < = P[*] 
avviene il contrario (590). 

7 se tutte le celle contengono Io stesso simbolo si conclude l’elabo¬ 
razione (380) 

8 stampa situazione, se necessario (380) 

9 rientro al punto 4. 

Per valori p< = 1/2, si va verso il "totalitarismo” e, prima o 
poi, il programma ha termine; per p> 1/2, il meccanismo fa si’ che 
gli elettori si distribuiscano equamente tra i due partiti (in media) 
e l’elaborazione prosegue senza concludersi (tendenzialmente). 

Naturalmente le simulazioni piu’ ”realistiche” si ottengono quan¬ 
do ”p” e’ 0 o prossimo allo zero. 

E’ interessante confrontare gli sviluppi nei casi p= 1/2 e p = 0 
(vedi fig.20). 



Fig. 20 L’uniformità viene raggiunta piu rapidamente per bassi valori di p 
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Quando p= 1/2, i vicini ”non hanno peso” e si assiste ad una 
lunga ed estenuante battaglia tra le due fazioni con capovolgimen¬ 
ti di fronte abbastanza frequenti. Lo schema si mantiene irregola¬ 
re per molto tempo e la convergenza e’ piuttosto lenta. 

Per p = 0, l’elettore "stabile” ha bisogno dell’approvazione de¬ 
gli altri: individui isolati cambiano presto bandiera e i sostenitori 
dei due partiti si radunano rapidamente nel tentativo di minimiz¬ 
zare la linea di contatto con gli avversari. 

Il programma OPINIONE e’ presente in due versioni a seconda 
che si desideri l’output su video (OPINIONE/V) o su stampante 
(OPINIONE/S). 

Soprattutto nel secondo caso, conviene scegliere con attenzione 
le dimensioni del rettangolo e il periodo di stampa: con rettangoli 
"grandi”, l’uniformità’ viene raggiunta dopo molto tempo ed e’ 
meglio evitare stampe troppo frequenti. 

Un consiglio: utilizzate OPINIONE/S solo dopo numerosi espe¬ 
rimenti con la versione video variando i parametri in ingresso (na¬ 
turalmente per p> 1/2 bisogna usare OPINIONE/V altrimenti si 
ottengono stampe a non finire). 

Il dischetto contiene, sganciata dai vari menù, la versione com¬ 
pilata del programma ‘opinione’: OPINtONV.EXE (MS DOS); 
c/opinione/v (C 64). 

ULTERIORI COMMENTI AL LISTATO 

Sub. 30-50 e 60-80: calcolo della posizione delle celle modulo L-l 
e C-l (trasformazione del rettangolo in toro) 

Sub. 100-210: stampa situazione 
220-280: richiesta dati 
310-360: schema iniziale 
380-440: scelta coppia elettori 
450-500: calcolo di ”x” 

590-600: risoluzione del confronto tra gli elettori scelti. 

Listato del programma OPINIONE 

10 G0 T0 220 

30 IF ML> LI THEN ML=ML-L: RETURN 
40 IF ML< Z THEN ML=ML+L: RETURN 
50 RETURN 

60 IF MOC1 THEN MC=MC-C:RETURN 
70 IF MC<Z THEN MC=MC+C:RETURN 
80 RETURN 

100 PRINT r 'UNITA' DI TEMPO" jT 
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120 FOR 1-Z TO LI:FOR J^Z TO CI 

140 IF A(I,J)=U THEN AS=A$+P$:GO TO 160 

150 A$-A$-fV$ 

160 NEXT J 

170 PRINT A$:A$='*' ! :NEXT I 

210 PRINT" # =’' ;NA; 1 ' .=' 1 ;TC-NA:RETURN 

220 INPUT''LINEE,COLONNE " ; L, C:L1=L-1:Cl= 

C-1:DIM A(L1,C1):TC=L*C 
280 INPUT P,V:K«(l-2*P)/10 
290 T=0:U=1:Z=0:D=2:MU=-1:F$=' '* ,! :VS=". “ 

310 FOR I=Z TO LI:FOR J=Z TO CI 

350 IF RND(l)<.5 THEN A(I,J)-U:NA=NA+U 

360 NEXT J,I 

370 IF NA=Z OR NA=TC THEN GOSUB 100: E® 

380 IF T-INT(T/W)*W ~Z THEN GOSUB 100 
390 T=T+U 

400 X=INT(L»RND(U)):Y=INT(C*RND(l)) 

410 IF RND(U)>.5 THEN 510 

420 MC=Y+U:GOSUB 60:XF A(X,Y)=A(X,KC) THEN 400 

440 XV=X:YV=MC:S=Z:ML=X-U:GOSUB 30 

450 FOR I=MU TO D:MC=Y+I:GOSUB 60:S=S+A(HL,MC):NEXT I 

460 ML=X+U:GOSUB 30 

470 FOR I=MU TO D:MC=Y+I:GOSUB 60:S=S+A(ML,HO):NEXT I 
480 MC=Y-U:GOSUB 60:S=S+A(X,MC) 

490 MC=Y+D:GOSUB 60:S=S+A(X,MC) 

500 GO TO 590 

510 ML-X+U:GOSUB 30:IF A(X,Y)=A(ML,Y) THEN 400 

530 XV=ML:YV=Y:S=Z:MC=Y-U:GOSUB 60 

540 FOR I=MU TO D:ML=X+I:GOSUB 30:S=S+A(ML,MC):NEXT T 

550 MC=Y+U:GOSUB 60 

560 FOR I=MU TO D:ML=X+I:GOSUB 30:S=S+A(ML,MC):NEXT I 
570 ML=X-U:GOSUB 30:S=S+A(ML,Y) 

580 ML=X+B:GOSUB 30:S=S+A(ML,Y) 

590 IF RND(U)<(K*S+P) THEN A(X,Y)=U:A(XV,YV)=U:NA= 
NA+U:GO TC 370 

600 A(X,Y)=Z:A(XV,YV)=Z:NA=NA-U:GO TO 370 


Esempio di output dei programma OPINIONE/S 

Parametri : L = 15, C = 20, P = 0, W = 100 
L’uniformità’ viene raggiunta in 831 unita’ di tempo. 
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UNITA’ DI TEMPO 600 



















UNITA’ DI TEMPO 800 



UNITA’ DI TEMPO 831 


*= 0 . 300 
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CAPITOLO 3 


ELEMENTI DI STATISTICA 


Nel capitolo precedente, problemi di calcolo delle probabilità’ 
di un certo interesse sono stati affrontati con l’ausilio di strumenti 
matematici relativamente semplici. Un’analoga operazione in re¬ 
lazione alla statistica risulta piu’ difficile in quanto problemi di sta¬ 
tistica si collocano a diversi e precisi livelli di difficolta’ c i piu’ 
interessanti (campionamenti, correlazioni ecc.) richiedono un ba¬ 
gaglio matematico avanzato. 

In queste pagine ci si limiterà’ pertanto agli argomenti di base; 
il lettore interessato ad approfondire questa parte troverà’ in ” Pro¬ 
babilità’ Statistica c Termodinamica” [3] il logico proseguimento 
e numerosi programmi. 


3.1 DATI 


Siete stati nominati responsabili della produzione del calzaturi¬ 
ficio X che inizierà’ l’attività’ nei prossimi mesi. Fra i moltissimi 
problemi da risolvere dovete naturalmente decidere, per ogni mo¬ 
dello di calzatura, quanti paia di scarpe delle varie misure intende¬ 
te produrre. 

Due sono le situazioni che dovete assolutamente evitare: che trop¬ 
pi potenziali clienti non trovino le scarpe adatte ai loro piedi e clic 
nei vostri magazzini si accumulino scarpe che nessuno vuole per¬ 
che’ troppo grandi o troppo piccole. 

Per organizzare un piano di produzione sensato (grossolani er¬ 
rori vi costerebbero il posto), vi servono delle informazioni e, se 
queste non sono già’ disponibili, dovete procurarvele, sintetizzarle 
ed interpretarle. 


97 


Universo 

statistico 


Carattere 


Rilevazioni 


Sintetizzare 
i dati 


Limitiamoci ad un tipo di calzatura destinata al pubblico ma¬ 
schile e supponiamo che siate interessati al solo mercato naziona¬ 
le: la popolazione maschile residente in Italia costituirà’ il vostro 
UNIVERSO STATISTICO (o COLLETTIVO STATISTICO); le 
singole persone saranno le UNITA’ STATISTICHE. 

Ciò’ che intendete esaminare in relazione alle varie unita’ stati¬ 
stiche, in questo caso il ”numero di scarpa”, e’ il CARATTERE 
che, in generale, può’ essere espresso da un numero ('risultato di 
un qualche tipo di misura) o dalla presenza o meno di certi attri¬ 
buti. Nel primo caso si parla di caratteri QUANTITATIVI (eta\ 
altezza, peso ecc.), nel secondo di caratteri QUALITATIVI (colo¬ 
re degli occhi, appartenenza ad un certo partito politico ecc.). 

Fissato l’universo statistico e il carattere, o i caratteri, oggetto 
di studio, e’ necessario procurarsi le informazioni ovvero RILE¬ 
VARE i DATI. La rilevazione può’ essere TOTALE, se riferita al¬ 
l’intero universo; o PARZIALE, se ri ferita ad un CAMPIONE ov¬ 
vero ad un sottoinsieme opportunamente scelto. 

Nella vostra situazione dovrete per forza accontentarvi di una 
rilevazione parziale. 11 campione dovrà’ essere RAPPRESENTA¬ 
TIVO dell’intero universo e di dimensioni tali da consentire l’esten¬ 
sione dei risultati a tutta ia popolazione entro ragionevoli margini 
di sicurezza. 

Poiché’ una rilevazione su vasta scala può’ essere assai costosa 
mentre una rilevazione troppo limitata può’ essere assai rischiosa, 
la determinazione di un campione adatto e’ un problema delicato; 
vi si accennerà’ nei quarto paragrafo di questo capitolo (Inferenza 
Statistica). 

Supponiamo risolto il problema del campionamento: ora dispo¬ 
nete di una consistente mole di dati; pagine fitte di numeri, risulta¬ 
ti della vostra indagine statistica. 

L’informazione che vi serve e’ contenuta nei dati in forma, per 
cosi’ dire, polverizzata: per ottenere qualcosa di utile e’ necessario 
un lavoro di sintesi. Si tratta dunque di elaborare i dati in modo 
da ottenere un insieme limitato di valori che riflettano caratteristi¬ 
che ritenute importanti. 

Abbandoniamo l’esempio del calzaturificio c consideriamo una 
situazione piu’ semplice: 

Nella TABELLA 3.1 sono riportati i voti di un compito di mate¬ 
matica svolto nella classe l’sez. A. Il professore ha raccolto que¬ 
sti dati per cercare una risposta ai seguenti quesiti: 

- il testo proposto era adeguato o eccessivamente difficile? 

- gli argomenti aggetto della prova sono stati recepiti allo stesso 
modo da lutti gli studenti? 
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Tabella 3.1 


6.8 

8.0 

7.1 

8.0 

5.8 

5.5 

4.5 

5.6 

5.1 ! 

6.6 

6.5 

5.5 

8.0 


3.8 

7.5 

5.9 

8.0 

5.0 

7.5 

7.0 

6.5 

6.8 

4.5 

6.1 

6.6 

3.9 


(Per i voti intermedi si e’ posto: 6+ = 6.1; 6 1/2 = 6.5; 6/7 = 
6.6 ; 7-- = 6.8; 7- = 6.9 ecc.) 


L’universo statitistico e’ la classe in esame e le unita’ statistiche 
sono i singoli studenti: la rilevazione e’ totale e il carattere esami¬ 
nato e’ il voto ottenuto, carattere quantitativo. 

Il numero di dati e’ talmente esiguo che, in realta’, si potrebbe 
tentare una risposta semplicemente analizzando la tabella, ma e’ 
facile immaginare situazioni in cui si considerano migliaia di dati 
la cui elencazione non dice granché’. Inoltre potrebbe risultare utile 
il confronto con un altro insieme di dati (ad esempio i voti ottenuti 
in una classe parallela in cui si e’ dedicato piu’ tempo agli esercizi) 
e la comparazione dei dati grezzi diverrebbe difficoltosa. 

Iniziamo il lavoro di sintesi che, grazie alle modeste dimensione 
della Tabella 3.1, potremo eseguire dettagliatamente senza ricor¬ 
rere al calcolatore. 

Una semplice riorganizzazione dei dati può’ già’ dire qualcosa; 
nella Tabella 3.2 i dati sono disposti, per righe, in ordine crescente: 


Tabella 3.2 


3.8 

3.9 

4.5 

4.5 

5.0 

5.0 

5.1 

5.5 

5.5 

5.6 

5.8 

5.9 

6.1 

6.5 

6.5 

6.6 

6.6 

6.8 

6.8 

7.0 

7.1 

7.5 

7.5 

8.0 

8.0 

8.0 

8.0 


Si nota ora che il voto piu’ basso e’ 3.8, quello piu’ alto e’ 8 
Se raggruppiamo i dati secondo la seguente 


Tabella 3.3 


i voti x : 

sono 

3^x<4 

2 

4<x<5 

2 

5^x<6 

8 

6^x<7 

7 

7^x<8 

4 

8^x<9 

4 
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Fig. 21 Rappresentazione grafica dei dati della Tabella 3.1 


e riportiamo i valori su un grafico (fig.21), saltano all’occhio altre 
caratteristiche: i voti estremi sono piu’ rari mentre la maggioranza 
si concentra tra il 5 e il 7. 

Vi sono molti modi per rappresentare un insieme di dati cd esi¬ 
stono in commercio programmi che forniscono rappresentazioni 
grafiche veramente suggestive (vedi fig. 22 e 23). 











Media 

aritmetica 


Consideriamo un insieme di dati relativi ad un carattere quanti¬ 
tativo; indichiamo con ”n” il loro numero c con Xj i= 1,2,...,n i 
valori espressi, in generale, da numeri reali (nell’esempio precedente 
e’ n = 27 = numero dei voti e x-, = voto). 

Al di la* delle possibili rappresentazioni grafiche, e’ molto utile 
calcolare degli INDICI numerici che evidenzino certe caratteristi¬ 
che dell’intero insieme; se un grafico fornisce una visione globale, 
un numero si presta ad essere confrontato con altri numeri e, co¬ 
me si vedrà’, la cosa può’ essere molto utile. 

Gli indici statistici piu’ frequentemente utilizzati si distinguono 
in indici di POSIZIONE e indici di DISPERSIONE; i primi forni¬ 
scono informazioni sulla "grandezza” dei dati, i secondi sulla lo¬ 
ro diversità’ reciproca. 

L’indice di posizione piu’ famoso e’ senz’altro la MEDIA ARIT¬ 
METICA: ”la media aritmetica di ”n” dati e’ quel valore che, so¬ 
stituito a ciascuno dei dati, ne lascia invariata la somma”. 

Detto M tale valore, e’, per definizione, 


x, + x 2 + 

da cui 


+ x n = M + M + ... + M = nM 


M - 


x, + x, + 


n 


+ X n = 


X, 
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ESERCIZIO 


Mostrare che, se x’ c x” sono rispettivamente il minimo e il mas¬ 
simo di Xj, e’ x’ ^ M ^ x”. 

Una proprietà’ caratteristica della media aritmetica e’: 

2 (* - M > = 0 

i 

ovvero ”la somma degli scarti dalla media aritmetica e’ zero”. 
In effetti: 

(x, - M) + (x 2 - M) +... + (x n - M) = 

= (Xj + x 2 + ... + x n ) - (M + M + ... + M) = nM - nM = 0. 
Per quanto riguarda la Tab. 3.1 si ha M = 167.1/27 = 6.189 

Media Quando i dati sono raggruppati e il dato x, compare nella stati- 

aritmetica stica f; volte, l’espressione di M e’ data da: 

ponderata 



f,x, + f 2 x 2 + ... + f k x k 

f] + f 2 + ... + f k 


dove ”k” e’ il numero di raggruppamenti (media aritmetica PON¬ 
DERATA). 

Volendo applicare la formula precedente ai dati della Tab. 3.3, 
potremo scegliere come Xj il valore centrale di ogni intervallo del¬ 
la suddivisione scelta (cioè’ considerare uguali a 8.5 tutti i valori 
compresi tra 8 e 9 ecc.), ottenendo: 

3.5x2 + 4.5x2 + 5.5x8 + 6.5x7 + 7.5x4+ 8.5x4 

M = - = 

2+2+8+7+4+4 

= 169.5/27 = 6.278 

La media ottenuta in questo modo e’ diversa dalla precedente 
poiché’, considerando il valore centrale, si commette un errore di 
approssimazione. 

Alla suddivisione in classi e alla media ponderata si ricorre so¬ 
prattutto quando il carattere quantitativo e’ CONTINUO ovvero 
quando il dato può’ assumere, almeno teoricamente, un qualsiasi 
valore reale di un certo intervallo. Cosi’, se viene misurata la sta¬ 
tura di un gruppo di persone, converrà’ contare il numero di indi- 
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Mediana 


vidui la cui altezza e’ compresa tra 175 e 180 em, tra 180 e 185 cm 
ecc. evitando, già’ in fase di raccolta dei dati, una eccessiva fran¬ 
tumazione. Diverso il discorso per un carattere DISCRETO come, 
ad esempio, il numero di persone che costituiscono un nucleo fa¬ 
miliare; una tabella del tipo: 


i nuclei familiari 
formati da 

sono 

1 persona 

10 

2 persone 

52 

3 persone 

60 

4 persone 

47 

ecc. 



contiene dati raggruppati per classi, ina, poiché’ si tratta di un ca¬ 
rattere discreto, la media aritmetica coincide con quella ponderata. 

Nel caso del voto, il carattere andrebbe considerato come discre¬ 
to, dato che esiste un numero finito di voti intermedi convenzio¬ 
nalmente assegnati (7 + , 7 1/2 ecc.), ma la suddivisione in classi 
di ampiezza ”1 voto” risulta comoda sia dal punto di vista grafico 
sia per gli scopi che si propone il professore. 


Un altro indice di posizione di uso molto frequente e’ la ME¬ 
DIANA (X) definita come ”valore centrale” dei dati quando que¬ 
sti siano ordinati in ordine non decrescente. 

Se x, < x 2 ^ x 3 ^ ... ^ x n e’ 

X = x (n+1V2 se n e’ dispari (valore centrale) 


mentre si assume 

X — ~~ [x n/2 + x^+J se n e’ pari (media aritmetica dei 2 
valori centrali) 


ESEMPIO 1 : 

DATI : 7; 2; 7; 10; 100 (n = 5) 
DATI ORDINATI : 2; 7; 7; 10; 100 

o 
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(n +1)/2 = 3; X = 7 




ESEMPIO 2 : 

DATI : 4; 1; 1; 6; 30; 5 (n = 6) 

DATI ORDINATI : 1; 1; 4; 5; 6; 30 

X = (4 + 5)/2 = 4.5 

Nell’esempio dei voti e’ n = 27 e la mediana sara’ il 14’ valore 
delle Tab. 3.2 : X = 6.5 il che’ vuol dire che una meta’ della classe 
ha ottenuto voti inferiori (o uguali) a 6.5* l’altra meta’ superiori 
(o uguali) a 6.5. 

Una grande differenza tra media e mediana sta ad indicare una 
distribuzione "sbilanciata”: un solo valore eccezionalmente gran¬ 
de o eccezionalmente piccolo influenza la media aritmetica ma non 
la mediana (Es. DATI 1; 2; 96 : M = 33; X = 2). 


Fig. 24 Grafico di 
= fJ lf + f2-tj + 

\!2t( 
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yuartili 


Moda 


Media aritmetica e mediana godono di notevoli proprietà’: la me¬ 
dia aritmetica rende minima la funzione f(t) = X (x r t) 2 ; la 

mediana minimizza la funzione g(t) = H |x r t|. Verifichiamole in 
un esempio: 


Sia n = 3 : x, = 1; x 2 = 2; x 3 = 12. La media aritmetica c’ 5, la me¬ 
diana e’ 2. 

La funzione 

f( t ) = (l-t) 2 + (2-t) 2 + (12-t) 2 = 3t 2 -30t + 149 
assume il valore minimo per t = 5, ascissa del vertice della parabola. 
La funzione 
g(t) - |l-t| + |2-t| 4- 12-t 


può’ essere riscritta come 


g(t) 


C 15-3t 
13-t 
t + 9 
V_3t-15 


per t< 1 
per 1 ^t<2 
per 2 ^ t < 12 
per t^ 12 


il minimo si ottiene (vedi fig.24) per t = 2. 

In analogia con la mediana, che divide i dati ordinati in due grup¬ 
pi ugualmente numerosi, si possono definire altri indici di posizio¬ 
ne. Per primo, secondo, terzo QUARTILE (q,, q 2 , q 3 ) si intendo¬ 
no i tre valori che dividono i dati in quattro tronconi di pari di¬ 
mensione, una volta che i dati siano disposti in ordine non decre¬ 
scente. 

In riferimento alla Tab. 3.2, il 1’, 2’, 3’ quartile corrispondono 
rispettivamente al 7’, 14’, 21’ valore dell’elenco: 

q, = 5.1 ; q 2 = 6.5 ; q 3 = 7.1 

il secondo quartile (q 2 ) e’ la mediana (A). 

Generalizzando ulteriormente, si possono definire analoghi in¬ 
dici di posizione (quintili, decili, percentili) di un certo interesse 
se i dati sono molto numerosi. 

Definiamo come MODA (o VALORE NORMALE) il dato che 
compare piu’ spesso nell’insieme dei dati. Se il carattere in esame 
e’ di tipo qualitativo o quantitativo discreto, la moda sara’ appun¬ 
to la caratteristica o il valore cui corrisponde la massima frequen¬ 
za; nel caso di un carattere continuo (o, comunque, registrato per 
raggruppamenti), avra’ senso parlare di CLASSE MODALE in ri¬ 
ferimento all’intervallo cui compete la maggior frequenza (suppo- 
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nendo di considerare intervalli della stessa ampiezza). Nell’esem¬ 
pio dei voti, la moda e’ 8. 

Oltre la media aritmetica, sono definite altre medie utilizzate per 
specifici problemi: 

si definisce MEDIA GEOMETRICA dei numeri positivi x l5 x 2 , .. 
,x n quel numero G che, sostituito ai dati, ne lascia invariato il 
prodotto. 

Per definizione: 

x^x, ... *x u = G-G ... G = G" 
da cui 


n . _ 

G = v'Xj-x,- ...-x„ 


Il calcolo di G viene normalmente eseguito utilizzando i logaritmi: 
1 

log G = — [log Xj 4- log x-, + ... + log xj 
n 

(il logaritmo delle media geometrica e’ uguale alla media aritmeti¬ 
ca dei logaritmi dei dati). 

Sempre in relazione ad un insieme di dati positivi x p x 2 , .. ,x n , 
ricordiamo la definizione di altre due medie (raramente utilizzate): 

- MEDIA QUADRATICA : 

q ■ V 

- MEDIA ARMONICA : 

n 

A = - 

1 1 1 

— + — + ... + — 

X, x 2 x n 

E’ possibile dimostrare che, se x,>0 v i, e’ 

A^G^M^Q. 

Esercizio. Dati i valori 3; 3; 5; 7; 10, verificare che A = 4.51; 
G=5.01; M = 5.6; Q = 6.2. 

Esercizio. Verificare che, se gli n dati sono tutti uguali ad uno 
stesso numero h>0, e’ A = G = M = Q = h. 

Esercizio. Dimostrare che A< G ^Q nel caso di due dati po¬ 
sitivi. 


Soluzione 

1. E’ (a - b) 2 ab ^ 0 

[(a + b) 2 - 4ab] ab ^ 0 
(a + b) 2 ab_^ 4a 2 b 2 
(a + b) \ ab ^ 2ab 


2ab 

v ab ^ - 

a + b 


2 


1 

-1- 

a 


1 

b 


cioè’ G(a,b) ^ A(a,b) 

2. E’ (Va - vb) 2 ^ 0 

a + b - 2 V ab_ ^ 0 
a + b ^ 2 Vab 
a + b > vab 


quindi M(a,b) ^ G(a,b) 
(a - b) 2 

3. E’ --- $5 0 


a 2 4- b 2 - 2ab 


> 0 


2a 2 + 2b 2 - a 2 - b 2 - 2ab ^ 0 


a 2 -t- b 2 a 2 + b 2 + 2ab 


2 

a 2 + b 2 




m 


o 


4 


a 2 + b 2 


a + b 


ovvero Q(a,b) ^ M(a,b) 


Chi si avvicina per la prima volta alla statistica sara’, a questo 
punto, un po’ disorientato: si voleva sintetizzare un insieme di da¬ 
ti e ci si trova di fronte ad una miriade di indici. Il fatto e’ che 
ciascuno degli indici introdotti evidenzia certe caratteristiche dei 
dati; alcune informazioni potranno risultare piu’ utili, altre meno. 
Se, tornando all’esempio del calzaturificio, decidete di produrre 
inizialmente uno stock di scarpe di un’unica misura (tanto per sag¬ 
giare il gradimento del pubblico), vi interesserà’ la moda piu’ della 
media o della mediana. 
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Un esempio 


Concludiamo questo paragrafo con un’applicazione della media 
geometrica (l’esempio e’ preso da [ 12 ], pag. 16): 

La popolazione di una citta’ negli anni dal 1980 al 1984 e’ 


Anno 

popolazione 

1980 

40000 

1981 

42000 

1982 

50400 

1983 

57960 

1984 

63756 

i tassi di 

incremento annuali sono: 

80-81 

42000/40000 = 1.05 (a,) 

81-82 

50400/42000 = 1.20 (a,) 

82-83 

57960/50400 = 1.15 (a~) 

83-84 

63756/57960 = 1.10 (a 4 ) 


Definiamo come ”tasso medio di crescita” nel quadriennio quei 
tasso ”t” che applicato ripetutamente a partire dal 1980 riproduca 
lo stesso ammontare della popolazione nel 1984. ”t” e’ la media 
geometrica dei tassi di incremento annuali: 

t = \ LÒ5 x 1.2 x 1.15 x ì.l = Vi.5939 = 1.1236091 


infatti: 40000 x t = 44944 
44944 x t = 50500 
50500 x t = 56742 
56742 x t = 63756 

In effetti. Pgi — ajPgg, Pg 2 = ^ 2 P«i> Ps?~^.iP 82 * Pg 4 = a 4 P 33 . Sosti¬ 
tuendo, P g4 = a 4 a 3 a ; a,P S0 . Per definizione di di ”t”, deve essere 
P s . = tP S0 ; p s 2 = tP 81 ecc. in modo che P 84 = t-t-t-t-P g0 . Quindi 


t 4 = a,a 2 a 3 a 4 e t = Va 1 a 2 a 3 a 4 . 


3.3 INDICI SINTETICI DI VARIABILI 
SEMPLICI; INDICI DI DISPERS I ONE _ 

I voti dei compiti di matematica di Massimo sono 6 , 8 , 6 ,4; quelli 
dei compiti di Elisabetta 6 , 6 , 6 , 6 . Per entrambi e’ M = 6 , Media¬ 
na = 6 , Moda = 6 ma, mentre Elisabetta e’ sempre sufficiente- 
mente preparata, il rendimento di Massimo e’ discontinuo. 
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Variabilità 
dei dati 


Scostamento 

semplice 

medio 


Scostamento 

quadratico 

medio 


La DISPERSIONE (o VARIABILITÀ’) di un insieme di dati 
e’ senz’altro un aspetto importante e risulterà’ utile calcolare dei 
valori che ne misurino l’entità’ (INDICI DI DISPERSIONE). 

Una prima informazione e’ data dal CAMPO DI VARIAZIO¬ 
NE (V) definito come differenza tra il maggiore e il minore dei dati; 
se 


x’ - min Xj e 
x” - max Xj e’ 
V = x” - x’ . 


Generalmente non si tratta di un valore molto indicativo (un so¬ 
lo dato eccezionalmente grande o eccezionalmente piccolo ”dila- 
ta” il campo di variazione); può’ rivelarsi utile se V risulta parti¬ 
colarmente ridotto rispetto all’ordine di grandezza dei dati. 

Il campo di variazione può’ essere espresso attraverso la coppia 
(x’,x”) in modo da fornire anche indicazioni ”di posizione”: tutti 
i dati sono compresi nell’intervallo (x’,x”). 

Piu’ interessante e’ lo SCOSTAMENTO SEMPLICE MEDIO 
definito come: 



1 * 


- M\ 


che dice di quanto, mediamente, i dati si discostano dalla loro me¬ 
dia aritmetica. 

In relazione all’ultimo esempio si ha: 

- Voti di Massimo (M = 6 ): 

(x’,x”) = (4,8) ; V = 8-4 = 4 

a’ = ~ [|6 - 6( +18 - 61 +16 - 6| + \4 - 6|] = 1 

- Voti di Elisabetta (M = 6 ): 

(x’,x”) = ( 6 , 6 ) ; V = 6-6 = 0 
ct ’ = 0 . 

Piu’ comodo da calcolare e di uso piu’ comune e’ lo SCOSTA¬ 
MENTO QUADRATICO MEDIO (o DEVIAZIONE 
STANDARD): 

a = Ay ■—- ^ ( Xi - 

i 

il cui significato e’ analogo a quello di o’ (anche se, nel computo 
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Varianza 



di a, grossi scarti hanno ”maggiore peso”). II quadrato dello sco¬ 
stamento quadratico medio (a 2 ) e’ la VARIANZA. 

Per quanto riguarda i voti di Massimo abbiamo: 



la somma dei valori deH’ultima colonna e’ 8 per cui 



= V 2 = 1.414. 


Ovviamente a (Voti di Elisabetta) = 0. 
Esercizio. Dimostrare che 


? (Xj - M) 2 = 2 x 2 - nM 2 


Soluzione: 


2 (X; - M) 2 = 2 (x 2 i - 2MXi + M 2 ) = 

= 2 x 2 j - 2M 2 Xj + nM 2 = 

i i 

= 2 x 2 ; - 2M(Mn) + nM 2 = 2 x^"- nM 2 . 


Questa identità’ consente di calcolare a come 


o 



— [2 x 2 j - nM 2 ] 
n ' 



- M 2 


una formula alternativa assai utile dato che non e’ necessario di¬ 
sporre preventivamente della media aritmetica (si evita di passare 
in rassegna i dati due volte e quindi di doverli memorizzare in un 
vettore). 


Ricalcoliamo o (Voti di Massimo) con la nuova formula: 


^ : 6 

8 

v — 

6 

4 

x 2 - : 36 

64 

36 

16 




X. x 2 j = 152 ; M 2 = 36 ; n = 4 per cui 



= V2 


X-; 

! x, - M | 

(x; - M ) 2 

3.8 

2.389 

5.707 

3.9 

2.289 

5.240 

4.5 

1.689 

2.853 

4.5 

1.689 

2,853 

5.0 

1.189 

1.414 

5.0 

1.189 

1.414 

5.1 

1.089 

1.186 

5.5 

0.6S9 

0.475 

5.5 

0.689 

0.475 

5.6 

0.589 

0.347 

5.8 

0.389 

0.151 

5.9 

0.289 

0.084 

6.1 

0.089 

0.008 

6.5 

0.311 

0.097 

6.5 

0.311 

0.097 

6.6 

0.411 

0.169 

6.6 

0.411 

0.169 

6.8 

0.611 

0.373 

6.8 

0.611 

0.373 

7.0 

0.811 

0.658 

7.1 

0.911 

0.830 

7.5 

1.311 

1.719 

7.5 

1.311 

1.719 

8.0 

1.811 

3.280 

8.0 

1.811 

3.280 

8.0 

1.811 

3.280 

8.0 

1.811 

3.280 

Totale 

_ 

28.511 

41.531 


Campo di variazione = 8 - 3.8 = 4.2 
Scostamento semplice medio a' = 28.511/27 = 1.056 
Deviazione standard a — V41.531/27 = 1.24 

Differenza interquartile = 7.1 - 5.1 = 2 







Ecco un breve programma Basic per il calcolo di M e a : 

10 INPUT''QUANTI DATI 1 ';N 
20 FOR I = 1 TO N 
30 INPUT-'DATO..' f ;X 
40 S=S+X:Q=Q+X*X 
50 NEXT I 

60 M=S/N:D=SQR(Q/N-M*M) 

70 PRINT M,D 

Differenza Altro indice di dispersione e’ la DIFFERENZA INTER- 

interguartile QUARTILE 

Qi - Qi 

che, come il campo di variazione, può* essere espressa dalla cop¬ 
pia Il significato e’ semplice: meta’ dei dati cadono in que¬ 

sto intervallo. 

Tornando all’esempio della Tab. 3.2 (n = 27 voti; M = 6.189), cal¬ 
coliamo gli indici di dispersione: 


Ricapitolando: 


CLASSE 

r SEZ. A 


M 

= 6.189 

media aritmetica 

X 

= 6.5 

mediana 

MD 

- 8 

moda 

(x\x”) 

= (3.8 , 8) 

campo di variazione 

V 

- 4.2 

campo di variazione 

(qi.q 3 ) 

= (5.1 , 7.1) 

differenza interquartile 

q 3 - qi 

= 2 

differenza interquartile 

a’ 

= 1.056 

scostamento semplice medio 

o 

= 1.24 

scostamento quadratico medio 

a 2 

= 1.5376 

varianza 


La statistica mostra che il compito proposto può’ ritenersi ade¬ 
guato (M e X oltre ”la sufficienza”); la classe e’ abbastanza omo¬ 
genea: anche se V c’ piuttosto elevato, a’ e a mostrano che, media¬ 
mente, ci si discosta di 1 voto dalla sufficienza. 

Nella classe parallela, la 1’ B, gli studenti hanno usufruito di 
un’ora di lezione in piu’ dedicata al ripasso. I risultati, relativi allo 
stesso compito valutato con gli stessi criteri, sono riportati nella 
Tabella 3.4 
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li programma 
STATISTICA 


Tabella 3.4 


3.8 

4.1 

4.5 

4.9 

5.5 

5.8 

6.0 

6.1 

6.5 

6.6 

6.6 

6.6 

6.8 

6.8 

6.8 

6.8 

6.9 

7.0 

7.0 

7.0 

7.0 

7.0 

7.0 

7.1 

7.1 

7.1 

7.6 


i- 

Per questa classe abbiamo: 

CLASSE 

1’ SEZ. B 


M 

= 6.37 

media aritmetica 

X 

= 6.8 

mediana 

MD 

= 7 

moda 

(x\x”) 

- (3.8 , 7.6) 

campo di variazione 

V 

= 3.8 

campo di variazione 

(qi>q 3 ) 

= (6,7) 

differenza intcrquartiìe 

q 3 - di 

= 1 

differenza intcrquartiìe 

< j ' 

= 0.76 

scostamento semplice medio 

a 

= 0.97 

scostamento quadratico medio 

2 

a- 

= 0.94 

varianza 


Dal confronto delle due statistiche si nota come la situazione della 
1 ’ B sia decisamente migliore: voti complessivamente piu’ alti e mag¬ 
giore omogeneità’. Sembra proprio che la lezione aggiuntiva sia 
servita. 

Naturalmente, trattandosi di valutazioni scolastiche, il ”6” ri¬ 
sulta un valore particolare poiché’ discrimina tra prove sufficienti 
ed insufficienti; nella 1’ A sono pienamente sufficienti (x, ^6) 15 
studenti su 27 (circa il 56%), nella 1’ B 21 su 27 (circa il 78%). 

Il programma STATISTICA, che si propone scopi puramente 
didattici, consente di calcolare gli indici statistici descritti in rela¬ 
zione ad un insieme di dati positivi che possono essere registrati 
su disco (o nastro) e da li’ richiamati. 

Esempio di Output del programma STATISTICA 

(Dati: 7, 3, 10, 10, 10, 9, 9) 


Risultati : 

Media Aritmetica . 8.28571429 

Media Geometrica . 7.76426976 

Media Quadratica . 8'.61891608 

Media Armonica. 7.01112878 

Moda . 10 
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Frequenza vai.modale . 3 

Mediana . 9 

Primo Quartile . 7 

Terzo Quartile . 10 

Campo di variazione . 7 

Minimo Xi . 3 

Massimo Xi . 10 

Differenza interquar. 3 

Scost. semplice medie . 1.87755102 

Scarto quadr. medio . 2.3733211 

Varianza . 5.63265305 


Campionamento 

casuale 

semplice 


Campioni 

stratificati 


3.4. INFERENZA STATISTICA 


Nei paragrafi precedenti ci siamo occupati della rappresentazio¬ 
ne e della sintesi di un insieme di dati (STATISTICA DESCRIT¬ 
TIVA) indipendentemente dal fatto che questi provenissero dall’in¬ 
tera popolazione o da una sua parte. 

E’ facile rendersi conto che uba rilevazione totale e’ molto spes¬ 
so da escludere: non appena l’universo statistico supera una certa 
dimensione, diviene praticamente impossibile rilevare i dati da tutte 
le unita’ statistiche. Inoltre una rilevazione su vasta scala può’ es¬ 
sere assai costosa sia in termini monetari (si pensi di valutare l’ef¬ 
ficienza di un farmaco) sia in termini di tempo (nei sondaggi pre¬ 
elettorali, i risultati servono subito); per non parlare dei rilevamenti 
”distruttivi” in cui l’unita’ esaminata risulta poi inservibile (per 
saggiare la resistenza dei tondini metallici utilizzali in edilizia, li 
si sottopone a ”prove a trazione fino a rottura”). 

Quando si e’ costretti ad un rilevamento parziale, si pone un pro¬ 
blema: fino a che punto i risultati dedotti dal campione (normal¬ 
mente assai ridotto rispetto all’universo statistico) possono esten¬ 
dersi alla totalità’? 

Di questi tipi di problemi si occupa l’INFERENZA STATISTI¬ 
CA, il cui studio richiede conoscenze matematiche piuttosto avan¬ 
zate; ci limiteremo a qualche considerazione. 

Un primo grosso problema consiste nella determinazione del cam¬ 
pione da utilizzare: come sceglierlo e quanto grande ? 

Per la scelta, si ricorre spesso a metodi probabilistici; nel CAM¬ 
PIONAMENTO CASUALE SEMPLICE le ”k” unita’ che for¬ 
mano il campione vengono estratte casualmente dalle ”n” che co¬ 
stituiscono l’universo. 

In questo modo si evita di privilegiare, magari inconsciamente, 
certe categorie e di lavorare su insiemi non rappresentativi. Se la 
popolazione e’ molto variabile, ma può’ essere suddivisa in sot- 
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toinsiemi singolarmente piu’ omogenei, si preferisce operare sepa¬ 
ratamente sui sottoinsiemi (con estrazioni casuali) e mediare poi 
i risultati (CAMPIONAMENTI PER STRATIFICAZIONE). 


Tre esempi 


ESEMPIO 

Vi sono 1000 persone di età’ compresa tra i 20 e i 21 anni resi¬ 
denti nella citta* X; sono interessato al loro peso medio e voglio 
limitare l'indagine ad un campione di 100 persone. 

PROCEDIMENTO 1: scelgo a caso 100 persone, rilevo i dati 
e medio i pesi 

PROCEDIMENTO 2: scelgo a caso 50 ragazzi e 50 ragazze, cal¬ 
colo la media dei due gruppi e medio i risultati. 

Che succede se la popolazione e ! di 600 ragazze e 400 ragazzi? 
Nel primo caso il campione conterrà’, probabilmente, piu’ femmi¬ 
ne che maschi ed il valore ottenuto con il procedimento 1 sara’ piu’ 
attendibile di quello ricavato con il procedimento 2. 

Se pero’ so, fin dall’inizio, che l’universo e’ formato da 600 ra¬ 
gazze e 400 ragazzi, posso considerare un campione stratificato di 
60 -t- 40 ottenendo un valore ancora piu’ attendibile (si sta lavoran¬ 
do nell’ipotesi che i due insiemi siano singolarmente piu’ omoge¬ 
nei della totalità’). 

Quando si utilizzano campioni casuali, l’inferenza statistica ior- 
nisce tutta una serie di strumenti che consentono di valutare con 
quale grado di fiducia le informazioni ricavate dal campione si pos¬ 
sano riferire a tutta la popolazione. Viceversa, fissato il margine 
eli errore che si e’ disposti a tollerare, si potrà’ determinare qual 
e’ la dimensione minima del campione necessario. 

Concludiamo con tre esempi (il 2’ e il 3’ sono presi da [14] pag. 
53 e pag. 58) che diano un’idea dei problemi tipici dell’inferenza 
statistica c del tipo di approccio per risolverli; gli studenti degli ul¬ 
timi anni della secondaria superiore troveranno in [4] (volume 3) 
una trattazione molto piu’ esauriente. 


ESEMPIO 1 

Vorrei sapere quanti, tra i 6 vicini di casa, possiedono un televi¬ 
sore a colori. 

Decido di chiedere a tre di loro e, dopo averli numerati da 1 a 
6 , lancio (re dadi ed intervisto le persone corrispondenti ai numeri 
ottenuti. Degli intervistati, due mi dicono di avere un apparecchio 
a colori, il terzo no. 

Come stanno le cose in relazione a tutto il vicinato? 


SOLUZIONE 

L’universo statistico c’ costituito da 6 unita’. Il campione, di di¬ 
mensione 3, viene determinato attraverso una ESTRAZIONE BER- 


115 



NOULLIANA (o CON RIPETIZIONE): dopo l’estrazione, l’uni¬ 
ta’ statistica sorteggiata torna a far parte della popolazione e può’ 
essere nuovamente estratta (con un universo statistico di dimen¬ 
sioni tanto ridotte, questo modo di procedere appare ridicolo poi¬ 
ché’ la stessa persona può’ essere intervistata due volte o addirit¬ 
tura tre). 

Identifichiamo i vicini con 6 palline contenute in una scatola: 
i possessori di televisori a colori corrispondano a palline bianche, 
gli altri a palline nere. 

Il problema diviene: 

”Se una scatola contiene 6 palline tra bianche c nere e ne vengono 
estratte 3 con reimbussolamento ottenendone 2 bianche e 1 nera, 
che si può’ dire del contenuto della scatola?” 

Indichiamo con x il numero di palline bianche presenti nella sca¬ 
tola e con Hx l’evento 

Hx = {la scatola contiene x bianche e 6-x nere). 

Prima di eseguire le estrazioni, sappiamo solo che x può’ essere 
un qualsiasi intero da 0 a 6: l’assenza di informazione equivale al¬ 
l’ipotesi P(Hx)=l/7 per ogni x. 

Posto 

A = {vengono estratte 2 bianche e 1 nera) 
e’ 


P(A/Hx) = 




6 - x 


x 2 (6 - x) 
36 


si tratta infatti della probabilità’ di avere 2 successi (successo = 
estrazione di una bianca) su 3 prove in uno s.p.r. di parametro in¬ 
cognito p = x/6 (vedi 4.5). Giustamente P(A/H 0 ) = P(A/H 6 ) = 0 
dato che A non può’ verificarsi se le palline sono tutte bianche o 
tutte nere. 

Il problema consiste nel valutare P(Hx/A) in quanto sappiamo 
che A si c’ verificato mentre non conosciamo x. Ci viene in aiuto 
la formula di Bayes (di cui si parlava nel problema n.7 del 2’ capi¬ 
tolo); poiché’ Hx x = 0,l,..,6 costituisce una partizione di 1), e’ 


P(Hx/A) = 


E (A/HJ-H (Hx) 

6 

Z P(A/H,).P(H,) 

j”0 


Calcoliamo il denominatore: 


6 


E 

j=0 


P(A/Hj) P(Hj) = ^ 

j=o 


r (6-j) ^ 

36 ' 7 



I 

36 


• y [0 + 5 + 


16 + 27 + 32 + 25] = — 

12 


quindi 


P(Hx/A) = 


x 2 (6-x) 

36 


7 


12 x 2 (6-x) 

5 ~ 105 


Per x = 0,...,6 abbiamo: 


X = 

P(Hx/A) = 

0 

0 

1 

5/105 = 0.05 

2 

16/105 = 0.15 

3 

27/105 = 0.26 

4 

32/105 = 0.30 

5 

25/105 = 0.24 

6 

0 


Il valore piu’ alto si ottiene per x = 4 (come ci si aspettava) an¬ 
che se la probabilità’ associata non e’ poi elevatissima. La tabella 
fornisce pero’ informazioni piu’ utili, ad esempio: 

P[3 ^ x 5] = P(3) + P(4) + P(5) = 0.8 ; 
posso ora affermare, con un margine di sicurezza dell’80%, che 
il numero di vicini che possiedono un televisore a colori e’ com¬ 
preso tra 3 e 5 mentre prima dell’esperimento statistico tale proba¬ 
bilità’ era pari a 3/7=43%. 


ESEMPIO 2 

Una fabbrica produce componenti elettronici in confezioni da 
100 pezzi che provengono da due macchine MI e M2 e precisamente: 
60 pezzi da MI, 40 da M2. 

A causa di un temporaneo malfunzionamento di M2, alcune con¬ 
fezioni contengono 40 circuiti difettosi. Ne controllo una verifi¬ 
cando 5 pezzi scelti a caso: risultano perfettamente funzionanti. 

Posso essere sicuro, al 95 %, della bontà’ di tutti i cento pezzi 
della confezione in esame? 

SOLUZIONE 

Interessa calcolare la probabilità’ che una confezione difettosa 
superi il controllo ovvero la probabilità’ che, nelle 5 estrazioni, siano 
sempre usciti pezzi funzionanti quando la confezione ne contiene 
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40 direttosi. 

La probabilità’ richiesta e’ (estrazioni SENZA reimbusso- 
lamento): 


60 59 58 57 56 

100 ’ 99 ' 98" ’ In ' ~96 


0.07254.. 


7.2 % 


P e’ la "probabilità’ di sbagliarsi” cioè’ la probabilità* di consi¬ 
derare valida una confezione difettosa. Poiché’ 7.2%>5%, non 
posso essere sicuro al 95% (ma solo al 92.8%). 

Se avessi eseguito le 5 estrazioni CON reimbussolamento, la pro¬ 
babilità’ di "errore” sarebbe stata: 

P’ = j-) = 0.078 = 7.8 % . 

\100 ) 


Supponendo di eseguire estrazioni con reimbussolamento, quanti 
pezzi devo controllare (come minimo) per essere sicuro della bori¬ 
ta’ della confezione al 95%? E al 99%? 

Sia ”n” il numero di circuiti da controllare; ”n” e’ il piu’ pic¬ 
colo intero per cui: 

(0.6) n ^ 0.05 per il 95% 
e (0.6) n ^ 0.01 per il 99% 


n 

(0.6)“ 

1 

0.60000 

2 

0.36000 

3 

0.21600 

4 

0.12960 

5 

0.07776 

6 

0.04666 *** 

7 

0.02799 

8 

0.01680 

9 

0.01008 

10 

0.00605 *** 


Quindi 6 pezzi per essere sicuri al 95% e 10 pezzi per il 99%. 

11 risultato si ottiene piu’ rapidamente passando ai logaritmi. Dove 
essere: 

log (0.6) n ^ log 0.05 per il 95 % 
e 

log (0.6)“ ^ log 0.01 per il 99 % 
da cui, rispettivamente: 

n • log 0.6 ^ log 0.05 




Campione 

distorto 


e n • log 0.6 < log 0.01 
poiché’ log 0.6 < 0 e’ 


n > - 


log 0.05 


= 5.864 


log 0.6 

per il 95% : il minimo intero e’ 6; mentre 


log 0.01 
log 0.6 


9.0151 


per il 99% : il minimo intero e’ 10. 

4 


ESEMPIO 3 

Si vuol sapere qual e’ staio il numero medio di giorni di degenza 
dei pazienti ricoverati presso un certo ospedale nel 1986. 

Marinella e Sergio, incaricati dell’indagine, procedono in modo 
diverso: Marinella, procuratasi un elenco delle persone ricoverate, 
ne estrae un campione su cui lavorare; Sergio sceglie a caso un cer¬ 
to numero di date del 1986 e prende in considerazione le persone 
che risultano ricoverate in quei giorni. Chi dei due otterrà’ risulta¬ 
ti piu’ attendibili? 

SOLUZIONE 

Il campione utilizzato da Sergio risulta essere un CAMPIONE 
DISTORTO; a parità’ di dimensioni, il valore da lui calcolato Sa¬ 
ra’, probabilmente, piu’ elevato di quello valutato da Marinella. 

Con il metodo di Sergio infatti, le persone costrette a lunghe de¬ 
genze hanno maggiore probabilità’ di entrare a far parte del cam¬ 
pione contribuendo ad aumentare il valore medio. 

Consideriamo un ”caso limite” in cui, di 100 pazienti, 50 siano 
stati ricoverati dal 1’ gennaio al 10 gennaio e i rimanenti 50 dall’ 
11 gennaio al 31 dicembre (una situazione davvero poco invidiabile). 

Mentre Marinella est.rarra’ un campione che, in linea di massi¬ 
ma, conterrà’ meta’ clementi del primo gruppo e meta’ del secon¬ 
do, quello di Sergio sara’ costituito quasi esclusivamente da perso¬ 
ne del secondo gruppo, polche’ le date scelte casualmente cadran¬ 
no quasi tutte oltre il 10 gennaio. 

Nella realta’, il carattere in esame risulta molto variabile e con¬ 
verrebbe far uso di campioni stratificati, studiando separatamente 
i vari reparti in diversi periodi dell’anno. 
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CAPITOLO 4 


CALCOLO COMBINATORIO 
E DIST RIBUZIONE BIN O MI ALE 


Il calcolo combinatorio si occupa della formazione e del conteg¬ 
gio dei raggruppamenti che si possono realizzare utilizzando gli ele¬ 
menti di un insieme finito (o di piu’ insiemi finiti). 

Poiché’ si tratta di una parte affrontata in quasi tutti i libri di 
lesto, ci limiteremo agli argomenti principali. Considereremo un 
insieme A formato da n elementi. 


4.1 DISPOSIZIONI SEMPLICI 


Si definiscono "disposizioni semplici di classe k” (l ^ k^n) tut¬ 
ti i raggruppamenti di k elementi distinti di A tali che ogni rag¬ 
gruppamento differisca dagli altri 

- per la natura degli elementi (cioè’ per avere almeno un elemento 
diverso) 

oppure 

- per l’ordine degli elementi (cioè' per avere gli elementi disposti 
diversamente). 

Es.l Sia 
A = {a,b,c,dj 

quindi n = 4; costruiamo le disposizioni semplici di classe 3 (k = 3): 


abe 

abd 

aed 

bed 

acb 

adb 

adc 

bdc 

bac 

bad 

cad 

ebd 

bea 

bda 

eda 

cdb 

cab 

dab 

dac 

dbc 

eba 

dba 

dea 

deb 
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le terne in colonne diverse differiscono per almeno un elemento, 
quelle nella medesima colonna differiscono per l’ordine. 

Il numero delle disposizioni semplici di n oggetti della classe k 
si indica con D n k ed e’ 

D nk = n-(n-l)-(n-2)-...-(n-k + 1) (*) 

cosi’ D 4 3 = 4-3-2 = 24. Si noti che il numero di fattori e’ k per 
cui D„ k e’ dato dal prodotto di k fattori naturali decrescenti a par¬ 
tire dan(D 52 = 5-4 = 20; D 53 — 5-4-3 = 60; D 73 — 7-6-5 = 
210 ). 

Dimostrazione della (*). Supponiamo di aver elencato le dispo¬ 
sizioni di classe k- 1 ; per ottenere quelle di classe k basterà’ aggiun¬ 
gere ad ogni raggruppamento uno degli n-(k-l) = n-k + 1 elementi 
non ancora utlizzati. Dunque D nk = (n-k + 1)-D n k _,; poiché’ D n , 
= n avremo; 

D nJ = n 

D n 2 = (n-2+ 1)- D n , = (n-l) -D n j 
D n 2 = (n-3 + 1)- D n2 = (n-2) • D n 2 

Dn.k = (n-k + 1) • D n k .j ; 

uguagliando il prodotto dei primi membri con quello dei secondi 
membri e sopprimendo i fattori comuni si ottiene la (*). 


4.2 DISPOSIZIONI CON RIPETIZIONE 


La definizione c’ analoga a quella data al punto precedente con 
la differenza che lo stesso oggetto può’ apparire nel raggruppamento 
piu’ volte cioè’ da 0 a k volte e che k può’ essere maggiore di n. 
Es.2 Sia 
A= {a,b,c,d} 

quindi n = 4. Le disposizioni con ripetizione di classe 2 (k = 2) sono: 


aa 

ab 

ac 

ad 

ba 

bb 

bc 

bd 

ca 

cb 

cc 

cd 

da 

db 

de 

dd 


Es.3 Sia A= {x,y ì quindi n = 2. Le disposizioni con ripetizione 
della classe 3 (k = 3) sono: 


xxx xxy xyx xyy 



yxx yxy yyx yyy 


(in questo esempio e’ k>n). 

Il numero delle disposizioni di n oggetto della classe k si indica 
con D’ nk e si ha: 

= n k - 

In effetti, supponendo di aver elencato le disposizioni con ripe¬ 
tizione di classe k-1, quelle di classe k si otterranno aggiungendo 
a ciascuna delle precedenti uno degli n oggetti. Quindi 


D’nje = n-D’ 


n,k-l > 


poiché’ 



- n, 


avremo 

D’ a2 = n-D’„, = n-n = n 2 , D’ n _ ? = n-D’ n2 = n-n 2 = n 3 
e, in generale, 

D\ k = n k . 

Es.4 Quante sono le diverse colonne della schedina del totocal¬ 
cio? Sono tante quante le disposizioni con ripetizione di 3 oggetti 
(1,X,2) delle classe 13: D\ 13 = 3 13 = 1594323. 

Es.5 Una scatola contiene 3 palline di colore differente. Si ese¬ 
guono 2 estrazioni CON reimbussolamento. Calcolare la probabi¬ 
lità’ dell’evento 

A= {le palline estratte sono di colore diverso}. 

Soluzione 1. Qualunque sia l’esito della prima estrazione, l’e¬ 
vento A si verifica se alla seconda si ottiene un colore diverso dalla 
prima. Vi sono 3 casi possibili di cui 2 favorevoli: P(A) = 2/3. 

Soluzione 2. I casi possibili sono le disposizioni con ripetizione 
(reimbussolamento) di 3 oggetti (le palline) a gruppi di 2. Il loro 
numero e’ D’ 3 , = 9. I casi favorevoli sono le disposizioni semplici 
(colore diverso) in numero di D 3 2 = 6. Allora P(A) = 6/9 = 2/3. 


43 PE RMUTAZIONI _ 

Le disposizioni semplici di n oggetti della classe n vengono dette 
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permutazioni; il loro numero si indica con P n . Per definizione e’ 

Pj, — D nn = n(n-l)(n-2)... (n-n+1) = n(n-l)(n-2) ... 3-2-1 = n! 

dove si e’ indicato con n! (”enne fattoriale”) il prodotto degli in¬ 
teri da 1 a n. 

Dato che nelle permutazioni compaiono tutti gli oggetti, una per¬ 
mutazione differisce da un’altra per l’ordine. 

Es.ó Sia 
A= {a,b,c}, 
le permutazioni sono 

abc bac cab acb bea eba , 
p 3 = 3! = 3-2-1 =6. 


4.4 ANAGRAMMI 


Se due diverse parole sono formate dalle stesse lettere alfabeti¬ 
che, si dice che una e’ un anagramma dell’altra (es. CASO CAOS 
COSA; TORTA TROTA ROTTA RATTO). Ci proponiamo di cal¬ 
colare il numero di anagrammi distinti (anche senza significato) di 
una parola di n lettere (compresa la parola stessa). 

Se le lettere sono tutte diverse, la risposta e’ P n = n! poiché’ 
gli anagrammi coincidono con le permutazioni delle n lettere. Co¬ 
si’ gli anagrammi della parola ALT sono 6: 

ALT, ATL, LAT, LTA, TAL, TLA. 

Se una lettera compare piu’ volte, il numero di anagrammi e’ 
inferiore a n! in quanto andranno identificate quelle parole che si 
ottengono dalle permutazioni della lettera ripetuta e questo per ogni 
posizione fissata nell’ambito della parola. Cosi’ in AMACA (n = 6), 
uno scambio tra la prima e la seconda ”A” non origina una nuova 
parola: poiché’ le permutazioni delle tre ”A” sono in numero di 
3! =6, il numero di anagrammi distinti sara’ 6!/3! = 120. 

Il ragionamento si estende in modo analogo al caso di piu’ lette¬ 
re ripetute: se in una parola di n lettere 
una lettera e’ ripetuta k, volte 
un’altra lettera e’ ripetuta k 2 volte 
una terza lettera e’ ripetuta k, volte 

il numero di anagrammi e’ n!/(k,!-k 2 !-k 3 ! ■...). 

Es.7 Nella parola ORO di 3 lettere, la lettera O compare 2 volte; 
il numero di anagrammi e’ 3!/2! = 3 e sono ROO, ORO, OOR. 
Es.8 Nelle parola PAPA di 4 lettere 

- la lettera P compare 2 volte 

- la lettera A compare 2 volte; 




il numero di anagrammi c’ 4!/(2! -2!) = 6 : AAPP, APAP, APPA, 
PPAA, PAPA, PAAP. 

Es. 9 "MATEMATICA”: n= 10; due M, tre A, due T. Il nu¬ 
mero di anagrammi e’ 10 !/( 2 ! -3!- 2 !) = 151200. 

Di particolare interesse e’ il caso di una parola di n lettere in cui 
una stessa lettera compare k volte e una seconda lettera le rima¬ 
nenti n-k volte (come in AAABB n = 5, k = 3). Il numero di ana¬ 
grammi e’: 


k! (n - k)! 


valore che viene detto COEFFICIENTE BINOMI ALE c che si in¬ 
dica con il simbolo da leggersi ”n su k” o ”n sopra k”. Se 

k = 0 o k = n, le n lettere sono tutte uguali c il numero di anagram¬ 
mi si riduce a 1; ponendo 01 = 1, la formula vale per O^k^n. 

Nello sviluppo di (a + b) n , il coefficiente di a k b u k risulta pari al 
numero di anagrammi di una parola formata da k lettere ”a” e 
n-k lettere ”b”; ad esempio: 

(a + b) 3 = (a + b)(a + b)(a + b) = (aa + ab-t-ba + bb)(a + b) = 

= aaa + aba + baa + bba + aab + abb + bab + bbb = 

= a 3 + 3a 2 b + 3ab 2 + b 3 ; 
in generale: 


(a + b) n = 


-£ (E) 


a k b n 


k. 


da qui la denominazione di coefficiente binomiale attribuita a 



4.5 DISTRIBUZIONE BINOMIALE 


In uno schema delle prove ripetute di parametro p (0<p< 1), 
la probabilità’ che si verifichino k successi in n prove ( 0 <k^n) 
e’ data da: 


b(k,n;p) = Q p k q nk 


dove q= 1 -p. 


In effetti, l’evento 

A= [k successi in n prove) 

si verifica se si realizza una qualsiasi sequenza del tipo SFFSS...FFS 
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dove S (successo) compare k volte e F (fallimento) n-k volte. Ogni 
sequenza ha probabilità’ p k q nk , mentre il numero delle diverse se¬ 
quenze c’ dato dagli anagrammi della ”paro!a” SSS..SFFF..F 



Per Padditivita’: 


p ( A ) = p k q nk - 

Es.10 Un dado regolare viene lanciato 10 volle; calcolare la pro¬ 
babilità’ che la faccia ”6” compaia esattamente 3 volte. 

Soluzione. Detto ''successo” l’uscita del 6, e’ p= 1/6, q = 5/6, 
n= 10, k = 3. La probabilità’ richiesta e’ 


Es.11 Una moneta regolare viene lanciata 5 volte; calcolare la 
probabilità’ che la faccia testa compaia 0 volte, 1 volta, 2 volte,..., 
5 volte. 

Soluzione. Detto ''successo” l’uscita ”T”, e’ p = q = 1/2 cn = 5. 
Per k = 0,l,..,5 sara’ 

™>-C) ■(!)'• (!)“-(() (F 


Calcolando: 

P(0) = 1/32; P(l) = 5/32; P(2) = 10/32; P(3) = 10/32; P(4) = 5/32; 
P(5) = 1/32. 

Si noti che P(0) + P(1 )+...+ P(5) = 1 e, in generale, 

n n 

b(k,n;p) = 

k = t> k=0 



p k q n k = (p + q) B = 1" = 1. 


4.6 COMBINAZIONI SEMPLICI 


Dato l’insieme A di n elementi, si definiscono COMBINAZIO¬ 
NI SEMPLICI di classe k i raggruppamenti di k elementi distinti 
di A tali che ogni raggruppamento differisca dagli altri per la na¬ 
tura degli elementi. 

Es/12 Sia 
A= (a,b,c,d,e}, 

quindi n = 5. Costruiamo le combinazioni di classe 3: 


abe abd abe aed acc 



ade bcd bcc bde cde. 


Due qualsiasi combinazioni differiscono per avere almeno una 
lettera diversa (”abc” c ”acb” costituiscono due disposizioni di¬ 
stinte ma la stessa combinazione). 

Il numero delle combinazioni semplici di n oggetti della classe 
k si indica con C„ ik (quindi C }-3 = 10). 

Vale la seguente uguaglianza: 

D nk n(n-l) • ... • (n-k+1) 


in effetti, permutando in tutti i modi possibili i k oggetti di ciascu¬ 
na combinazione si ottengono le disposizioni semplici per cui 
C rMc P k = D nk (si consideri l’insieme D delle disposizioni sempli¬ 
ci e si introduca la relazione di equivalenza d~d’ se d c d’ sono 
formate dagli stessi oggetti. Ogni classe di equivalenza contiene esat¬ 
tamente P k disposizioni e l’insieme quoziente sara’ formato da 
D 1 Jk /P k classi di equivalenza). 

Applicando la (**) si ha: 



D s ,3 5-4-3 


P 3 3-2-1 

7-6-S-4-3 

- =21 

5 • 4 • 3 - 2 ■ 1 


10 


Es.13 Qual e’ il numero di diagonali di un poligono non intrec¬ 
ciato di n lati? 

Soluzione. Dopo aver numerato i vertici da 1 a n, le diagonali 
saranno tante quante le coppie non ordinate del tipo {h,k ] con h 5 * k 
eccettuate le n coppie che identificano i lati. 11 numero di diagona¬ 
li e’ C„ 2 -n = n(n-3)/2. 


Proprietà’ 
dei coeffieienli 
binomiali 


Moltiplicando per (n-k)! entrambi i termini della frazione nella 
(**) si ottiene: 

c - n! - ( n ) 

”' k k!(n-k)! W 

I coefficienti binomiali godono di diverse proprietà’, tra cui: 
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n-k 


Triangolo 
di Pascal 



La prima, di verifica immediata, si spiega intuitivamente: per 
ogni combinazione di k elementi ne rimane una costituita dagli n- 
k non utilizzati per cui C n k = C n n . k . 

Per quanto riguarda la seconda abbiamo: 

/ n-l\ + / n-l\ = (n-1)! (n-1)! = 

V k / Vk-iy k! (n-l-k)! + (k-l).'(n-k)! 

= (n-1)! I" J_1_ 1 

(k-l)J(n-l-k)! L k n-k J 

(n-1)! n n! /n\ 

~ (k-l)!(n-l-k)! ' k(n-k) = k!(n-k)! ~~ \k/ 
dove si e’ tenuto conto che m! = m(m-l)L 

I C n k anagrammi della parola AA..ABB..B formata da k lettere 
A e n-k lettere B iniziano con A oppure con B. 

Quelli che iniziano con A sono tanti quanti gli anagrammi di una 
parola di n-1 lettere di cui k-1 sono A; il loro numero e’ C n _, k _,. 

Quelli che iniziano con B sono tanti quanti gli anagrammi di una 
parola di n-1 lettere di cui k sono A; il loro numero e’ C n _, k . 
Quindi C n>k = e*,*., + C n . 1>k . 

Della seconda proprietà’ si fa uso nella costruzione del cosid¬ 
detto "triangolo di Pascal” (o "triangolo di Tartaglia”). 

Di facile verifica la terza proprietà’ che fornisce un’utilissima 
relazione ricorrente per il calcolo dei termini della distribuzione bi- 
nomiale: 

/ n \_ n! n! (n-k) /n\ n-k 

\k+l/ ” (k+l)!(n-k-l)! (k+l)-k!-(n-k)! \k/ k+1 


da cui 


b(k + l,n;p) = b(k,n;p) • 


(n-k)p 

(k+l)q 


Il programma 

CALCOLO 

COMMINATORIO 


Questa relazione e’ stata utilizzata nel programma BINOMI ALE. 

Il programma CALCOLO COMBINATORIO consente il cal¬ 
colo del numero dei vari raggruppamenti e l’elencazione degli stessi. 

Anche per modesti valori di n e k, il numero di raggruppamenti 
può’ essere molto elevato (ad esempio, le permutazioni di 10 og¬ 
getti sono 3628800), per cui conviene evitare la stampa di elenchi 
molto lunghi (soprattutto su carta) che, inoltre, richiedono tempi 


128 



non indifferenti per essere elaborati. Per le routine di generazione, 
si veda [10] e [13]. 

Esempi di Output del programma CALCOLO COMBI¬ 
NATORIO: 

1. Disposizioni semplici: n=4; k = 2 
D (4,2) = 12 

ab ac ad ba bc bd ca cb cd da db de 

2. Disposizioni con ripetizione: n = 2; k = 3 
D’(2,3) = 8 

aaa aab aba abb baa bab bba bbb 

3. Permutazioni: n = 3 
P (3) = 6 

abe acb bac bea cab eba 

4. Combinazioni semplici: n = 5; le = 3 
C(5,3) = 10 

abe abd abe aed ace adc bed bee bde ede 
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P.F. Tramontano 

ALGEBRA 

Teoria, esercizi e simulazioni 

Il software interviene a completare il testo 
laddove il calcolatore può arricchire in mo¬ 
do corretto la materia rendendola didatti¬ 
camente più valida e fornendo materiale 
per integrare i corsi istituzionali di Algebra 
Elementare. 

120 pagine 
Libro + disco 
versione software per: 

MS-DOS 
Cod. SD2921 

M20 • 

Cod. 292DOL 
C64 • 

Cod. 292DCM 
• a richiesta 

M. Beretta 

IDRAULICA 

Teoria, esercizi e simulazioni 


L. 50.000 

L. 50.000 

L 50 000 


P.F. Tramontano, B. Rinateti 

ANALISI 

Teoria, esercizi e simulazioni 

Il testo si rivolge in particolar modo a stu¬ 
denti delle scuole superiori ed il software 
allegato interviene a completare il testo 
laddove il calcolatore può arricchire in mo¬ 
do concreto la materia rendendola didatti¬ 
camente più valida e fornendo materiale 
per integrare i corsi istituzionali della di¬ 
sciplina 

224 pagine 
Libro + disco 
versione software per: 

MS-DOS 
Cod. SD293I 


HI L. 50.000 


M20 • 

Cod. 293DOL 


H L. 50.000 


C64 • 

Cod. 293DCM 
• a richiesta 

A. Conti, M Rosa - dot 

MECCANICA 


9 L. 50.000 


Con il binomio iibro-software vengono ri¬ 
presi alcuni aspetti classici della materia 
sia dal punto di vista teorico sia dal punto 
di vista applicativo attraverso le simulazio¬ 
ni consentite dal computer. 


Teoria, esercizi e simulazioni 

Il software allegato al testo affronta in 30 
problemi i principi base della cinematica e 
della dinamica, coprendo l'intero pro¬ 
gramma ministeriale di Fisica. 


144 pagine 

Libro + cassetta/disco 



300 pagine 

Libro + cassetta/disco 



versione software per: 



versione software per: 



C64 

Cod. SD296 


] L. 40.000 

C64 

Cod. SD290 

E 

| L. 40.000 

MS-DOS 

Cod. SD296I 

a 

L 50.000 

MS-DOS 

Cod. SD290I 

a 

1 

L. 50.000 

Apple • 

Cod. 296DAP 

a 

L. 50.000 

Apple • 

Cod. 290DAP 

a 

L. 50.000 

C64 • | 

Cod. 296DCM ! 

a 

L. 50.000 

C64 • 

Cod. 290DCM 

a 

L. 50.000 


• a richiesta 


• a richiesta 


A. Martinelli, M. Persico 

CHIMICA 

Teoria, esercizi e simulazioni 

La teoria molecolare e la stechiometria 
possono essere apprese attraverso eser¬ 
cizi e simulazioni mediante il software al¬ 
legato al libro. 


192 pagine 

Libro + cassetta/disco 
versione software per: 

C64 > 1 

Cod. SD295 L9BJ L. 40.000 


MS-DOS 
Cod. SD295I 



50.000 


E. Boni, S. Schacheri 

FISICA ATOMICA 

Teoria, esercizi e simulazioni 

Il dischetto insieme al libro si propone di 
presentare alcuni tra gli esperimenti più 
importanti che hanno portato alla "costru¬ 
zione”, oggi accettata, del modello di ato¬ 
mo 


168 pagine 

Libro + cassetta/disco 
versione software per: 

C64 —— 

Cod SD300 LgSj L. 40.000 


MS-DOS 
Cod. SD300I 


H L. 50.000 


C64 • 

Cod. 295DCM 
• a richiesta 


S L. 50 000 


Apple e 
Cod. 300DAP 


3 L. 50.000 


L. Fronzoni, S. Faetii 

COMPORTAMENTO DEI SISTEMI 
COMPLESSI: ORDINE E CAOS 

Teoria, esercizi e simulazioni 

Il disco, o la cassetta, insieme al libro han¬ 
no come scopo la descrizione logico-mate¬ 
matica dei sistemi complessi. Vengono di¬ 
scussi alcuni modelli per la descrizione 
delle instabilità ecodinamiche, laser, dina¬ 
mica delie popolazioni. Ciascun sistema 
presentato viene analizzato secondo le li¬ 
nee seguenti: 

- descrizione 

- parametrizzazione 

-comportamento dinamico (simulazione 
attraverso il calcolatore) 

160 pagine 

Libro + cassetta/disco 
versione software per: 


C64 • 

Cod. 300DCM 
• a richiesta 

L. Danusso, S. Schacheri 

CAMPI ELETTRICI 
E MAGNETICI 

Teoria, esercizi e simulazioni 

Il calcolatore permette la simulazione di 
esperimenti altrimenti difficili da elabora¬ 
re e preparare nella pratica di laboratorio, 
per una migliore comprensione di questa 
materia. Gli argomenti trattati coprono il 
programma ministeriale di elettromagne¬ 
tismo per le medie superiori. 

96 pagine 

Libro + cassetta/disco 
versione software per: 

C64 r—-, 

Cod. SD301 L. 40.000 


K3 L. 50.000 


C64 

Cod. SD297 



40.000 


MS-DOS 
Cod. SD301I 



L 50.000 


MS-DOS 
Cod. SD297I 



L. 50.000 


Apple • 

Cod 301DAP 



L 50.000 


C64 • 

Cod. 297DCM 



L. 50.000 


C64 • 

Cod. 301DCM 



L. 50.000 


• a richiesta 


• a richiesta 




F. De Michele, M. Rosa-dot 

PROBABILITÀ STATISTICA 
E TERMODINAMICA 

Teoria, esercizi e simulazioni 

Il software allegato al testo dà luogo ad un 
binomio particolarmente indicato per la 
simulazione di casi edi esperimenti anche 
sofisticati di statistica e termodinamica. 

160 pagine • 

Libro + cassetta/disco 
versione software per: 

C64 

Cod. SD291 

MS-DOS 
Cod SD291I 

Apple • 

Cod. 291DAP 

C64 e 

Cod. 291DCM 
• a richiesta 


L. 40.000 


L. 50.000 

L. 50.000 

L. 50.000 


M. Fusilli, M. Poivani 

L’ELETTRONICA 
DEL COMPUTER 

Teoria, esercizi e simulazioni 

I concetti fondamentali dell'elettronica 
vengono illustrati con numerose simula¬ 
zioni proposte dal software allegato. 

Volume 1 - 236 pagine 
Volume 2 -168 pagine 
Libro + cassetta/disco 
versione software per: 

MS-DOS 

Cod. SD298I Voi. 1 UH L. 50.000 

Cod. SD299I Voi. 2 ■■ L. 50.000 

C64 • 

Cod. 298DCM Voi 1 
Cod. 299DCM Voi. 2 
• a richiesta 
in preparazione 


Marco Rosa-Clot 

LA FISICA CON IL COMPUTER: 
LA DINAMICA 


H L. 50.000 
L. 50.000 


P.P. Tramontano, B. Rinaldi 

GEOMETRIA ANALITICA 

Teoria, esercizi e simulazioni 

Lo studio della retta, della circonferenza, 
della parabola e delle coniche avviene at¬ 
traverso l’utilizzo del binomio dischetto- 
libro in una logica di completasimulazione 


Il libro ed il software ad esso collegato 
permettono ad uno studente di esplorare 
problemi anche non elementari di dinami¬ 
ca del punto, dei sistemi e del corpo rigido. 
Per ogni esercizio vi sono 5-6 simulazioni, 
con possibilità di costruirne altre, per af¬ 
frontare aspetti particolari del problema in 
esame. 


su video. 


Il software integra i corsi istituzionali rela¬ 
tivi a questa disciplina 

270 pagine 

Libro + cassetta/disco 



160 pagine 

Libro + disco 
versione software per: 


versione software per: 

Apple 

Cod.550A 

a 

L. 50.000 

MS-DOS 

Cod. SD294I 

L. 50.000 

C64 

Cod. SD225 

HU 

L. 40.000 

M20 • 

Cod. 294DOL 

H L. 50.000 

MS-DOS • 

Cod. SD225I 

a 

L. 50.000 

C64 • 

Cod. 294DCM 

L. 50.000 

C64 • 

Cod. 225DCM 

a 

L. 50.000 

• a richiesta 


» a richiesta 






C. Bacchetti. L. Accomazzi 

L’INSEGNANTE 
E IL CALCOLATORE 

Come e perché il computer nella 

scuola 

Cod. SF226 

176 pagine L. 16.500 

La trasformazione della scuola in seguito 
alla sua progressiva informatizzazione, lo 
sviluppo in parallelo di nuovi strumenti e 
tecniche di insegnamento, costituiscono il 
nucleo centrale di questo libro 

C. Tirittico, R Traversini 

IL CAI 

Un manuale per l’uso 
dell’elaboratore nella didattica 
Cod. SF227 

128 pagine L. 13.000 

Partendo da un'analisi dettagliata delle 
tecniche e degli strumenti per l’uso dell'e¬ 
laboratore nella didattica, il testo fornisce 
metodi e suggerimenti per affrontare la 
scelta e il dimensionamento di sistemi per 
l'insegnamento con l'ausilio del Computer 
(CAI). 

G. Mauri 

COMUNICARE 
CON IL COMPUTER: 

I LINGUAGGI 

Cod. 5D284 

148 pagine L. 15.000 

II libro fornisce un utile strumento che con¬ 
sente all’Insegnante la formulazione di i- 
potesi su cui operare in un contesto non 
rigidamente predeterminato. 


C. Poma 

IL COMPUTER 
NELLA SCUOLA MEDIA 

Cod. 508A 

200 pagine L. 18 000 

Un libro adatto a chiunque si sta avvicinan¬ 
do per la prima volta al mondo dell'infor¬ 
matica e della programmazione; indicato 
in modo specifico per studenti e insegnanti 
delle scuole medie. 

G. Tonfom 

LA COMUNICAZIONE 
CAMBIATA 

Cod. 538P 

116 pagine L. 10.500 

Un libro per scoprire, analizzare, poten¬ 
ziare le capacità comunicative di ciascuno 
grazie ai nuovi strumenti della tecnologia. 

L. Ortolani. 0. Cepelli 

VALUTAZIONE 
SCOLASTICA CON L’AIUTO 
DEL COMPUTER 

Cod. SD232 

136 pagine L. 34.000 

Supporto DISCO PER APPLE 

Elementi fondamentali della statistica su 
cui fondare un metodo di valutazione sco¬ 
lastica semplice ed oggettivamente valido. 

A.I.C.A. 

SOFTWARE DIDATTICO 

Cod. 804C 

408 pagine L. 45.000 

Di ogni programma sono descritti gli obiet¬ 
tivi, gli elementi caratterizzanti dal punto 
di vista didattico e degli allievi a cui si 
rivolge, gli strumenti informatici su cui fun¬ 
ziona e la documentazione di cui è corre¬ 
dato. 

E. Tonti 

DIDATTICA CON 
IL PERSONAL COMPUTER 

Cod. 400A 

156 pagine L. 24.000 

Un libro per docenti orientati all’impiego 
del calcolatore in classe, con programmi di 
immediato utilizzo e spiegazioni sulla loro 
funzione didattica. 




Finito di stampare 
nei Febbraio 1988 
presso Rotolilo Cotogno Milano 




“Ogni scatola di detersivo contiene una figurina che trova 
posto su un album che ne ospita 100. Quante scatole di 
detersivo si dovranno acquistare mediamente per comple¬ 
tare la raccolta?”. Ecco un bel problema di calcolo delle 
probabilità: un solo dato, una situazione abbastanza reali¬ 
stica, una soluzione non banale. 

È questo uno dei 17 problemi proposti e risolti nel secondo 
capitolo di questo libro; la maggior parte dei programmi 
associati agli esercizi realizzano delle simulazioni il cui 
scopo è quello di verificare “sul campo" la validità dei 
risultati ottenuti teoricamente. Ecco dunque il calcolatore 
gestire gli acquisti di una nutrita schiera di collezionisti e 
sfornare, in pochi minuti, risultati che richiederebbero inte¬ 
re giornate di simulazione “manuale" (a proposito, per un 
album di 100 figurine, si dovranno acquistare mediamente 
519 scatole di detersivo). 

Gli esercizi sono preceduti da un capitolo dedicato alla 
teoria del calcolo delle probabilità; di statistica si parla, a 
livello molto elementare, nel terzo capitolo. Chi è interessa¬ 
to ad approfondire questa parte troverà, nella stessa colla¬ 
na, il volume "Probabilità, Statistica e Termodinamica" 
ulteriore materiale e il relativo software. 
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